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RP303 


NOTE ON CONTRACTION ae ICIENTS OF JETS OF 


By Edgar Buckingham 


ABSTRACT 


1. A jet of liquid from a sharp-edged orifice contracts more than a jet of gas 
‘because the gas expands after leaving the orifice; and theoretical expressions for 
the ratio of the two contraction coefficients have hitherto been deduced from the 
assumption that the expansion is the same in all directions. In the present 
paper this is replaced by the assumption that the force exerted on the upstream 
face of an orifice plate on the end of a pipe depends only on the static pressure 
and the momentum of the approaching stream. 
| 2. The resulting expression is compared with the experiments of Witte on 
| water, of Bachmann on air, and recent unpublished experiments on natural gas. 
The conditions presupposed by the theory are strictly satisfied in Bachmann’s 
experiments but in the others they are satisfied only when the ratio, 8, of orifice 
to pipe diameter is small. 
3. The conclusions are as follows: 
(2) Bachmann’s experiments on air at low differentials agree with Witte’s on 
water within 0.5 per cent. 
; (b) The measurements on natural gas agree with Bachmann’s on air within 
5 0.5 per cent. 
c) As far as the critical pressure ratio, the new formula represents the results 
on natural gas within 0.5 per cent for 8=0.27; within 0.3 per cent for 8=0.38; and 
iwithin 0.8 per cent for 8=0.5. For 8=0.7 the departure increases to nearly 4 
B per cent. 


CONTENTS 


[. Introduction 

Il. Development of a new formula 

Ran gc oe ee, nant nid eis a ire 
Notation 

Remarks 

Basis for a new treatment of the problem _____-___--.------- 
Isentropic flow from (pi, p1) to po 

Flow at constant density from (py, p1) to ps 

. Application of the fundamental assumption 

Ill, Examination of the experimental data 

lV. Comparison with the new formula 

¥. CO cdot dnnidmamede nA eladnatthieh oe aedbbeins 


I. INTRODUCTION 


NOOR Whe 


The pressure in a jet of fluid issuing from a square-edged orifice in 
thin plate does not fall to equality with the outside pressure until 
the filaments have become parallel. Between the plane of the orifice 
and the vena contracta there is a radial pressure gradient outward, 
and if the fluid is a gas it expands transversely as well as longi- 
tudinally, whereas in the case of a liquid there is no expansion. The 
cross section of the vena contracta is therefore larger with a gas 
ithan with a liquid. 

The ratio of the two sectional areas would be given by a complete 
olution of the equations of motion, but until the mathematical 
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difficulties of the problem have been overcome, any attempt to com- 
pute the ratio a priori must rest on simplifying, approximate assump- 
tions. Two such computations have rec ently been made,!? both 
based on the assumption that the gas in the jet expands uniformly j in 
all directions; but while the formulas obtained agree well with the 
experimental ‘data with which they have been compared, the funda- 
mental assumption has been criticized * as very improbable under the 
conditions of acceleration and distortion prevailing in the jet. | 
have therefore thought it worth while to start from a different initia] 
hypothesis and to compare the result obtained with certain additional 
experimental data which have not yet been published and so have 
not been available to previous writers. 


II. DEVELOPMENT OF A NEW FORMULA 
1. CONDITIONS 


The orifice is to be a cylindrical hole in the center of a thin, flat 
diaphragm across the end of a straight pipe. The outside space into 
which the discharge takes place is to be so large that the static pres- 
sure in it is sensibly uniform everywhere outside the jet from the 
orifice; and the jet will be regarded as having a definite boundary, as 
far as the vena contracta, the frictional drag at the boundary being 
negligible. We consider only steady states of flow. 

The fluid shall be some gas, such as air, which may be treated as 
subject to the ordinary ideal gas equations over the short range of 
pressure and density involved in flow through an orifice at jet speeds 
below the speed of sound in the jet. 

The Reynolds number shall be high enough that the motion is no 
longer appreciably affected by viscosity, and that the velocity and 
static pressure are sensibly uniform all over the section at the vena 
contracta, the pressure there being the same in the jet as in the 


surrounding space. 
2. NOTATION 


The following notation will be adopted, all values being expressed in 
terms of a set of normal units, such as the cgs or the “British absolute” 
system. 

A=the area of cross section of the pipe; 
mA=the area of the orifice, m being the area ratio and ~m=8 the 

diameter ratio; 

p: =the static pressure in the pipe ahead of the orifice; 

p: = the density of the gas at that place; 

S, = the linear speed of flow along the pipe; 

p2= the pressure in the outside space; 

p2=the density of the gas after isentropic expansion from (p,,) 

tO Po} 
S;=the speed at the vena contracta; 

a~ the contraction coefficient, so that the area of the vena contracta 

is uamA; 
M =the mass flow; evidently 


M = pamAp,.S, (1 





1R. Witte, Techn. Mech. u. Thermodyn., 1, p. 118; March, 1930. 
2G. Ruppel, Techn. Mech. u. Thermodyn., 1, p. 151; April, 1930. 
3A. Busemann, Techn. Mech. u. Thermodyn., 1, p. 338; September, 1930. 
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the outside pressure that would be required to produce this same 
mass flow, from the same initial pressure and density, if the 
density were kept constant at the value p,, which could, in 
principle, be accomplished by withdrawing heat at a suitable 
rate instead of allowing the flow to be adiabatic; 

the speed at the vena contracta in this case; and 

the contraction coefficient. By the definition of S; and by (1) 
we have 


umAp,S; ~ M pen PamAprS, (2) 


C,/C, =the specific heat ratio of the gas; 

the force exerted by the gas inside the pipe tending to push the 
orifice plate downstream off the end of the pipe. To simplify 
the appearance of the equations the following abbreviations 
will also be employed 


= = 


a ny a 


3. REMARKS 


If the density does not change during a fall of pressure, compress- 
ibility does not come into play, so that its value is immaterial and the 
shape of the jet depends only on the value of the Reynolds number 

2, »in which D is the diameter of the orifice and 7 the viscosity of 
the Guid. 

Since R depends on p as well as 7, if its values are all so high that 
changing from one to another does not affect the shape of the jet, it 
follows that the shape is unaffected, not only by a change of viscosity 
caused by substituting one fluid for another, but also by a simultaneous 
change from one constant density to another. Hence at high values 
of R the contraction coefficient will be the same in the supposed case 
of discharge of a gas at the constant density p; as it would be in the 
discharge of a liquid of that or of any other constant density; and the 
value of 1 may therefore be determined by experiments on water at 
hich values of R, 


4. BASIS FOR A NEW TREATMENT OF THE PROBLEM 


We assume that with given initial conditions (p;, p:), and with a 
given mass flow 1/—which implies a particular linear speed S, along 
the pipe—the force Fy exerted by the gas on the upstream face of the 
orifice plate is the same when the subsequent discharge through the 
orifice takes place isentropically, as if it went on with artificial cooling 
at constant density. This enables us to dispense with the assumption 
of uniform expansion beyond the plane of the orifice, and with several 
auxiliary approximations of which the importance is difficult to 
estimate. 

The approximation to reality will be closer for small than for large 

values of (p,—py)/p, and, probably, for small than for large values of 
m; and it is hardly to be expected that the results obtained from the 
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for pressure ratios near the critical. 


assumption will agree well with experiment for large values of m o; 


{ Vol. ¢ 


In so far as the assumption is valid, it permits of applying the 


momentum principle, and we proceed as follows: 


5. ISENTROPIC FLOW FROM (p;, ;)) TO pz 


The equations of energy and continuity are 


S—S,= 2y ; P( = y7 )= 7B 
y-1 p j 
1y 


2 |= 


> 


and 


1 
M= ApS; _ baMA pS = MamApyy Y S2 
whence 


S.= / mare 
2 i 5 
| py (1—pemtys) 


and 


1 + ee ae 
M= pamApyy 7 rae? Pie = 
pry 7 (G — p*m*y 7) 


The time rate of increase of momentum is 


9 
M(S,— 8,) =2HsMAPit 
1+pamy7 
and the force acting is 
F= A(p,— mp2) — Fy 
whence, equating (8) with (9) and dividing by Ap,, we have 


1 P Fy << 2uame 
— my i” emoree | 
' 1+pamy> 
6. FLOW AT CONSTANT DENSITY FROM (p,;, 1) TO pz 


The equations of energy and continuity are 


S2—S2=2"! (1-w) 
Pi 


and 
M= ApS, rae umAp,S; 


g =f 270 —w) 
“SV pi (1 — pm?) 


in 2p —w)_ 
a pi(1 — ym’) 


whence 


and 





(10) 


(11) 
(12) 


(13) 


(14) 









| Bu 


Sal 


» W 
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| The rate of increase of momentum is 


2umAp,(1— 
M(S;—-S,)= “EmAP oo) 
, and the force acting is 
F=A(p.—mps)— Fr 
| whence by (15) 
— may _ fo 2Hm (1 — w) 
[~~ Ap, 1+nm 


7. APPLICATION OF THE FUNDAMENTAL ASSUMPTION 


The assumption is that if equations (7) and (14) give the same 
value of M, the value of F, will be the same in (10) as in (17). 
Accordingly, we may eliminate the unknown force F, by subtracting 
(17) from (10) and after dividing by m we have 


2Hab —_ 2u(1—w) 


it cot 1+pm 


i (18) 
1+ pamy 7 


Upon equating (7) and (14), removing the common factors, and 
squaring, we have 
1 
2 ai Qap aw 
#0") — _Be yr re (1 9) 
—~Bm Senos 
1— pa m*y 7 


and after eliminating w between (18) and (19) we have a quadratic 
equation of which the solution may be expressed in the form 


i mae 
ree 1-41-49 (20) 


where 


1 
A=y7| n+ wem?—1)$— pem’y7 (1—y) | 
B= po 
Q=1- 


y 
wt 
= Gy(1-y7) 
Pa 


y 





The result embodied in equations (20) and (21) is to be compared 
with experiment and we have next to examine the available data. 


III. EXAMINATION OF THE EXPERIMENTAL DATA 


In the absence of accurate measurements of the diameter of the 
vena contracta, we have recourse to the values of discharge coefficients 
obtained from measurements of mass flow. 

The “adiabatic” discharge coefficient C, for a gas is defined by 


M= CMe (22) 
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where M is the observed mass flow and M,., is the “‘theoreticg]” 
mass flow found from equation (7) by ignoring the phenomena of cop. 
traction and setting u4z=1. Since C, is not constant, we may write 


Ca=CYa (23) 


where C is the limiting value which C, approaches as p2.= py, and the 

‘“‘expansion factor”’ Y, is a function of p/p, which reduces to unity 

when p.= pi. ' 
The discharge coefficient C, for a liquid is defined by 


M=C\M.. (24) 


in which M,.; is the “‘ theoretical” mass flow found by setting u=1 in 
equation (14). But if we set wa=u=1 and ps=po, equation (7) 
approaches (14) as a limiting form when p.=p,. Hence OC, and (' 
are identical, and instead of attempting to find the value of C by 
extrapolating to very small pressure differences from measurements 
of C,, we may equally well find it from experiments on a liquid. 
The stipulation that the Reynolds number shall be so high that its 

precise value is no longer important, implies that the dissipative 
resistances due to viscosity have become insignificant. The contrac- 
tion coefficients Ha and w» may therefore be regarded as sensibly iden- 
tical with the discharge coefficients C, and C, and for the purpose oj 
comparison with equation (20) we may set 


Ma Ca _ Y (25) 
aa. wa 

The experimental data suitable for a direct test of the formula are 
very few; for it was postulated that the discharge should take place 
into a space that was large enough for the static pressure in it to be 
uniform, and the only satisfactory experiments of this sort that hay 
come to my attention are those of Bachmann‘ on the discharge 0! 
air into the atmosphere through an orifice on the end of a pipe. 

The internal diameter of the pipe was 8.25 cm and the diameter of 
the orifice was 2.003 cm, so that the diameter ratio was 8 = 0.2428 and 
the area ratio was 8?=m=0.0590. Of the 18 values of OC, obtained 
at speeds below the critical, the one that was measured at the lowes! 
rate of discharge and is therefore the least reliable is far out of line 
with the others. The remaining 17 are represented without system- 
atic departures by the equation 


C,=0.5989 (1 +0.2431 PB) (26 
2 


the greatest departure being less than 0.5 per cent and the meal 
departure less than 0.2 per cent. 
The important thing for our present purposes is the expansioi 
iactor 
, oe PTD? ™ 
Y,=1+0.2431 oe (27) 
2 


4H. Bachmann, ‘ ‘nelees zur Messung von Luftmengen,’’ Dissetaien, Darmstadt; 1911, 





ee | 


a 
a 


i 
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ral F but it is interesting to compare the value of C with values obtained 
‘O0- B som orifices installed in pipes, with p. taken as the minimum pressure 
nite & 


’ observed on the downstream side of the orifice. 

' Witte’s ° experiments on water flowing through perfectly square- 

» edged orifices installed in smooth pipes, with the pressures taken at 

' the face of the plate, gave values of C which, for diameter ratios 

| between 6=0.2 and B=0.75, or 0.04<.m< 0.56, are represented with- 
out systematic error by the equation 


ee = = K=0.5983 + 0.395m? (28) 
vyl—m 












The maximum departure is about 0.25 per cent and the mean depar- 
ture for the 17 orifices tested is less than 0.1 per cent, so that the repre- 
sentation by equation (28) is very good. 

For Bachmann’s area ratio, (28) gives C=0.5986, but this must be 
corrected to obtain the corresponding value which would have been 
found if the pressure difference had been measured between the 
upstream and downstream minima instead of between the two faces 
of the plate. Experiments on air which will be published in the near 
future show that this correction factor is about 0.9983, so that we have 
from Witte’s experiments on water C=0.5976 as compared with the 
value 0.5989 from Bachmann’s work—an excellent agreement. 

The equation adopted for representing the results of the experi- 
ments on compressed air, described in Bureau of Standards Research 
Paper No, 49,° namely 
































C=0.5970+0.09 m? (29) 






gives C=0.5973 for Bachmann’s value of m. 
Next to Bachmann’s data, the best at my disposal were obtained 
from experiments on natural gas, performed at Los Angeles in 1929 
under the direction of H. S. Bean of the Bureau of Standards, as part 
of an investigation being carried on by the Natural Gas Division of 
the American Gas Association, in cooperation with the Bureau of 
Standards. The circumstances were more favorable as regards steadi- 
ness of flow than in the experiments on compressed air just referred 
to, so that the precision was higher and the results are preferable. 
The orifices were installed in straight pipes, as in Witte’s experi- 
ments on water, so that the discharge space was restricted and the 
static pressure in it not entirely uniform. It is to be presumed, 
however, that if the area ratio m is small, and if the minimum pressure, 
observed at the wall of the pipe on the discharge side is taken as the 
value of p., the arrangement will correspond fairly closely with that 
presupposed in developing equation (20); and the agreement noted 
_ above between Witte’s and Bachmann’s results lends support to this 
presumption. 

New commercial steel pipes of 4, 8, and 16 inches nominal diameter 
were employed, and a honeycomb of smaller pipes was placed in each, 
from 10 to 15 pipe diameters ahead of the orifice. The orifices (23 in 
number) had diameter ratios ranging from 0.124 to 0.869, but the 


f 





5 See footnote 1, p. 766. 
‘B.S. Jour. Research, 2, p. 561; March, 1929. 
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data to be considered here are only those obtained with the 14 orifices 
for which the diameter ratios were between 0.3 and 0.75, or 0,09 
<m<0.56. For, as is the case with Witte’s values of C for water, 
the variations of C can be more accurately represented by simple 
empirical equations within these limits than for 8<0.2 or >0.75. 

The results are expressed in the form 


y= 
3 Los 


C,=f(z) 
where x= (p,—p2)/p, and C, is defined by the equation 
M= Cima] 2p: (ps — Pa) (31 
1—m 


M being the actual mass flow. The coefficient C,; is connected with 
C, by the relation 








0; — OZ,én (32) 
where 
2 ca! 
= yw pat ed ) 
iy—t) i-9) 
and 
peer 
23> 2 
1— my? 





Since 2,+2,.+1 when p.=p, or y=1, C, and C, have the same 
limiting value C and we may write 


C,=CY, (34 


in which Yj, like Y,, is a function of z which reduces to unity when 
z=0. By equations (32) and (34) we have 


C, Y, 
0,=CY,=—1 =c — 
2122 2122 
whence 
— 58 
* 2122 


and if Y, is known, the value of Y, for use in equation (25) may be 
found from it. 

When the values of C, for each of the 14 orifices are plotted against 
the values of z at which they were measured, the points lie along 8 
straight line, within the experimental errors. The slope of the line 
varies with m, and the average result for ail the orifices 1s represente¢ 
by the equation 

Y, =1— (0.32 + 0.3m?)z (37 


In order to use this result in connection with Bachmann’s values 0! 
Y, for air, an allowance must be made for the difference betwee! 
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4 y=1.40 for air and y=1.284, computed from the composition of the 
| Los Angeles gas. For this purpose I have assumed the relation 


cs?) 


| proposed by E. S. Smith,’ which appears to be accurately true for low 
| values of z, and to be a good approximation even as far as the critical 
pressure ratio. 

On the basis of Smith’s relation, equation (37) is to be regarded as 
' a special form of the more general equation 


Y,=1—(0.41+0.39m?) : (38) 


and for the case of air, with y=1.40, this reduces to 
Y;=1— (0.293 + 0.28m?) x (39) 


For the purpose of comparison with Bachmann’s data on air and 
with equation (20), equation (39) will be assumed to represent the 
results that would have been obtained at Los Angeles if the experi- 
ments hed been made with air instead of natural gas. 


IV. COMPARISON WITH THE NEW FORMULA 


Since the computations required are tedious, they have been carried 
out only for a few simple values of m? and of p,/p;, and the experi- 
mental values of Y, have not been treated individually but taken 
from the empirical equations which represent them. 

Bachmann’s results have been taken from equation (27) or 


2 


Y,~=1+0.2431 - i (40) 


which represents them within the experimental errors. In accordance 
with (36) and (39), the equation 


Y fs. 1 — (0.293 + 0.28m?) x 





(41) 


2122 


we ng taken to represent the results of the experiments at Los 
Angeles. 

| The values of »=C required for substitution in equations (20, 21) 
were computed from equation (28), which represents Witte’s experi- 
' Iental results, and then corrected to what they would have been if 
the pressures had been observed at the upstream and downstream 
' minima, by means of factors determined from the experiments on air 
_ already referred to. 





‘The proposal was made in a paper presented to the World Engineering Congress, Tokio, 1929, and 


| ‘subsequently to the Annual Meeting of the American Society of Mechanical Engineers, New York, De- 


cember, 1929. The introduction of the “expansion factor” Y is also due to Mr. Smith. 
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In addition to computing values of y,./u from (20), I have also made 
corresponding computations by means of the equation 


9 


Me oN ae cit a 


m + —m 2 


deduced by Witte * from the assumption of uniform expansion, 
The results are summarized in the following table. 




















| se 
| | ‘ oo | Y. Ma Ma | 
. ; 2 | Bac hms ann | Los Angeles “ mn | 
| “s v¥m=B # | pr equation | equation | equation | equation | 
| (40) | (41) (20) (42) 
. oat 2 3 ‘ 5 6 7 8 
| 
.§ 1. 028 1. 030 1 
0. 005 0. 266 0. 5970 5: 4 sor ie : 
. \ , ene ° ° ( v4 
(. 00348) a ae: 6 1 167 1 170 i 
eet 1. 222 1. 221 1 
| 
, 9 1. 028 1. 030 1 
.8 1. 064 1. 067 1 
- 02 376 | . 5975 x 1. 109 1.112 1 
s(n 1. 170 1. 170 1 
| | o & 1.225 | 1. 222 1 
| | 
| .9 1.030 | 1.030 1. 
| 8 1. 067 1. 067 1. 
. 06 . 495 . 6003 a 1.115 1.113 ‘4 
.6 | 1.177 1.172 | a 
. 53 | 1. 233 1.223 | 1. 
| 1 
9 1. 042 1. 033 bet hi 
8 | 1. 091 Ree 
. 24 . 700 6038 | ~ fa 1. 149 + @saa 
| oO Bisiil 1. 220 ie) essa. Sd. 
: yp oe | 4,282 | 1.235 Sete 
| | | 











Columns 1 and 2 give the values of m? for which the — 
were carried out, and the corresponding values of the diameter 
ratio 8. The figures in parentheses refer to Bachmann’s orifice 
they do not differ enough from the figures above them to affect th 
values in columns 6 to 8. 

Column 3 gives the values assumed for yu in equation (20) and col- 
umn 4 shows the values used for the pressure ratio p./p,. The 
remaining columns of the table are sufficiently explained by thei 
headings. 


V. CONCLUSIONS 


1. Comparison of columns 5 and 6 of the table shows that the 
empirical equation adopted for representing the average results of the 
Los Angeles experiments on natural gas, when reduced by means 0! 
Smith’ s assumption that C,=f(z/y), agrees with the results of Bach- 
mann’s measurements on air within 0.5 per cent, the systemaii 
difference being within the precision of the Los Angeles experiments. 





§ See footnote 1, p. 766. 


Nace 


Contraction Coefficients of Jeis of Gas 


The assumption proposed by Hodgson,* namely 
e D> ] 


Ned Lo) | 


while not so simple as Smith’s relation seems better justified. If the 
reduction of the Los Angeles results from y=1.284 to y=1.400 is 
made in accordance with this equation, the greatest difference be- 
tween columns 5 and 6 is decreased from 0.006 to 0.002. 

2. Comparison of columns 6 and 7 shows that the new formula 
agrees with the Los Angeles results to within 0.5 per cent up to 
6=0.376, and to within 0.8 per cent at 8=0.495. For 6=0.7 the 
agreeme nt is poor throughout, as might be expected; for with so 
large a diameter ratio the condition of disc harge into a pipe is extremely 
different from the free discharge postulated in deducing the formula. 

3. Comparison of columns 6 and 8 shows that Witte’s formula 
(42) represents the facts, as far as they are known, nearly as well as 
the new formula, for pressure ratios down to about 0.7, but that it is 
inferior for still lower values of p2/p,. 

Wa ASHINGTON i, February 16, 1931. 


* Discussion tom apers on Fluid Flow seiamated at the annual athe 4 of the A. M. E. .» New York, 
December, 1929. 
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THE EFFECT OF SMALL VARIATIONS IN PITCH UPON 
THE INDUCTANCE OF A STANDARD SOLENOID 


By Chester Snow 


ABSTRACT 


\ method is described by which the variations in pitch of a single-layer solenoid 
may be compared with those of the precision lathe screw used to wind it. Assum- 
ag that the pitch of the latter is constant, but not known with precision, it is 
shown how the precision measurements of the length of the windings may be 
combined with the former measurements to obtain the variation in pitch of the 
windings relative to an ideal helix whose beginning and end points coincide with 
the actual one. A mathematical formula for computing the correction 6L to the 
inductance due to this variation in pitch is then developed. This formula 
requires the graphical integration of a certain function times the deformation 
and practical methods of evaluating it are described. Application is made to 
observations on two standard solenoids and shows that the correction 6L to the 
inductance L may amount to 4 or 5 parts in 100,000 even on the most carefully 
constructed solenoids. In absolute measurements the correction is, therefore, by 
no means to be neglected, although it seems hitherto not to have been taken 
into account in determinations of the ohm. 


CONTENTS 


. Introduction 
. Formula for 6L_--- 
. Practical method of evaluating the integral for 6L 
1. First method--~_~-_- 
2. Second method 
’. Numerical examples 
1. First example 
2. Second example 
. Summary : ox x 
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I. INTRODUCTION 


In computing the inductance of a single-layer solenoid, it is at first 
assumed that the turns of wire are uniformly spaced. Measurements, 
however, show that this is never the case, no matter how carefully 
they are constructed. It has been possible by several different 
methods to measure the very small distance 6(y), by which a turn at 
mean distance y from one end is displaced from its ideal position. 
The ideal position of the nth turn is the position of the nth turn of a 
uniform solenoid of the same total number of turns, N, whose initial 
and final turns are at the mean position of the initial and final turns 
of the actual solenoid. The constant pitch of the windings which is 
used in the formula is defined as the quotient of the measured length 
of the solenoid by the total number of turns. To the value L of the 
inductance, computed by any formula which assumes a uniform pitch, 
we must add a small correction 6Z, which represents the increase that 
would be made in the inductance of the ideal solenoid when its turns 
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are displaced by the amounts 6(y) which deform it into the actual 
solenoid. These displacements 5(y) of each circular turn of wire are 
only in the direction of the axis of the solenoid so far as concerns the 
present problems. Variations of diameter of the solenoid, as meas- 
ured over the windings or between turns, which would be due to the 
elastic deformation of the wires or of the cylindrical form, or of both, 
have been measured and their influence gives rise to other corrections 
which do not concern us here. However, for the sake of other possible 
applications, we formulate the first order correction 5L caused by 
any infinitesmal deformation of a circuit. 


II. FORMULA FOR 56L 


Let LZ, be the inductance of any sort of winding and let L,+51 
be the inductance when each element of the wire is given a small dis- 
placement which is represented by the vector 5, and which varies 
from point to point on the wire. Let 7, and 7',+6T be the electro- 
kinetic energy in the two cases, the current being one electromagnetic 
unit ineach case. Let q, be the geometric coordinates which determine 
the configuration of the wire. They are infinite in number and the 
wire will be treated as a linear circuit. The electromagnetic force 
(in the generalized a sense) tending to increase any such 

re) 


geometric coordinate is ay and the work 6W, which is done by 


t 
all electromagnetic forces during the displacement 6, is therefore 


> 5q:= 6T where 6q; is the 7 component of the displacement 4 in 
i 


the Lagrangian sense, and where the current is constant during the 
displacement. Therefore 67’=6W (with constant current unity). 


The electrokinetic energy is always given by T=5LP=5L (since 


the current J=1) so that 67’- SL = 6W, or 
6L=26W (1 


That is, the increase, 6, in the self-inductance of any circuit which is 
produced by giving it any infinitesimal deformation is equal to twice 
the work 6W, which is done by its own electromagnetic forces during 
the deformation with constant current unity. If the vector H is the 
magnetic field at a point where the element of wire is the vector ds 
(treating the wire as a linear element), which field is produced by 
unit current in the wire in its original configuration (omitting the 
contribution of the element itself), then the electromagnetic force 
upon ds is the vector product [ds x H] and the scalar product of this 
by the vector displacement ([ds x H]- 6) = ({H x 6]-ds) is the work done 
upon the element so that the total work done is 


bW= JS ((ds x H]-6) = S ((HX 6] - ds) (2) 


integrated with respect to ds, along the length of the wire in its initial 
or unvaried configuration. Combining (1) and (2) gives 


6L=2/S ((H x 6]-ds) = —2f ((H x ds] - 5) (3) 
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The integral last written represents the total lines of force cut by the 
displacement. The penalty for idealizing the wire as a linear circuit 
(with zero cross section) is the same as that incurred in the use of 
point charges. Since the element’s own field does not move it, this 
singular part of the field must be omitted from H. This may be 
done when the wire is a closely wound coil by taking H as that due 
to the equivalent uniform distribution of current, such as the equiva- 
lent current sheet in the case of a single layer solenoid or helix. In 
extreme cases, such as that of a very thick wire, the formulation (3) 
as a line integral would be inadmissible and must be replaced by an 
analogous volume integral over the volume of the wire. 

For the case in hand we replace the single-layer helix by N equal 
and parallel circular (linear) turns whose centers are on the y-axis 
with the origin y=0 at one end. In the initial configuration the circles 
are equally spaced and the inductance is L,. It is then deformed by 
displacing each turn (center at y) a distance 6(y) in the positive 
y-direction where 6(y) is derived from the observations. It may be 
positive or negative, and is not only small compared to the distance 
between two adjacent turns, but also varies slightly from turn to 
turn. Instead of integrating along the wire, we may then safely take 


: dy (where 1/=length of windings so that }-=number of turns per 


centimeter) as the number of turns in the interval dy. Each of these 
turns is acted upon by a force —H,(a, y) 27a in the +y direction, 
where a is the radius of the circular turn and H,(a, y) is the component 
of magnetic field perpendicular to the axis, which would be produced 
by the equivalent current sheet around which a current of strength 
y per centimeter along y is circulating right-handedly (each circular 
turn carries unit current). The equation (3) then gives 


ON Ll 
37, = ) He y) 2ra 6(y)dy (4) 
0 


If r, 6, and y are cylindrical coordinates, the r-component of H is a 
function of r and y only; that is, H,(r, y) and H,(a, y) denotes its value 
on the sheet where r=a. It is continuous at r=a@. We may compute 
H,(r, y) by recalling that H is the curl of a vector potential A, whose 
cylindrical components A, and A, are zero, while the component Ag 
isa function of rand y only. H, is given by 


_OAd(r, y) 


H,(r, y)= oy 


(5) 


where 


aN (' Qn cos 6d 6 
Ag(r y)= dy’ 7 = ———= ; 
L Jo 0 +yr’?—2ar cos 6+a’?+ (y—y’)? 


so that by (5) 
aN(',,a(%___—s—cos 0d 8 a 
H,(r, y) = — 4 fev =f, \r?—2ar cos 6+a?+ (y—y’)* 


aN (* dhy' 0 ( cos 0d 0 
i f, y yl, r?—2ar cos 6+a?+ (y—y’)? 


49527°—31——2 
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Integrating with respect to y’ gives 


1 
H,(r, y= tl cos 0d 0 > 


dr -2ar cos 6+ a2 + (I—y) 
Te Ahead, edness, ts 
re Vr?—2ar cos 6+ a?+ y”) (7) 


Placing r=a gives 





QaN 
-H,(a, y=") le) - ex’) ( 
where 
uit oii 
in y \and p?=_ (l-y (9 
1+(¥.) : (| ) 
and 
rei (u)=— mu F pe. =| KW) a «) |- uK(u) (10 
w ¥l—yp* sin’ 0 a 


where K and E are the complete elliptic integrals of the first and 
second kind of modulus yn. The formula (4), therefore, becomes bj 
the use of (8) 


82’aN? (* , ( 
Pa . 6(y)Le(u) — er(u’)ldy (11 
where yu and y’ are functions of y given by (9), and of course 6(y) is 


the observed function of y. We may replace the latter by a function 
Fi(y) defined by 


fily) = 6(y) — 6—-y) (12 
The equation (11) then transforms into 
8 
6L=- — “f fueuluddy (13 


It is next convenient to introduce a new variable of integration by 
the transformation 
=? : (14 


"74 i 


If we define f(A) as f;(2ad) and ¥() by 
1 


V(\)=r7¢,(u) where p?= i+? (15) 
then (13) finally becomes 
N2 (>= 
L=4x( 5 y{ f(A)W(d)Ed (16) 
m 0 
where 
Xn = =the maximum value of X (17) 
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The deformation of the helix is represented in (16) by the small 
> observed function f(A) which may be plotted as a function of \ from 
> \=0 to A=Am. The function ¥(d) defined by (10) and (15) has 
' been computed by the use of tables of the elliptic integrals and 
F plotted in Figure 1 from A=0.07 to A=5.0. Since it becomes 
logarithmically infinite at A\=0, the equation (16) needs some re- 
| shaping in order to become practicable for the evaluation of 6L by 
| craphical integration. 


éd 





4 J 


Ficure 1.—The function y(a) 


_ Ill. PRACTICAL METHOD OF EVALUATING THE INTEGRAL 
FOR 6L 


1. FIRST METHOD 


Evidently a positive value « may always be so chosen that the 
curve of observation f(A) may be represented by 


S(\)=f(o)+f’(o). » in the range o=A=e (18) 


so that (16) may be written 


L an() | flo) f ‘Var +f’ f, “Aw(ndd) + f he fond) (19) 


| For the last integral the integrand is easily plotted and is finite, and 
_ may be evaluated by graphical integration. The first two integrals 
| are finite although the integrand becomes infinite in the first, but both 
| tay be evaluated analytically by means of two integral formulas 
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which are given in Bureau of Standards Scientific Paper No. 537 op 
page 449, which are 


(udu _ Vi—# 
foe VIT# (KW) EW) (2) 


— 5 feeyee_ Pole Ki.) — E(u)|— E(u) (21) 


3p 


Since 7¢;(u) =y(A), and since the dependence of » upon d becomes 


1 
v-T+s (2) 
it is evident that =. Add and — oe ==d). The left sides of (20) 
we wey l— pe 


and (21) become /y(A)ddX and fAy(A)dX, respectively. Taking 
the definite integrals from 4=0 to A=e (noting that for the partic. 
ular function considered {(0)=0), we find from (19) after placi ing 


1 
= 93) 
v1l+é ( 


the final formula for computing the increase 6Z in inductance due to 
the deformation. 


sL—ar(3-) (3[ 14225 e@-Ee))-22 |y0) 
+ f¥oosondan| (24 


Assuming that the straight line representation (18) of f(A) wil 
generally be possible for the range 0< A<e=0.1, this becomes 


L= ar(5-) | 0.011f’ (0) + f v NFO)da} (24a 


To compute 6L by this formula, it will be recalled that N is the total 


, ev? e-% 
number of turns each of radius a and that d is - and dA,» is - where the 


length of the solenoid is / cm and y is the distance in cm from one end 
to a variable point and the method of constructing and plotting the 
small function f(d) is as follows: The observed displacement in cm 
relative to the ideal helix of a turn at distance y from one end being 
the small function of y, 6(y), (which is reckoned as positive whet 
the displacement is in the direction of increasing y), we define /,() 

by fi(y)=45(y)—6—y) and then make the transformation y=2a\ 
representing the function f,(y)=f,;(2a\) by the new symbol /()). 
This may be plotted as a function of \ and the ordinates when multi- 
plied by the ordinates of the curve Figure 1 for ~(A) give the inte 
grand of the integral in (24) or (24a). If this integrand, y(A)f(A), be 
plotted as a function of i, the area under the curve can be found by 
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graphical integration and represents the integral desired in (24). 
Since \is dimensionless and 6 (y) and f(A) and /’(A) are given in cen- 
timeters, the change of inductance 6L given by formulas (16), (24) 
or (24a) will be measured in absolute c. g. s. electromagnetic units 
or centimeters. (1 cm=10~° henries.) 


2. SECOND METHOD 


The procedure outlined above is applicable when f(A) is any func- 
tion which is finite and has a finite slope at A=0. If f(0) is not zero, 
an improper integral is involved and this method of evaluating the 
integral analytically over a small initial range must be used as graph- 
ical integration could give no certainty in the case of a curve going 
to infinity at X=0. Since, however, the function {(\) always vanishes 
at \=0 and has a finite initial slope, the integrand ¥(A)f(A) is always 
finite (in fact zero) at A\=0O, and there is no practical difficulty in 
evaluating it graphically over the entire range. The practical diffi- 
culty is in plotting it mght down to the point \=0. This may be 
overcome by obtaining an asymptotic expansion for y¥(d) in the 
neighborhood of }=0. Such an expansion is given in Bureau of 
Standards Scientific Paper No, 537 on page 449, formula (56) by which 


3 , 4 3 2 
r#,(u)=[ 1 | a.) |In v1 pees 2- 4 (1 — p’) 


in which additive terms of the order of (1— y?) log - = have been 


1—- 
2 


r F 
inh this becomes 


neglected. Since 1—y?= 


vay=(14+9")ing-2-3™ (25) 


which may be used for values of \ less than 0.1 with an absolute error 
which is then always less than 0.14/n 10 =0.00023. 

If we use (25) for computing the integrand ¥(A)f(A) of equation 
(16) in the small range, say, 0<(\<0.1 then formula (16) may be 
readily evaluated graphically. 

_ In the first example below, the curve for f(A) is practically a straight 

line through the origin for 0< \<<0.1 and the formula (24a) has been 
et, convenient. In the second example, the second method 
is used. 

Before taking up these numerical examples, the method should be 
described by which the value of 5(y) was inferred from the observa- 
tions, As this method was somewhat indirect, new measurements 
are being considered by which 6(y) may be measured directly with 
relerence to a standard of length by means of a cathetometer. 

We define the ideal solenoid of constant pitch as one of the same 
number N of turns whose beginning and end coincide with the actual 
one. Taking the y-axis as the axis of the cylinder with the origin 
at one end, let 5(y) be the small displacement in the y-direction which 
deforms the element of the ideal winding, which is distant y from the 
end of the coil into its actual position. Evidently 5(y) is zero at both 
ends; that is 

(0) = 6(l)=0 (26) 
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Since the cylindrical form is of quartz, it may be assumed to have 
the same length / when mounted on the lathe for winding and com. 
parison with the screw as it was when finally set up for measure- 
ment of its length and diameter. 

The nearest substitute for the ideal coil in terms of which the actual 
one may be calibrated is the screw. The displacement A(y) which 
would deform an element of the screw thread into the position of the 
corresponding element of the actual winding was determined as a 
function of y by means of an optical lever on the carriage which made 
contact with the side of the wire. As the coil was slowly rotated and 
the lever carried forward on the screw, four observations were taken 
for each turn—4,000 readings in the case of the quartz cylinder with 
1,000 turns. The scale deflections, multiplied by the calibration 
constant of the lever, gave A(y). These measurements were taken 
immediately after winding, with the cylinder geared to the screw 
exactly as in the windings, so that the element of winding was com- 
pared with practically the same part of the screw that was used in 
winding it. Hence, if the wire had not rolled in going on to the 
cylinder its pitch should have possessed practically the same variations 
as the screw and no deflection of the optical lever should have been 
observed in moving from one end to the other, even though the screw 
were variable in pitch. Deflections A(y) were observed, however, 
and these were reasonably concordant with a similar set of observa- 
tions (which were taken on one coil) in which the optical lever was 
replaced by a micrometer microscope The value of A(y) was some- 
times positive and sometimes negative, but never more than 2 per 
cent of the pitch or distance between two turns, which was 0.1 cm 
in the case of the quartz cylinder of 1,000 turns, whose length was 
very approximately 100 cm. The displacements 6(y) and A(y) are 
reckoned as positive when they are in the direction of increasing y. 

If we let A,y denote the displacement (in the +y direction) which 
would deform the ideal helix into the screw, then 


5(y) = A(y) + A,(y) (27) 


Now if the screw is assumed to have a constant but unknown pitch 
at the time of comparison, then 


A,(y) =A+ By where A and B are constants (28) 


the latter depending upon the temperature. A change in the value 
of A merely shifts the screw as a rigid body and is irrelevant. Com- 
bining (27) and (28) gives 


5(y) = A(y) + A+ By (29) 


The constants A and B are determined by the two conditions (26) 
that 6(0)=6(/)=0. They are 


A=-—A(o) 
B=7{A(0)- AQ) 
so that (29) becomes 


ity) =a) —| S952 


y+A0) | 
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The observations are given in the form of a curve in which is 
plotted as a function of y from y=0 to y=/, the small observed func- 
tion A(y), which is the displacement of the winding relative to the 
screw. The equation (30) then gives the desired function 46(y) 
representing the displacement which deforms the ideal solenoid into 
the actual one. The curve for 5(y) may be plotted without computa- 
tion by drawing a straight line through the end points of the curve 
for 6(y). This is the straight line whose ordinates are represented 
by the bracket which is subtracted from A(y) in equation (30). 
Hence, the vertical distance of the curve for A(y) above this straight 
line represents 5(y). (See fig. 2.) 


Figure 2.—The observed variations in pitch of coil No. 1 


IV. NUMERICAL EXAMPLES 
1. FIRST EXAMPLE 


This is the case of a single layer helix of 1,000 turns of length 
/=100 em wound upon a quartz cylinder whose radius a is 14 cm, so 


‘ ) 
that N=1,000 and \»= oy 3.57 and formula (24a) becomes 


5L= aa 0.011 f’(0) + i 4(N)f(A)d | 


The observed variation in pitch is shown in Figure 2 by the curve 
which gives A(y) the displacement which transforms the screw thread 
into the actual winding. The curve in this figure for 5(y) was obtained 
graphically as explained since 6(y) is the height of the A(y) curve 
above the straight line connecting the ends of the latter. The 6(y) 
curve represents the displacement which would deform the ideal 
helix into the actual windings. From 6(y) the function f,(y)= 
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5(y) —-6(/—y) was constructed and plotted in Figure 3 as f(A) where 
y 
X a > ~* 
28 
y being taken from the curve of Figure 1. The area under the latter 
curve of Figure 3 was found by use of the integraph to be 


In Figure 3 is also shown the integrand ¥(A) f(A), the factor 


3.57 
i} v(A)f(A)dd = —4.7(10)~* with error < 1 per cent of itself, 
0.1 


The unit is here the centimeter of inductance (10~°), not area on the 
paper. The initial slope of the curve f(A) in Figure 3 is 


Ff’ (O) = —105(10-*) em 


| 


TR Tle SN RS 


z BIN (ian oa 
Figure 3.—The functions f(r) and ¥(A)f(A) for coil No. 1 


Hence formula (24) becomes 
_ 4007 
(3.57)? 


which must be subtracted from the value L, which is computed on the 
assumption of constant pitch. Since L, is approximately 68 x 10° cm, 
this correction is about 8 parts in 1,000,000. 


6L= {1.15+4.7} = —580 cm 


2. SECOND EXAMPLE 


This is the case of a single-layer helix of 270 turns with the same 
nominal pitch (0.1 cm) as in the preceding case, but wound upon a 
cylinder of pyrex glass, the mean diameter of the coil being 2a = 44.56 

7 
-—“"__( 605. In this case there were 

44.56 
two types of observations whose discrepancy served to indicate the 


cm, so that since 1=27 cm, },,= 
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probable total error of the method. The first was by means of the 
optical lever; the second was by means of readings taken with a 
micrometer microscope. The two sets of data are shown in Figure 4, 


No. 2 





~ 


} 
~f0A) OHO LEVER 


} 


| 
f hy Aaceo sper} 


Fiaure 5.—The functions 6(A) and f(A) for coil No. 2 

which gives Ay in the two cases. The corresponding 6(y)’s were 
constructed graphically and shown in the upper curves of Figure 5 
and the f(A)’s in the lower curves. In Figure 6 are shown the two 
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curves for ¥(A)f(\), these being constructed by the use of formula 
(25) for values of \ less than 0.1. The area under these two curves 
was found by the use of an integraph and gave 


0.605 


¥(A)f(A)dA= + 4.65 (10)~* cm for optical lever data 


= +4.47 (10)~* cm for microscope data 
Hence 
eee Al 


2 0.605 
ae) f ¥fd=1,160 cm for optical lever data 


= 1,120 cm for microscope data 
(The units here are in centimeters of inductance.) Since L, is approxi- 
mately 26(10)° cm, the difference between these two corrections js 
% \ 


> 
S 


50 


Figure 6.—The function ¥(d)f(d) for coil No. 2 


about 2 parts in 1,000,000, which gives an idea of the probable error 
in the correction. Giving the optical lever data the weight 2 and the 
microscope data the weight 1, the best value of the correction due to 
variation in pitch seems to be 6L= +1,150 cm (to be added to L,). 
This correction amounts to about 4.4 parts in 100,000. 


V. SUMMARY 


A general formula (16) is derived for computing the increase in self- 
inductance of the helix when its variation in pitch is observed. The 
practical form of evaluating it graphically is explained under (24) and 
(24a), and also in connection with equation (25). The directions for 
constructing the observational function are summarized after formula 
(24a). The application shows that this correction may amount to 
4 or 5 parts in 100,000, which shows that it must be taken account of 
in absolute measurements which use the coils as standards whose induc- 
tance may be computed from their geometric dimensions, 





VI. APPENDIX 


q 1, AN INDEPENDENT METHOD OF DERIVING THE EFFECT OF SMALL 
' VARIATIONS IN PITCH ON THE INDUCTANCE OF A STANDARD 
SOLENOID 


(The following method was derived by F. W. Grover during the editorial reading of the 
preceding paper) 


The inductance of a coil may be obtained by summing the self-inductances of 
the turns and the mutual inductances of all the pairs of turns of which the coil is 
composed. We will assume, as has been done in the preceding paper, that the 
coil is made up of parallel, coaxial, circular turns of wire; that is, the helicity will 
be neglected. Thus the effect of irregularities i in the winding pitch will be taken 
into account by supposing each turn to be displaced in the axial direction by an 
amount equal to the average of the displacements of all the elements of the turn 
from their ideal positions. 

It is evident that the self-inductance of an individual turn will not be altered 
by such a displacement. The mutual inductances of the various pairs of turns 
are, on the contrary, changed by the displacements of the wires from their ideal 
positions, and the algebraic sum of the changes of mutual inductance of all the 
pairs of turns of the coil will give the change of inductance of the coil, due to the 
assumed irregularities of winding. 

If two turns, whose axial coordinates are y and (y-+ Ay) in the ideal coil are given 
displacements 6; and 62, respectively, their mutual inductance suffers a change 
equal to 

dM 


6M=— (81-83) "Gy 


(1) 


: dM : : : 
Furthermore, di will be the same for all pairs of turns which are equally spaced 


in their ideal positions. Applying these considerations to the observed displace- 
ments 6, 52, 63,—.. 5, of all the turns of a coil of N turns, wound with an average 
pitch g, the inductance of the coil is changed by an amount 


sL=2eM=—2[ (GE) {(— a) +d) +++ + Gedy) 


(BY 1 (5; — 83) + (52— 54) + +++ + (Sn-2—5n)} 


dy 
dM e 
doses +( Se S tndaid: 
or, rearranging 


igh x a aM dM l 


7 ((dM wet 4 (a 
+ (é: “i ‘a dy ad 


i CAEP RRR b's - 


: : M 
in which (2 =). is the rate of change of the mutual inductance with respect to 
their axial pa al apart, of two turns separated by a distance pg. The final 


term in (3) is 
dM 
Ga—dn1)(4 ) , if N is even 
2 3 \dy/y 


2 


(Sn-1 — dnas,( SO acts if N is odd. 
7) 3 2 
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dh 

The calculation of 6L is thus seen to depend upon the calculation of (a) f 
for pairs of coaxial circular filaments, with axial separations g, 2g, 3g, etc. ai) 

dM 
the force between two filaments carrying unit current is -(<), the terms of 

( 
equations (2) and (3) are seen to represent work done by the electromagnetic 
forces during the deformation from the ideal winding, a result agreeing with 
equation (1) of the text. 

dM 


Except for turns very close together, the values of the (= dy ) can not be ex- 


pressed by a single series expression. However, the well-known elliptic integral 
formula of Maxwell, Elect. and Mag. volume 2, page 339, is universally appli- 
C: able, or more conve niently, numerical values may be inte rpolated from the tables 
of Nagaoka and Sakurai, Institute of Physical and Chemical Research, Table 2 3 
Tokyo, 1927. 

The calculation of 6L for short coils, or coils of only a few turns may readily 
be carried out by the use of formula (3) and such tables. For a coil of many 
turns, however, this process would become very laborious. In —_ cases the 


winding pitch is small compared with the dimensions of the coil, and (3 dy *) changes 
slowly with the spacing of the turns. Thus it is possible to transform the sun- 
mation of equation (3) into an integration over the length of the coil, and it will 
be shown that this leads to Doctor Snow’s equation (18). 

To prove this, we note that the displacements of the turns enter in (3) in the 
form 6(y) — 6(l— y) =fi(y), so that on summing along the coil, under the assumption 

: aa L es / 

that Ay=g is very small, each term of (3) takes the form = [6(y) —6(l—y)](M,- 
MI—,), where M, and M,—, are the mutual inductances of equal coaxial circles, 
separated by axial distances y and (l—y), respectively. 

Thus, if we remember that the winding density is 7 turns per cm. of axial 


length, equation (3) becomes 


Ee Nd 
6L= i, TA) (M,—M_r,)dy (4) 


It is, however, to be noted that the function 7¢,(u), of equation (15), is a factor 
in Maxwell’s elliptic integral formula for the mutual inductance of equal coaxial 
circles (loc. cit.), so that we may write 


M, _ 
“we A) 
and, therefore, 
N 
BL = SE LW) Wu) Idy 


The second integral in this equation may be transformed by changing the variable 
to y’=l—y, so that it takes the form fin ¥(A)dy. Thus (7) becomes finally 


SxaN? ¢ / 7 
6L= |, fiw (Ady (6 
which is exactly the equation (13) of the paper, since w¢;(u) =y(A). 


The alternative method is thus seen to lead to exactly the same correction for- 
mula for the inductance change due to variations in winding pitch as was found 
by Snow’s more general solution of the problem in the preceding paper. 


WASHINGTON, September 15, 1930. 


DEAT CONTENT VALUES FOR AQUEOUS SOLUTIONS OF 
THE CHLORIDES, NITRATES, AND HYDROXIDES OF 
HYDROGEN, LITHIUM, SODIUM, AND POTASSIUM AT 
18° C. 

By Frederick D. Rossini 


ABSTRACT 


The existing data on the heats of dilution of the chlorides, nitrates, and hydrox- 
ides of hydrogen, lithium, sodium, and potassium at various temperatures and 
co neentrations have been converted to one temperature, and a series of values 
the concentration range from infinite dilution to about 2 molal have been 
obtained. The extrapolation to infinite dilution has been made with the aid of 
» Debye-Hiickel theory of strong electrolytes. 
“The data, exhibited in tabular and graphic form, give values for ®,—#®,°, 


the relative apparent molal heat content of the solute; H-—i,° , the relative 


partial molal heat content of the solute, and H,-—i,°, the relative partial molal 
heat content of the H,O. The procedure employed in calculating the above 
quantities is given in detail. 


CONTENTS 


I. De finitions NEA Ree et gs 2 RIE SS ORE ee eh ee ee LTE ee cs 
II. 
III. 
IV. The ‘ 
’. The cahanladed meson Oa bi. 
VI. Discussion 


I. DEFINITIONS 


In this discussion the nomenclature of Lewis and Randall! will, in 
the main, be followed: 

m is the molality in moles of solute per 1,000 g H,O. 

The subscripts ! and 2 refer to the solvent and ‘solute, respectively. 
_ The superscript zero on any symbol refers to the given property at 
infinite dilution. 

H (printed in italic caps) is the symbol for the total heat content of 
a given solution. 

H, and FH, (printed in italic caps with bar) are the partial molal 
heat content * of the solvent and of the solute, respectively, and are 
defined by the following equations: 


BE 
Fi,= ($7 )n, 2,7 (1) 


-- 61 
Ha=(Gy, Jo? *, 


Lewis and Randall, Thermodynamics, Chap. VIII. McGraw-Hill Book Co., New York; 1923. 
, the partial molal symbols, the use of the italic cap letter indicates a molal property while the overhead 
if is the differential mark. 
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where m, and nz are the number of moles of the solvent and solute 
respectively. 

®, is the apparent molal heat content of the solute. 

MX-n,H,O designates a solution composed of 1 mole of th 
solute MX and n, moles of the solvent H,O. 

The heat content of any solution is given by the equation:* 


H=n,M, + nHe 
or by definition of ®,, oT 
H=mH2° +nPy, 
For the reaction 


mM X -n,H,O + (co = n,)H,0 =n.MX - co HO 


AH = [noHT2° + © H,°|—[noHs +m, + (@ —n)H,°] - n2(H2° — H)) 
+n,(H,° — A) (0 


or, using the apparent molal heat content function, for the same 
reaction, 


AH =[n2®,° + © Ff,°| — [nb +A? + (ao — n)Ay°| = No(#,° — P,) 
Combination of equations (6) and (7) gives 


N2(P, — &,°) = no(He— He”) + nm (Hi; — H,°) 


},—®,°= H,—H,° +=, — Fh’) 


®,—®,° will be called the relative apparent molal heat content’ 
of the solute; H:—H,° °, the relative partial molal heat content o/ 


the solute®, ond H,—H,° , the relative partial molal heat content 
of the solvent. 


II. INTRODUCTION 


Values for 6,—®,°, H2—H,°, and H,—H1° are needed in many 
thermodynamic calculations involving aqueous solutions. If value: 
are known for the relative apparent molal heat content of the solute 
and for the heat change in a given reaction at one concentration, ou! 
can calculate the heat change for the given reaction at any desired 
concentration. 

The partial molal heat content of the solute and of the H,O serve 8s 
a measure of the change of the respective partial molal free energie: 
with temperature, according to the thermodynamic formula 


dq ta _H-& (10 
d’ ae T? 


’ See Chap. IV of reference 1. 
‘ See footnote 3. 


5 Lewisfand Randall (see footnote 1, p. 791) use the symbol Za for this property, but in order to better kee 
in mind the physical significance of it, we shall use the complete designation 
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) One application of this formula is in the calculation of activity 


efficients from freezing point data, where the value for H.-H,’ is 
needed. 


x In order to evaluate ®,-®,°, H.-H,,°, and H,-H,°, there must be 


‘available data on heats of dilution, or heats of solution, extending 


‘to very small concentrations from which extrapolation to infinite 
dilution will introduce no appreciable uncertainty. Until recently, 
"the only extensive data on the heats of dilution of uni-univalent 


electrolytes in aqueous solution were those of Richards and Rowe,° 


| who measured heats of dilution from MX-25 H,O to MX-400 H,0O, 


and those of Pratt,” who covered about the same range of concentra- 


"tion but at many more temperatures. Without additional informa- 
' tion, the extrapolation from MX -400 H,O to MX- ~ H,O is, however, 
' over too large a gap to bridge with any measure of certainty. Later, 
' Richards and coworkers extended the data on HCl to 800 H,O and 
on NaOH to 3,200 H,O. Following 1927, there began to appear the 
© data of Lange * and coworkers, who measured the heats of dilution of 


salts at extremely small concentrations. 

However, the Richards data were obtained at 16° to 20° C., while 
the Lange measurements were made at 25° C. Inasmuch as the 
temperature coefficient of the heat of dilution is rather large, corre- 
lation of the two sets of data with those of Pratt would be impossible 
without sufficiently accurate data on heat capacities in the given 
solutions. These data are available from values which have been 
compiled by the present author,’ chiefly from the experimental data 
of Randall and Rossini,’ and Richards and coworkers." 

For extrapolating the data to infinite dilution, a new tool has been 
made available with the extension to heats of dilution of the Debye- 
Hiickel * theory of strong electrolytes. This extension “ involves 
simply taking the temperature coefficient of the original expression 
given by Debye and Hiickel. That the heat of dilution of strong 
electrolytes in the region of very small concentrations follows more or 
less closely the predictions of this theory was first shown by the work 
of Nernst and Orthmann,' and Lange and Messner," and later by 
Lange '° and others. 

A short derivation of the equations expressing the change of the 
apparent and partial molal heat content with concentration, as derived 
theoretically from the Debye-Hiickel theory, is given here. 

Consider the two systems A and B. A is an infinitely dilute 


' solution and B is a dilute solution of molality m. The difference in 


the free energy of 1 mole of solute in the two systems is 


AF=F,-F,° (11) 


* Richards and Rowe, J. Am. Chem. Soc., 43, p. 770; 1921. 
Pratt, J. Frank. Inst., 185, p. 663; 1918. 
* Lange, Fortschritte der Chemie, Physik, und physikalische Chemie, 19, No. 6, p. 1; 1928. 
These data will appear in B. 8S. Jour. Research, 1931 
" Randall and Rossini, J. Am. Chem. Soc., 51, p. 323; 1929. 
. Richards and Hall, Mair, Gucker, Dole, J. Am. Chem. Soc., 51, pp. 507, 712, 727, 731, 740, and 794; 
1929. Richards and Rowe, J. Am. Chem. Soc., 43, p. 770; 1921. 
. Debye and Hiickel, Physikalische Zeit., 24, p. 185; 1923. Debye, ibid., 25, p. 97; 1924. 
re Gross and Halpern, Physikalische Zeit., 26 p. 403; 1925. Bjerrum, Z. physik. Chem. 119, p. 156; 1926, 
‘toss, Monatshefte f. Chemie, 48, p. 243; 1927. Lewis and Randall, “Thermodynamics,” German trans- 
lation by Redlich, Berlin; 1927. 
Nernst and Orthmann, Sitzb. Preuss. Akad. Wiss., p. 51; 1926. 
Lange and Messner, Z. Elektrochem., 33, p. 431; 1927. 
* For a complete list of this work see Lange, Z. Elektrochem., 36, p. 772; 1930. 
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Following Debye and Hiickel, this free energy difference may }, 
separated into two parts: (1) The “classical” or that due to the free Jie 
energy change with concentration where the molecules of solute ar Hy < 
considered as particles of an “‘ideal’’ solution, and (2) the “electri. 
cal’’ or that due to the free energy change resulting from the electrica| HB , 
forces of the ions. Then 5 





























ee 


= s 


a 8. 5 a : =] 
F,- F,° r [F’, (class.) — B$  cotnesy) T [PF (elec.) —H5° (elec.)] (12 


“ER ; ; 
For very low concentrations, one may write 


F-Fe=rTl | - N® é/ «x y 2 a 
TT pe ra (7000 k} (2y,22)?m'? Us 


where the last term is that given by Debye and Hiickel for the chang 
in ‘“‘electrical”’ free energy for the very dilute range with the excep. 
tion that m has been substituted for the concentration. In this 
expression N is Avogadro’s number, e is the electronic charge, Dis & 
the dielectric constant of the medium, *& is Boltzmann’s constant, 7 
is the absolute temperature, z is the charge on the ith kind of ions 
whose number per molecule of solute is v, and m is the molality as 
already defined. 

By operating upon equation (13) according to the thermodynami 
relation expressed by equation (10) 


F,— F, A,- ~ H,° 7 
mes -)= a 


one obtains the expression 


ines ir, N* e* TdD - “ae 
4,— Hi’ =—h pets raa0k) (4 Da?) nad)? 4 


New data on the dielectric constant of water and its temperatur 
coefficient have recently been published.’ These have been re- 
viewed and compared with the existing data * by Lange and Robin- 
son.” Following the suggestion of these latter authors, the values 


dD : P ; 
for D and qT used in the present work were obtained by averaging & | 
the values of Drude, Wyman, and Drake, Pierce, and Dow. Thi 
gives the following values for 18° C. | 


ia A dD ai . 
D=81.15 qT" 0.37 a 
Then, for a uni-univalent electrolyte at 18° C., eveluation of equs- 
tion (14) gives 


H,— H,° = 653 m* calories per mole of solute (15) 





17 Drake, Pierce, and Dow, Phys. Rev., 35, p. 613; 1930. Wyman, Phys. Rev., 35, p. 623; 1930. Cutt- 
bertson and Maass, J. Am. Chem. Soc., 52, p. 483; 1930. 

18 Drude, Wied. Ann., 59, p. 48; 1896. Kockel, Ann. Physik., 77, p. 430; 1925. 

i9 Lange and Robinson, J. Am. Chem. Soc., 52, p. 2811; 1930. 
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| This expression gives the variation with concentration of the relative 
§ partial molal heat content of the solute from the region of very low 
» oncentration to infinite dilution. 


The expression for the like property of the solvent, H,O, can be 


F deduced from the partial molal equation.” 


55.508 dH, +m dH, = (16) 


On substituting in equation (16) the value for dH, obtained by dif- 


ferentiating equation (15), and then integrating, one obtains for an 
aqueous solution of a uni-univalent elec trolyte at 18° C, 


=. 7 _(A- H,*) m_ _—653 m% 
A, — A 3 (55.508)  — 3. (55.508) (17) 


= —3.92 m”* colitis per mole of H,O 


For a solution composed of m moles of solute and 1,000 g or 55.508 
moles of H,O, equation (9) may be written 


55.508 


&— 4,2 = H,— Ae + 55508 (A, — 9) (18) 


(H,—H,°) m 


3 (55.508) for (H,— H,°) (from 


Substituting into equation (18), - 


equation (17) ) gives 
9 = aie 
}, — &,° =3 (H,— H,°) (19) 


This relation always holds whenever a partial molal property of the 
solute is directly proportional to m”. 
Evaluating equation (19), one obtains 


&, — &,° = 435 m” calories per mole of solute (20) 


This equation gives the theoretical variation of the relative apparent 
molal heat content with concentration in the very dilute region for a 
ni-univalent electrolyte at 18° C., and is the device us sed in this 
aie to extrapolate the experimental data from regions of low con- 
centration to infinite dilution. 
The uncertainty in the value of the numerical ceefficient in equa- 


ae ; dD 
tion (20) is in large part due to the uncertainty in the value of aT 


1 


used in evaluating equation (14). If the uncertainty * 


dp. 
of qT is taken as about 3 per cent, the resulting uncertainty in the 


in the value 


value of the numerical ccefficient in equation (20) is about 50 calories. 
However, since this equation is used in extrapolating only over a 
small range, the uncertainty due to lack of knowledge of the exact 
value of this coefficient introduces a possible error in the extrapolated 
values of only several calories in most cases. As can be seen from an 
inspec tion of Figures 1,2 Y and : 3, W here are plotted the experime ntally 


“20 See f ootnote 3, p. 792. 21 See footnote 19, p. 794. 
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observed values for the heats of dilution, the error in locating th BT 
intercept at infinite dilution is due more to lack of knowledge of th. He: 
exact position of the curve in the region of low concentration than t) & re 
the selection of the slope used in extrapolation. 
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Figure 1.—Plot of the experimental data on the heats of dilution of HCl, 
LiCl, and OH 


The data have been calculated to 18°C., in calories per mole of solute. 

The ordinate scale measures the apparent molal heat content in calories per mole while the abscissa seal 
gives the square root of the molality. Each unit of the ordinate scale is for HCI, 200 calories, and for Li(! 
and LiOH, 100 calories. 

The data of the various investigators are: 


@ Richards and Rowe. O Wrewsky and Savaritzky. 
0) Richards, Mair, and Hall. @ Nernst and Orthmann. 
@ Muller. X Nernst and Naudé. 

A Steinwehr. A Lange and Durr, 


III. METHODS OF CALCULATION AND TREATMENT 
OF DATA 


The starting point for the evaluation of ,—,°, H.— H,°, and 
H,— H,° lies in equatior. (7) where the quantity 4H can be deter- 
mined by proper extrapolation of the experimental data to infinite 
dilution. This may be illustrated by the following procedure: _ 

_ Given the solution HCl:25H,0. The heat content of this solution 
is 

H=77,+ 25, (21 
or 

H=%,+25H,° (22) 
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© The heat change which occurs when this solution is diluted, for 
© example with an equal amount of H,O, is measured. We have the 
a reaction 


HCl-25H,0 + 25H,0 = HC1-50H,O (23) 


The heat content of the newly formed solution is H’. 
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Fiaure 2.—Plot of the experimental data on the heats of dilution of NaOH, 
KOH, HNO;, and LiNO, 
The data have been calculated to 18°C., in calories per mole of solute. ‘ 
lhe ordinate scale measures the apparent molal heat content in calories per mole while the abscissa scale 
gives the square root of the molality. Each unit of the ordinate scale is 100 calories. 
lhe data of the various investigators are: 
@ Richards and Rowe. @ Richards and Gucker. 
A Pratt. 


Then 
AH = ©’, + 507,° — [(&, + 25H,°) + 2577,°] = &,’ — , (24) 


One then proceeds to measure the heat change when this new solution 
is further diluted: 


HCl1-50H,0 + 50H,O = HCI- 100H,O 


This gives 
AH = ©,/’ — ®,’ 
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Combining with equation (24), one obtains for the reaction 
HCl1-25H,0 + 75H,O = HCI-100H,O 


the heat change 
AH = ®,,’’ —~ P;, 
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Ficure 3.—Plot of the experimental data on the heats of dilution of NaCl, KCi, 
NaNO, and K NO; 


The data have been calculated to 18° C., in calories per mole of solute. 

The ordinate scale measures the apparent molal heat content in calories per mole while the abscissa scai 
gives the square root of the molality. Each unit of the ordinate scale is for NaCl and KCl, 100 calories, and 
for NaN Ox3 and KN 03, 200 caiories. 

The data of the various investigators are: 


@ Richards and Rowe. @ Nernst and Naudé. 
@ Lange et al. X Nernst and Orthmann. 
A Pratt. 


Proceeding in this manner, one can measure, with sufficiently accurate 
and sensitive apparatus, jheats of dilution to very small concen- 
trations. When a proper extrapolation to infinite dilution is made, 
then, for the reaction 


HCl-25H,O+ (o —25) H,O=HCl- oH,O 
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' the heat change becomes 
AH = H°—H=%,°—&, (30) 
' or for the reverse change 
&, — &,° = H— H° (31) 


| where 1 mole of solute is concerned. 

| When the data obtained for &’,—%,, &’’,—®,, &,’’’ —®,, etc., are 
| plotted against the square root of the molality and extrapolated by 
“means of equation (20) to infinite dilution, values of 6,—®,° can 
be read off. 

Where sufficiently accurate values of the apparent molal heat 
| capacity of the solute are available, the calculation to 18°C. of the 
| heat of dilution data obtained at other temperatures introduces no 
uncertainty. The temperature coefficient of the apparent molal 
heat content is equal to the apparent molal heat capacity: 


d& 
aT = & (32) 


7 (},—,°)=%,— 6,° (33) 


Having obtained values for &,—®,°, the values for H,—H,° and 
H,—H,° can be derived by a procedure similar to that employed by 
Randall and Rossini” and Rossini” for heat capacities. 

From equation (28), one may write, for m moles of solute, 


H—H°=m (®,-—®;,°) (34) 


Differentiating with respect to m, 


1 d 
Jy HLH?) = (®n— © x°) + my (H,— #)°) 


But fi (H— H°) = ff, — H,° 


Hence 


d 
H,—1 2° =b,—O,° + ma (,—,°) 


> 


substituting 2m d m™ for dm, we obtain 


, 14 d € > 
H,° =©,— 4° + m4 7 (@a— 2°) (38) 


From a plot of values of 6,—®,° against m™, one obtains values 
for H,—H,° from equation (38). 


# See footnote 10, p. 793. 23 Rossini, B. S. Jour. Research, 4, p. 313; 1930. 
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Now by equation (18) 


%,—4,°=H,— H+ (H,— H,°) (18) 





Combination of equations (38) and (18) gives 


3% 
H,—H,°= “55 508) ge (f2— ®,") (39) 





Values of the slope of 6,—,° against m™ have already been obtained 
in evaluating H,—H,°. There remains only the multiplication of 
m% 


—3(65.508)’ to obtain values for 


this slope by the coefficient, 


H. 1 —H, oft 
IV. THE AVAILABLE DATA 


The data which have been used in this paper to carry out the 
procedure and calculations described in the foregoing section are 
plotted, for the various solutes, in Figures 1, 2, and 3. The plotted 
data have, of course, all been calculated to 18° C. from whatever 
temperature they were measured. Only those data which give heats 
of dilution in the range from MX .-25H,0 to infinite dilution were 
considered. 

Richards and Rowe “ measured the heats of dilution of HCl at 16° 
and 20° C. over the range HCl-25H,O to HCl-400H,O. These same 
authors “ obtained data at the same temperatures and for the same 
range of concentration on HNO;, LiNO;, NaNO;, KNOs, LiCl, 
NaCl, KCl, LiOH, NaOH, and KOH. For the purposes of the 
present work, the measurements at 16° and 20° C. were averaged to 
give the heat of dilution at 18° C. 

Richards, Mair, and Hall * repeated the measurements on HCl and 
extended them to HCl-800 H,O. Richards and Gucker * reported 
new measurements on NaOH, extending the data to NaOH -3200H,0. 

Pratt * measured the heats of dilution of NaCl, KCl, NalO,, 
NaOH and KOH over the range of concentration MX-25H,0 to 
about MX -400H,O at many temperatures from 0° to 33.5° C. For 
NaCl and KCl, data at 18° C. were available, but for the other elec- 
trolytes data for 18° C. were obtained by interpolating the measure- 
ments made by Pratt at 15° and 20° C. 

Data on HCl were available from the measurements of Muller” 
at 16° C., of Steinwehr™ at 18° C., and the data of Wrewskii and 
Savaritzky,*! who measured the heat of solution of HCl gas in water 
at 21.5° C. for various concentrations, 

In the very dilute range of concentrations there are the measure- 
ments of Lange and coworkers at 25° C. and of Nernst and coworkers 





* Richards and Rowe, J. Am. Chem. Soc., 42, p. 1621; 1920. 

2% Richards and Rowe, J. Am. Chem. Soc., 48, p. 770; 1921. 

% Richards, Mair, and Hall, J. Am. Chem. Soc., 61, p. 727; 1929. 
27 Richards and Gucker, J. Am. Chem. Soc., 51, p. 712; 1929. 

% Pratt, J. Frank. Inst., 185, p. 663; 1918. 

29 Muller, Bull. Soc. Chim., 13, p. 1053; 1913. 

8 Steinwehr, Z. physik. Chem., 38, p. 185; 1903. 

st Wrewskii and Savaritzky, Z. physik. Chem., 112, p. 90; 1924, 
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et 18° C. Lange and Robinson *® obtained data on NaNO; and 
NaCl; Lange and Leighton,* on KCl; and Lange and Monheim,* 
on KNOs. The original work of Lange and Messner® contained a 
few data on LiCl, NaCl, KCl, and KNO;. The measurements of 
Nernst and Orthmann a include data on LiCl and KCl, while Nernst 
and Naudé * record some measurements on LiCl, NaCl, KCl, NaNO,, 
} and KNO3. 

All of the foregoing data, where obtained at other temperatures, 
were calculated to 18° C, by means of the apparent molal heat capac- 
ity values already mentioned.* The experimental data for the 
heats of dilution have been plotted against m” in Figures 1, 2, and 3. 
The extrapolation to infinite dilution has been made by the use of 
equation (20). 


V. THE CALCULATED RESULTS 


By a shift of the ordinate scale so that the intercept of the curve 
for the heat of dilution on the ordinate axis becomes the origin, 
m=0 and &,—%,°=0, the curves in Figures 1, 2, and 3 give immedi- 
ately values for &,— ,° , by equation (30). 

In Table 1 are given values of &,—%,°, the relative apparent molal 
heat content of the solute, for the 11 different solutions at dilutions 
expressed in terms of the number of moles of H,O associated with 1 
mole of solute. 

From the curves of ®,—,° against m™%, values of H.—H,2°, the 
relative partial molal heat content of the solute, were determined by 
means of equation (38). Then by application of equation (39) 
values for H,—H,°, the relative partial molal heat content of the 
H,0, were obtained. In Table 2 are given values of &,—,°, H2—H 


and ,—H,°, for the 11 different solutions at concentrations from 
m=0 to m=1.96. 


TABLE 1.—b,—,°, relative apparent molal heat content, at 18° C. in g-cal.is per 
mole of solute + 


Solution 


| | 
MX nc | LiCl |NaCl) KCI NO; Nal =N a Si 








MX: 0 0 0 
MX. 6,400 H20 38 35| 3 : : é 38 
MX: 3,2 ; ; 45 ‘ g 2g 5 50 
MX- 1,600 H20 : 55 f f : { 77 63 

* 800 H20 | 115 61 i 7 5 3 7 76 


61 3 , ¢ 86 113 
47 5 —86 x 78 130 
—5 | q — 246 —465 39 137 
—120 |—128 —520 | —885 —45 137 
450 |—330 |—318 28 5 | —960 |—1, 530 3 | —173 162 
































e duu composition of the solution is given in terms of the number of moles of water associated with 1 mole 
f solute, 


* Lange and Robinson, Z. 2 pot Chem., ee 148, P. 197; 1930. 
‘ Lange and Leighton, Z lektrochem., 34, p. 566; 1928. 
‘ Lange and Monheim, Z. physik. Cnn ,~ 150, Pp. 349; 1930. 
’ Lange and Messner, Naturwissenschaften, 15, p. 521; 1927. 
* Nernst and Orthmann, sitzb. preuss. Akad. Wi iss., D. 51; 1926; p. 136; 1927. 
Nernst and Naudé, Z. Elektrochem. , 33, p. 532; 1927 “ph 

* See footnote 9, p. 793. 
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TABLE 2.—Apparent and partial molal heat content of solute and partial molq| 
heat content of H,O in aqueous solutions of HCl, LiCl, NaCl, KCl, HNO, LiNo,, 
NaNO;, KNO;, LiOQH, NaOH, and KOH at 18° C : 





m | | ™m 
moles Dp-Pp° | Fe-FI9° f,-Hi,° moles Pa-Pa° | A-F7,° | T,-Fi 
of so- | g-cal.is| g—cal.is g—cal.is of so- g-cal.is} g-cal.is | g-cal, 
lute per per r Solute lute per pe per 


pe: 
per mole of| mole of | mole of per mole of mole of 
solute | solute H30 yw solute H;0 
g He 


1,000 
g H:0 
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The data for #,—®,° and H,— H,° are plotted against m* in Figures 
sand 5, respectively, while Figure 6 shows values of 4: — H,° plotted 
. 36 
pagainst m”. 








a | 
| 

















CALORIES 


= 
tome 

















NN 
\ 
< 
































2 4 6 ° 
m2 
Figure 4.—Plot of ®)—®;,°, the relative apparent molal heat content of the 
solute, in aqueous solutions at 18° C. 


lhe ordinate scale gives ®,-®4° in calories per mole of solute while the abscissa scale gives the square root 
{the molality. 


VI. DISCUSSION 


I'he character of the family of curves shown in Figure 4, in which 
the relative apparent molal heat content is plotted against the square 
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root of the molality, has already been discussed, in the very dilute 
region by Lange® and coworkers, and in the more concentrated range 
by Richards and Rowe.” The curves in Figures 5 and 6, showing the 
relative partial molal heat content of the solute and of the H,0, ar 
self-explanatory. 
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FravureE 5.—Plot of H:-Hi2°, the relative partial molal heat content of the 
solute, in aqueous solutions at 18° C. 

The ordinate scale gives Hr-Fi’ in calories per mole of solute while the abscisss scale gives the square 
root of the molality. 

. . 0 

The usefulness of a complete compilation of values of ®,—®,, 

H,— H2°, H, — H,°, for the given concentration range will be apparent 

to whoever has need of such data in the calculation of heats of reac- 

tion at various concentrations from a heat of reaction measured 4 





8 See footnote 16, p. 793. # See footnote 6, p. 793. 
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Ficure 6.—Plot of Hi-Hi°, the relative partial molal heat content of the 


H20, in aqueous solutions at 18° C. 


The ordinate scale gives #f:-Fi:° in calories per mole of H30 while the abscissa scale gives the 3/2 power of 


the molality. 











806 Bureau of Standards Journal of Research [Vol 





one concentration, or in the calculation of the temperature coefficien; 
of the partial molal free energy of the solute or of the H,O. The heat 
content values given here for 18° C. can be converted to other tem. 
peratures by the appropriate use of the apparent and partial molal 
heat capacities of the substances involved. 


WASHINGTON, January 23, 1931. 
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ACCURATE MEASUREMENT OF SMALL ELECTRIC 
CHARGES BY A NULL METHOD ! 


By Lauriston S. Taylor 


ABSTRACT 


In the use of an electrostatic system for measuring charges and currents, it is 
necessary to know the electrostatic capacity of the system. For small capacities 
the error in this measurement may easily be 1 per cent. There is here described a 
new method for calibrating a null system in such a manner that the capacity of 
the leads does not enter and which, therefore, permits a reduction in the calibra- 
tion error to one-tenth. When a system is once calibrated in the manner de- 
scribed, any unknown capacity whatever may be added to the leads without 
affecting the measurement of the desired quantities. Expressions are given for 
the sensitivity of the system in terms of readily measured quantities. Applica- 
tions to the measurements of current, charge and capacity are discussed. 


CONTENTS 


I. Introduction 

Il. Theory of method 
1. Measurement of charge 
2. Calibration of condenser - 
3. Sensitivity of system 
4. Method of operation-_-- -_- 

III. Experimental study 
1. Description of system as ; used. 
2. Calibration 
3. Factors affecting se ensitiv ity. 

Application to current measurements... --- 


I. INTRODUCTION 


When making accurate measurements of small charges or currents 
by electrostatic means one of the greatest difficulties encountered is 
that of correctly determining the capacity of the system. This is due 
largely to the inability, without taking elaborate precautions, to 
determine corrections for the capacity of the leads from the capacity 
bridge to the electrometer system. ‘These errors have been discussed 
in detail by Rosa and Dorsey ? who showed that for capacities of the 
i x of 500 cm they might be as much as 1.0 per cent. This diffi- 

ty has been enc ‘ountered i in endeavoring to make accurate measure- 
ments of the ionization produced in air by X rays, in connection with 
the international unit of X- -ray quantity—the roentgen. 

Investigators have used several different means of measuring the 
small currents involved. The factors favoring the use of an electro- 

tatic measuring system have been discussed and it has been shown 3 * 

t the null electrostatic method, 1 in that it does not introduce field 





) footnote , Rac ei 18, p. 9; 1931. 
B 2osa and N. E. Dor: ey, B.S. Bull., 3, p. 433. 
silla, Am. J. Roent., 21, p. 47; 192% 
Taylor, B. S. Jour, Research, 2 (R Pp 56), Pp. 771; 1929. 
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distortion in the ionization chamber, has an advantage over deflection 
methods. An additional difficulty in deflection electrostatic systems 
of small capacity arises from the motion of the leaf or fiber altering 
the capacity and thereby the calibration. 

In null electrostatic measuring systems the potential of the electro- 
meter and collecting system is maintained at some constant value— 
usually zero with respect to ground—by one of three methods. One 
method of accomplishing this is by communicating to the collector 
system an evaluated charge or current opposite in sign to that meas- 
ured. Failla® and Jaeger * have devised sources of current for this 
purpose which serve as working standards. For absolute measure- 
ments this method is impractical because the working standard 
requires a calibration in kind. In the second method, the charge 
accumulated on the collector system is compensated by increasing its 
capacity to maintain the potential constant. Aside from the mechan- 
ical difficulties involved in a circuit of this kind, it has the disadvantage 
necessitating an elaborate calibration. In a third method, sometimes 
known as the Townsend method,’ the communicated charge is entirely 
localized in a condenser of constant known capacity, but varied 
potential difference as determined by a voltmeter. In this, the 
greatest source of uncertainty is the voltmeter while in the second 
method it is the variable condenser, Since, in general, a variable 
condenser suitable for such purposes is a less dependable instrument 
than a voltmeter, the third method is from this standpoint distinctly 
preferable to the second. 


II. THEORY OF METHOD 
1. MEASUREMENT OF CHARGE 


Given a circuit as represented in Figure 1, an unknown charge is 
communicated to the system containing an accurately known capacity 
C; and the unknown remainder C (which represents the capacity of 
the leads, ionization chamber, and electrometer, together with a cali- 
brated variable capacity which will be described below). One plate 
of Cis connected to earth, while the corresponding plate of C, is con- 
nected to a potentiometer circuit so as to bring it to any desired 
potential above or below that of the earth. As indicated by the 
deflection of the electroscope # this communicated charge raises the 
insulated side to a calibrated potential v relative to the earth as zero. 
Calibration of the electroscope is readily carried out by closing the 
key K and shifting the potentiometer contact to give any desired 
potential v which is read directly on the voltmeter. (The electroscope 
can at the same time be adjusted so that the deflection is at the most 
sensitive part of the scale.) With the electroscope at v, K is opened 
and an unknown charge Q communicated to the system through, say, 
the ionization chamber J. This requires the potentiometer contact 
to be shifted to V in order to bring back to, or maintain the potential 
of the system, at the calibrated value v, as indicated by the electro- 
scope. 

Since v is unchanged by this operation the charge on the part of the 
system of capacity C remains unchanged, hence the entire charge 





5 G. Failla, Radiology, 15, p. 437; 1930. 
6 R. Jaeger, Strahlentherapie, 33, p. 542; 1929. 
7J.8. Townsend, Phil. Mag., 6, p. 598; 1903. 
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is accumulated on C;. Its magnitude is readily obtained from the 
measured values; that is 
Q=C, (v—V) (1) 


2. CALIBRATION OF CONDENSER 


The difficult part of this method lies in the accurate calibration of 
(,. To present an adequate method of calibration is the prime object 
of the present communication. 

Cis here made up, in part, by an accurately calibrated variable con- 
denser connected into the system as dotted in the figure. The total 
capacity in situ of this condenser is not accurately known because the 
leads, etc., introduce an uncertainty. However, its capacity differ- 
ences are very accurately known. 









































FigurE 1.—WNull electrostatic system for measuring small charges or currents 


The part within the heavy lines is referred to as the ‘‘isolated system” and is so constructed that all 
capacities remain fixed. 


Starting now with the system uncharged (that is, when V is zero, 
K is first closed then opened again), the potentiometer contact is 
shifted to a value V such that the electroscope indicates accurately 
any convenient potential v. Since no charge has been communicated 
to the system by this operation the induced charge (—Q) on the 
unknown capacity C, must be exactly equal and opposite to that (Q) 
on the unknown remainder of capacity C. 


Since 
—-Q=C, (— V) 
Q=Cv 
C,(V —v) = Co 


and 


it follows that 
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The experiment is now repeated with C changed by an accurately 
known increment to a new value C’, and V, consequently to a new 
value } 

In this case 


CQ, (V’—v) =C'v (3 


Subtracting equation (3) from equation (2), C;, the desired quantity, is 
obtained in terms of the accurately known magnitude (C—(’), 
namely 


C,= (C—C¢ ) y- |’ (4) 


from which expression all extraneous capacity, such as the ionization 
chamber, electroscope, and leads are wholly eliminated. Since the 
capacity of the system varies with the deflection of the electrometer, 
the same deflection must be used throughout any one calibration. 
After having determined C, other capacities of unknown value may 
be inserted in the collector system without in any way affecting its 
calibration. This is of great practical importance, as shown later 


3. SENSITIVITY OF SYSTEM 


Although, as seen in equation (1), the quantity of the charge meas- 
ured does not depend on the magnitude of the distributed capacity (, 
yet the sensitivity of the measuring system must obviously decrease as 
C increases. 

8 ° : ° ° . 
0 of this type of null circuit can be expressed in 
terms of the sensitivity of the integral parts; that is 


ee ee 
The sensitivity A 


As_ As Av 
AQ Av AQ 


where As is the incrementinscale divisions of the electrometer deflection 
for a given increment AQ of the imparted charge; and at the same time 
Av is the increment in volts corresponding to the deflection increment 
As. 

As 


The electrometer sensitivity, — »is a magnitude which is obtainable 


separately, is in no way dependent upon the remainder of the system, 


and may be considered as known. As to the other factor he 


Av 

bd rt 
AQ 
sensitivity of the compensating system, we have (from the general 
law Q=CV) AQ=CAV. Any error Av in the adjustment of the 
potential v introduces an error AQ, which is given in this case by 


AQ=(Q,+ C)Av (6) 


That is, the error in the determined magnitude Q is proportional to 
the total capacity of the system. 

The uncertainty, Av, may readily be obtained by comparing the 
electroscope independently with an adequate voltmeter; the value 
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of C, is supposed to be accurately known; leaving C, which includes 
all capacity of the circuit external to C, to be evaluated in terms of 
observed magnitudes. 

From equation (2) 


O+0_V 


Qe 


hence equation (6) may be rewritten 


y 


AQ = C; i Av (8) 


a 


which in turn, since from equation (7) : is a constant (say k) of the 
system as given, may be written 
AQ=CkAv 
The sensitivity of the compensating system is then 
fe § 
AQ Gk 


in which l/k is a sensitivity factor of the system as given. Since C; 
is constant and known from (4) we have now expressed the sen- 


(9) 


Av 
AQ 
sitivity of equation (5) then becomes 


As_ As 1 
AQ Av Gk 
which, since & increases with the stray capacity C, shows that the 
sensitivity decreases as the stray capacity increases. 
By considering all capacities outside the heavy lines of Figure 1 to 
be removed from the system, we may determine a value of the sen- 
sitivity factor k, which is characteristic of the isolated part. Thus 


J y 
ar 0 


Vo 


sitivity without C being explicitly involved. The working sen- 


(10) 


ky (11) 


serves in practice as a ready control of the system. 
4. METHOD OF OPERATION 


Referring to equation (4), only two capacities whose difference is 
accurately known are required for the calibration of C;. If, however, 

variable capacity C having a number of accurately known capacity 
differences be used for the calibration, the following graphical anal- 
ysis which gives equal weight to all observations shows readily how 
— the magnitude of any of the capacities in the system, in- 
cluding C. 


Remembering that i i k, we have from equation (7) 


C=Ok—-C, 
49527°—31——_-4 
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showing a linear relationship between C and k. Here it is seen that 
, may be obtained either from the slope or from the C intercept; 
also that k=1 when C=0; and finally that for k=k,, C takes on q 
derived value which will later be found serviceable in the analysis, 
In plotting C against k, as in Figure 2, the capacity differences for 
various settings of C being accurately known while the absolute 
capacity is yet undetermined, we are unable to assign C its correct 
coordinate position at the start but merely a place, say a. Then, 
for a C scale, may be chosen coordinates representing capacity dif- 
ferences employed without regard to the location C=0. 
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Figure 2.—Schematic method of evaluating fixed capacity C, 


The next point b corresponding to a known capacity change from 
a is then plotted with reference to a, etc. A straight line drawn 
through the plotted points must also pass through the point k=1, 
’=0 and thus locate the origin. Having determined the position 
C=0 it is now possible to read off the graph the value of the other 
capacities in the circuit. For example, the isolated system has 4 
determined sensitivity factor k, so that referring to (12) the corre- 
sponding ordinate C, gives what will be called the stray capacity 0! 
the isolated system. Similarly, the difference between the ordinate 
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Figure 3.—Photograph of “isolated system” as used in the study 
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C, and that C, at any other point P on the curve gives the corre- 
sponding capacity of the external system. 

Obviously such a system is applicable to the accurate measure- 
ment of any inserted capacity; for example, (C,—C,). We will 
consider C as made up of the stray capacity C, of the isolated system 
and the total external capacity C.,,. Then equation (12) becomes 


nk =CQ,+C,+ Cost (13) 


In this C, is determined by equation (4), C.x, is removed, and Q, is 
determined from the relation 


Cho =(C,+ C, (14) 
Combining these 


Cort = C, (k— ke) (15) 


For any given measuring system, C,; and C, and k, remain fixed 
so that the measurement of any capacity, such as C,,:, requires but a 
determination of k with the external capacity inserted. 


III. EXPERIMENTAL STUDY 
1. DESCRIPTION OF SYSTEM AS USED 


We will now examine a specific null circuit represented in general 
by the diagram in Figure 1 and for the isolated part in Figure 3. 
In this latter, His a string electrometer of the general Edelmann type 
in which the deflection of a very fine platinum (Wollaston) wire, 
under an adjustable tension between the charged knife-edges, is 
observed by a suitable microscope. The potential on these knife- 
edges is supplied, through very high protective resistances R, by two 
22.5 volt batteries. For any given potential on the knife-edges, the 

As 
Av 
tension on the fiber by the micrometer head, 7. The sensitivity of 
this particuler instrument can be readily changed from 100 divisions 
per volt to 0.01 division per volt. 

The resistances of the compensating potentiometer are operated by 
the knobs “coarse” and “fine” on the front of the aluminum box 
container. Potentiometer, charging batteries, compensating con- 
denser C, and the necessary leads are all contained in this box, the 
electrometer circuit being led out through an amber bushing LZ. All 
parts are carefully shielded electrostatically. G is a grounding key 
and AK a switch to short circuit C, (for the purpose of calibrating the 
electrometer). Leads to the potentiometer battery and the volt- 
meter are in the form of jack and plug at M and N, respectively. 


voltage sensitivity of the electrometer is adjusted by varying the 
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TABLE 1] 


VARIABLE CONDENSER ALONE 

















| 
— . Average | neviation 
apacity k | deviation from mean 
C’ 4 from me an | ‘ 
| 
1 2 a 3 | 4 
| | Per cent 
10.0 1.057 | 0.0015 | 0.14 
43.8 | 1. 120 .0016 | 14 
97.6 | 1. 185 0018 | ~15 
175. 9 1. 279 . 0008 . 06 
54. 1, 378 . 0021 15 
| 13.5 1.470 . 0014 . 09 
112.0 j 1. 569 - 0018 12 
490. 6 | 1. 666 . 0009 . OF 
| | — 
| | Average oi} 
' 
i Vv 


ARIABLE CONDENSER + RUBBER 
CABLE 





1. 968 | 0. OOS 0. 36 | 
| 2. O75 | . 005 . 24 
2,232 | . 006 32 
| 
| a ees | 


| Average... -3l 
| 


VARIABLE CONDENSER, HIGH 
VOLTAGE 





| 117.1 | 1.207 | 0. 005 0. 41 

195.6 | 1. 303 . 003 2 

1 ital 1. 397 . 002 14 

| 352.8 1. 495 . 005 ae 

| 431.7 | 41.590 - 004 . 25 
| _——— 

Average... ont 


1 Variable condenser actually removed and L covered. 
? Capacity of variable condenser only. Full value of C is this plus the unknown capacity of the cable 


2. CALIBRATION 


The capacity used for the calibration of C, was a precision variabl 
condenser of well-known make; the voltmeter was of the laboraton 
standard type; the capacity differences and voltmeter were cali 
brated with an accuracy of one-tenth per cent by the electrical divi- 
sion of this bureau. 

The accuracy of the results obtainable under working conditions 
is best brought out graphically. The factor k (equation ‘(9)) may be 
obtained from the slope of the curve V plotted against »—a series 0! 
which curves are given in Figure 4 for several different settings of the 
variable condenser (C. 

The first part of Table 1 gives results of a number of determinations 
of k for different nominal v ralues of the external capacity C, in w hich h 
the variable condenser was connected to the terminal L by a wel 
nal ree and shielded conductor of unknown capacity. The ~~ 
age deviation of five determinations of each & is seen to be within the 
one-tenth per cent accuracy of the measuring instruments. 

Plotting (circles, fig. 5, Curve I) these values of k as abscissae and 
the corresponding variable condenser values as ordinates, it is found 
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Taylor] 


that all of the points lie close on a straight line. According to (12) 
AC ita ; in 
the slope AR’ which is obtainable in terms of the accurate differ- 


ences in capacity, gives C, as 809.2+04 uyf. 

Having removed the condenser C and covered the amber bushing 
L with a cap, the sensitivity factor /, for the electrometer-condenser 
system alone was found to be 1.057 + 0.0005. 
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FiaurE 4.—Curves showing the evaluation of the sensitivity factor k over a range 
of potentials 


This curve might be given its coordinate position on a C-k plot by 
extending it to the point k=1; hence, also, the ordinate value C=0 
located as seen at —48 on the variable condenser scale. The value of 
the internal lead capacity C, may be read from the curve at the point 














$16 


k = 1.057, thus giving C,=46 wyf. Likewise the value of C,, at the 
point P for example, is found to be 460 u wf. 
These capacity values may be calculated more accurately from 
equations (14) and (15). Using equation 14 for k,=1.057 we find 
1+ C,=855.3 uw uf where, since C,= 809.2 uw uf, we have C,=46,1 
uu. The capacity at C,, for which k=1.569, is obtained by substi- 
tuting the values of k and k, in equation (15) giving C,=459.8 u yj. 
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3. FACTORS AFFECTING SENSITIVITY 


The second part of Table 1 gives values of & for three nominal 
values of the variable condenser, to which was connected a shielded, 
rubber, single-wire cable of unknown capacity, with the shield 
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Figure 5.—Curves showing the evaluation of the fixed capacity C, by means 
of an accurate variable external capacity 


Curve I (circles) using low voltage, (crosses) using high voltage. Curve II using a rubber cable in the 
external system. 


grounded to serve as a condenser having a poor dielectric. The 
fourth column of this part compared with the part above shows that 
the precision of measurement was not so good—the average deviation 
changing to 0.31 from 0.11 per cent which, as already stated, was 
the error of calibration. The greatest error (0.31 per cent), however, 
is due to the experimental difficulty in determining k; for, due to 
dielectric absorption in the rubber cable, it was difficult to fix with 
accuracy upon the V necessary to apply to C, to induce a definite 
potential v 0 on the electrometer.’ 

Reference to the corresponding Curve II in Figure 5 shows the same 
effect in that the plotted points determined with the cable in the sys- 





§ This absorption is not operative when using, as in current measuring, the system at »=0. 
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tem do not fall closely on the straight line drawn parallel to Curve I. 
The superiority of conditions when all questionable external capacity 
is removed while calibrating C,° is obvious—a conclusion supported 
also by the fact that, when the cable was replaced by a highly insu- 
lated condenser of about the same capacity, the points fell well along 
a straight line, the slope of which represented satisfactorily the ca- 
pacity as checked by other means. Furthermore, the average devia- 
tion from the mean for the values of k made with the well insulated 
unknown was 0.12 per cent as against 0.31 per cent for the cable. 
This small average deviation, being again of the order of the error 
of the instrument calibration, indicates that the larger deviation 
(0.31 per cent) for the rubber cable was due to absorption and not to 
the magnitude of C. 

In the measurements for the first and second parts of Table 1, the 
electrometer sensitivities were, respectively, about 13 and 17 divisions 
per volt; » was never greater than 1.5 volts, and, consequently, the 
external condenser offered no difficulties due to leakage. However, 
with v increased to 80 volts and the sensitivity decreased to about 
0.3 division per volt, leakage became appreciable and k was subject 
to error as indicated in the third part of Table 1 where the average 
deviation was 0.27 per cent as against 0.11 per cent in the first part. 
The values obtained are plotted as crosses on Curve I (fig. 5). Such 
conditions should, of course, be avoided when determining C,. 


IV. APPLICATION TO CURRENT MEASUREMENTS 


In the measurement of current by this method, a time measure- 
ment, of course, enters. In the particular case of measuring X-ray 
intensities, a shutter in the X-ray beam between tube and ionization 
chamber is so controlled that it remains open for a definite interval 
of time during which the charge is wholly compensated, holding v=0, 
by varying the potential on C.?° 

In practice it is extremely advantageous to use a shielded rubber 
cable for connecting the electrometer to the source of current. As 
shown above, the presence of unknown lead capacities does not affect 
the measured value of the charge, although it does affect the accuracy 
of the measurement. 

The effect of poor insulation may be eliminated by choosing such 
conditions that »=0 with respect to the shielding. This confines 
leakage to the condenser C;, and it is a relatively simple matter to 
construct this fixed capacity so that leakage between its plates will be 
negligible. In the isolated part of the system here described, a 
20-volt charge on the condenser leaked to ground at the rate of only 
3 per cent in 8 hours; but, with a rubber cable attached, the same 
leakage occurred in about 10 minutes. 

To obtain the desired sensitivity with any length of cable added 
to the system, the electrometer sensitivity may be changed in accord 
with & as brought out in equation (9). 

To test the effect of a poor dielectric in the external condenser, 
two sets of intensity measurements were made of the same X-ray 





‘It is of interest to note that using equation (15), the capacity of the cable was found to be 522 uuf as 
compared with 526 yuf obtained with a capacity bridge calibrated with an accuracy of 1 per cent. 

“5 D. L. Webster and A. E, Henning, Phys. Rev., 21, p. 301; 1923. 

1 L. 8. Taylor, B, 8, Jour. Research, 2 (RP56), p. 771; 1929. 
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beam, using a standard X-ray ionization chamber.” In the firs; 
set, the electrometer and ionization chamber were connected by an 
amber insulated conductor having a capacity of about 206 uf and 
negligible leakage. In the second set, the connection was by the 
shielded rubber cable having a capacity of about 525 puf. For this 
the electrometer sensitivity was adjusted to accord with the change 
in k, approximately double that of the first case. The respective 
average deviations from the mean of 10 observations for each of the 
two cases were found to be 0.33 and 0.32 per cent, respectively, the 
averages of the two runs differing by 0.13 per cent. These deviations 
from the mean were no more than should be expected from the X-ray 
intensity variations. 

In a null circuit previously described by the author ® the expres- 
sion for the magnitude of the ionization current involved the total 
capacity Cr of the system instead of simply C, as above. The fact 
that the total capacity Cy; can not be measured accurately (with a 
bridge, as previously done) renders the other method less satisfactory 
than the present one. Intensity measurements of the same X-ray 
beam by the two outfits showed a difference of nearly 1 per cent which 
is principally attributed to the error in the measurement of C;. 
Compared with the average deviation above of one-third this magni- 
tude for both the highly insulated and the cable connected systems, 
the superiority of the method here presented is evident. 


Wasuineton, December 30, 1930. 
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LIGHT SENSITIVITY OF ROSIN PAPER-SIZING 
MATERIALS? 


By Arthur E. Kimberly” and J. F. G. Hicks? 


ABSTRACT 


The effect of light upon rosin, and compounds of rosin formed during the rosin 
sizing of paper has been studied, both in air and in an atmosphere of nitrogen, a 
relatively inert gas. The source of light used was an inclosed carbon are. Com- 
mercial rosin and bleached rosin as well as resinates of sodium, aluminum, and 
iron, and rosin size contaminated with ferrous and ferric iron were exposed to the 
rays of the are in quartz tubes. The fact that changes occur in an inert atmos- 
pl here, such as nitrogen, indicates that the presence of oxygen is not necessary to 
the chemical alteration of rosin when exposed to the action of light. 

The more important conclusions of this investigation follow: Rosin, both 
bleached and unbleached, is sensitive to light, darkening under the action of the rays. 
Ferrous resinate, and rosin size contaminated with ferrous iron are also light 
sensitive probably because of the ease with which the ferrous iron is oxidized by 
light. Ferric resinate, aluminum resinate, and sodium resinate have no effect 
upon rosin size in so far as light sensitivity is concerned. The yellowing of paper 
with age may be explained in terms of the light sensitivity of rosin and of ferrous 
ron. 


CONTENTS 


I. Introduction 
II. Experiments performed 
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2. The effect of light upon resinates of sodium, aluminum, ferrous 
iron, and ferric iron 
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3. The effect of light upon rosin size contaminated with ferrous 
and ferric iron 
IV, Summary 
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>. 


I. INTRODUCTION 


Light, according to the literature, is an important factor in the 
deterioration of high-grade papers, as is most frequently noticed in 
libraries or similar depositories of records. The term ‘‘high grade’ 
is used in this instance to distinguish papers containing no ground 
wood fibers from those containing ground wood fibers. The change i in 
color occasioned by the action of light upon papers containing ground 
wood has been attributed * to the “oxidation of the lignins contained 
in such fibers. In order t to ascertain the effect of light upon the ma- 





' This series of ‘Geietieiians is being ms vie at the request of the American Library Association, and is 
anced by a grant from the Carnegie Institution. The experimental work is being performed at the 
National Bureau of Standards. 


fir 
\ 


* Research associate, National Research Council. 


Ka on, P.: Svensk P appierstidning. See translation by C. E. Bandelin, Paper Trade J., 57, pp. 17, 
1913, Compare Cross, C. F., Roy. Soc. Arts J., 45, p. 684; 1897, 
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terials comprising the rosin sizing used in the manufacture of papers 
containing no ground wood, as distinguished from the effect of light 
upon the paper fibers themselves,‘ this phase of an investigation of 
the impermanency of records was inaugurated. 

Klemm ° indicated that the yellow discoloration in paper was due 
to what he called “rosin-iron soaps,” and that the intensity of the 
discoloration was generally proportional to the quantity of iron 
present. Schwalbe*® found that rosin itself was very sensitive to 
both light and air, and that a rosin-sized paper, upon exposure to 
direct sunlight, underwent a rapid loss of sizing. Schoeller’ indi- 
cated that all papers turned yellow to a greater or lesser extent, with 
passage of time, more especially those which had been sized with 
rosin. He further indicated that the darkening of rosin itself was 
due to polymerization or auto-oxidation, and that the discoloring of 
papers would be prevented by “bleaching” the rosins used in size 
making with an oxidizing agent.2 Zschokke,’ working independently 
of Schoeller, showed results which partially agreed with those of 
Schoeller, but asserted that the color changes in paper are due to 
the effects of heat and light. He further disagrees with Klemm” 
as to the direct relation between degree of discoloration and iron 
content of paper. At about the same time, Klason "= showed that 
rosin rapidly takes up atmospheric oxygen in the presence of sun- 
light, and only very slowly, if at all, in the dark. Haas ” recom- 
mended the use of ‘‘alum free from iron,” and of ‘bleached rosin” 
to prevent the deterioration of paper. Hitchins * showed that all 
rosin-sized papers will yellow with passage of time, that mere traces 
of iron will influence the depth of color produced, and that the iron 
content of a rosin-sized paper is therefore a very important factor. 
In addition to this somewhat contradictory collection of data, one 
finds in the journals relating to the paper industry several articles 
which appear to be résumés, translations or rewritings of the work 
just cited.’* % 1617 Jt appears that but little attempt has been made 
to connect the various color changes of rosins just described with 
the actual discoloring of paper," or to correlate the discoloration of 
sizes and of paper.’ ” One does not find any attempt to distinguish 
between the behavior of ferrous compounds, ferric compounds, or 
mixtures of the two, in connection with the term “‘iron”’ as related 
to the discolorations discussed. 

Because of the indefinite state of the knowledge concerning the 
cause of the discoloration of rosins and rosin sizes, and the relation 


‘ Trans. Inst. Rubber Ind., 5, p. 31; 1929; also Rubber Chem. and Tech., 3, pp. 1, 57; 1930. See also 
References 1, 3, and 4. 
5’ Klemm Papier, Ztg. Jahrg. 27, 1902, pp, 961-963. 
6 Schwalbe Wochbl, Papierfabr., 38, pp. 1303-1306; 1907. 
7 Schoeller Wochbl, Papierfabr., 48, p. 3222 et seq.; 1912 (serial). 
§ Schoeller’s German Patent, Bleaching of Rosin, No. 257158; 1912. 
® Zschokke Wochbl, Papierfabr., 44, p. 2976 et seq.; 1913 (serial). 
10 See footnote 5. : 
11 See translation by O. E. Bandelin, Paper Trade J., 5%, pp. 17, 46; 1913. Compare also Briggs, Chem- 
ical World 2, p. 119; 1913. 
12 Haas, B., Papierfabrikant, 12, p. 891 et seq.; 1914 (serial). 
13 Hitchins, A. L., ‘‘ Paper’’ pp. 11-15, July 24, 1918. 
14 Pulp and Paper Magazine of Canada, 11, p. 686; 1913. 
15 Paper Trade J., 58, pp. 18, 38; April 30, 1914, 
16 Papermakers’ Monthly J., 56, p. 40; 1918; also ‘‘ Paper,”’ 22, pp. 1, 11; 1918. 
17 “* Paper’’ 30, pp. 7, 70-72; 1922; also Paper Trade J., 74, pp. 15, 313-315; 1922. 
18 “‘ Paper,’’ 30, pp, 7, 70-72; 1922; also, Paper Trade J., 74, pp. 15, 313-315; 1922. 
19 Paper Trade J., 76, pp. 11, 50-52; 1922. 
» Geseboff, L. P., ‘La Papeterie’; March 23, 1923; also ‘‘ Paper,” 31, pp, 26, 143; 1923, 
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between these discolorations and the discoloration of papers, it was 
thought advisable to study the behavior of rosin, bleached * and 
unbleached, when exposed to the action of light, the samples of 
rosins to be placed in separate atmospheres, both chemically inert 
and chemically active with respect to the rosins placed in them. 
Since the rosins, sizes, chemicals, and papers were and always have 
been made and stored in air, these artificial atmospheres must rep- 
resent conditions based upon that fact. The ordinary components of 
the atmosphere were, therefore, used singly and together. 

For the same reason, it was thought advisable to study the behavior 
of rosin and rosin sizes, when intimately mixed with small quantities 
of different types of “iron compounds”’* and placed under the influ- 
ence of light. The usual commercial process of sizing paper with 
rosin is, briefly, as follows: The rosin is heated with a solution of 
“soda ash’? (commercial sodium carbonate) to such a degree as to 
hold the remaining free rosin in colloidal suspension. The size and 
the pulp are well mixed in the beating machine, and a calculated 
amount of so-called “paper makers’ alum” (aluminum sulphate) 
is added. This precipitates the rosin upon the paper fibers, which 
are then run through the machine and formed into sheets. There 
are many modifications of this general procedure, but, for the most 
part, they are minor deviations. By reason of the method of manu- 
facture, rosin sizes, as used in a paper mill, contain the resinates of 
sodium, aluminum, ferrous and ferric iron, in addition to colloidal 
rosin. The precipitating agent for these rosin sizes, aluminum sul- 
phate, nearly always contains traces of both ferrous and ferric iron, 
detectable by the ordinary reagents; hence, the presence, in such a 
size, of ferrous and ferric resinates, is readily accounted for. These 
facts made it necessary to add known small quantities of “iron com- 
pounds”’ to “‘iron free” sizes* and study the behavior of these under 
the influence of light. The addition of traces of solutions of iron 
salts to the liquid ‘‘size milk’? makes possible a very intimate mix- 
ture, similar to that existing on the surface of a commercial rosin 
sized paper. 

The effect of light upon rosin, alone and in combination with other 
materials, was studied. A study of bleached rosin was included in 
the program, since it has been claimed* that bleached rosin does not 
darken upon exposure to light. The source of light in these experi- 
ments was a carbon-are lamp, operated with a direct current of 13 
amperes at 220 volts, at a distance of 10 cm from the samples, and 
inclosed in a globe of special glass. These special glass globes are 
described by their manufacturers as deteriorating at the end of 
100 light-hours’ service; for this reason they were replaced at the 
end of each total 75 light-hours’ service, as an additional precaution. 





‘| Although the use of *‘bleached” rosin has been strongly advocated by some of the German investi- 
gators, it is not known that their recommendations were ever put into practice. So far as is known, no 
‘bleached”’ rosin has been used in the United States. 

* To avoid needless repetition “iron compounds” as used in this discussion will be taken to mean “ fer- 
rous and ferric compounds singly and together.”? Details of making such “iron compounds” will be given 
under the head of “ Experiments Performed.” 

* The term “iron free’? here and elsewhere in this publication means “tested for both Fett and Fet++ 
With negative results.” Further details under the head of “Experiments Performed.” 

* Schoeller, V,, see footnotes 7 and 8, p, 820; Haas, B., Papierfabrikant, 12, pp. 891 and 919; 1914, 
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The spectral distribution of radiation from such a glass-inclosed are 
as this, as compared with sunlight, is as follows:* 





Percentage of total 
radiation 
Spectral range in uu 











{ 
Arc | Sun 
ae rn 7 eee one SSS he 
Re ee Joa ee 0.0 2.0 
| SSR re es a 2.0 2.8 
360 to 480___- : Kaba 18. 5 12.6 
480 to 600____ man «ane 9.3 21.9 
1,400 to 4,200- . - : a 22.1 21.4 
4,200 to 12,000_- | 31.6 4 
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II. EXPERIMENTS PERFORMED 
1. EFFECT OF LIGHT UPON ROSIN 


Four varieties of rosin were used in this phase of the work: 

(a) Unbleached rosin, grade F, such as is commonly used for sizing 
paper. 

(b) Unbleached rosin, grade C, lighter in color than grade F, 
Used in papers when a lighter color is particularly important. 

(c) Rosin, grade F, bleached as follows: The rosin was melted on 
the steam bath, and chlorine gas passed through the melt until the 
color of the rosin no longer lightened on continued passage of chlorine. 
The resulting mass was then converted into sodium resinate by the 
addition of an excess of sodium-hydroxide solution, and the mixture 
heated on the steam bath until a clear liquid resulted. This sodium- 
resinate solution was then acidified by adding an excess of strong 
hydrochloric acid, which precipitated free rosin therefrom. This 
free rosin was filtered off, washed twice with distilled water, and dried 
in contact with air. This is Schoeller’s method. 

(d) Rosin, grade F, bleached as follows: The rosin was mixed with 
a slight excess of sodium carbonate solution, and the mixture heated 
on the steam bath until a clear liquid was obtained. Chlorine gas 
was then passed through this liquid until no further lightening of 
color was observed upon continued passage of the gas. Free rosin 
was precipitated from this liquid by adding an excess of strong hydro- 
chloric acid, and the free rosin was then filtered off, washed twice 
with distilled water, and dried in contact with air. 

Specimens of each of these rosins were then placed in quartz tubes, 
one specimen of rosin in each quartz tube, and the tubes connected 
with gas trains in such a manner that each separate sample of each 
variety of rosin was surrounded by one of the following gases: * 
Purified nitrogen, from which the admixed oxygen and carbon dioxide had been 

removed by passing the gas through a train consisting of gas wash botiles 
containing a mixed solution of potassium pyrogallate and potassium hydroxide. 

Purified oxygen, from which the admixed carbon dioxide had been removed by 
pasing the gas through a train of gas wash bottles containing potassium hy 
drcxide solution. 

Purified carbon dioxide, from which the admixed oxygen had been removed by 
passing the gas through a train of wash bottles containing cuprous chloride 
(CuCl) dissolved in hydrochloric acid, to which solution metallic copper had 
been added. The acid fumes were caught in a saturated solution of sodium 
bicarbonate, thereby preventing removal of carbon dioxide by the train itsel. 

% Measurements by W. W. Coblentz, Bureau of Standards. 


*% The composition of these artificial atmospheres is based on the composition of the air, as has been pre 
viously mentioned, with reasons for choice, in the Introduction. 
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Purified air, which had been washed free from dust by means of a train of bottles 
containing water, after which the excess moisture was removed by means of a 
' train of calcium chloride towers. 
' Air with carbon dioxide removed. ‘This was brought about by passing the air 
' through a train of wash bottles containing a solution of potassium hydroxide. 
Before coming into contact with the rosin samples, each gas was 
‘dried to such a degree that it contained 0.14 to 0.36 mg of residual 
water-vapor per liter at 25° C. 
Duplicate samples of these rosins were sealed in glass tubes, the 
air in which had been displaced by each of the purified gases just 
described. These sealed tubes were preserved in total darkness. In 
} this way there was always on hand a set of reference samples repre- 
senting the color of the rosin in the dark at the start of the experiment, 
Fone sample for every gas used. 
| The test samples of rosin, bleached and unbleached, each in its 
separate quartz tube and each surrounded by its own purified and 
uniformly dried gas—nitrogen, oxygen, carbon dioxide, air, and air 
less carbon dioxide, as the case might be—were then exposed to the 
action of the adopted source of light (q. v.) for a total period of 288 
hours. At the end of that elapsed time the test samples were removed 
from the source of light and compared with the reference samples. 
EF Results are collected and shown in Table 1. 


TaBLe 1.—EHffect of light upon rosin 12 





Sample number and description 
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‘olor is expressed by a comparative numerical scale: 1, darkest, to 5, lightest. Before exposure to light, 
riginal sample was lighter in color than ‘‘5.’’ 
\t the end of each 24-hour period each set of tubes, (a), (6), (c), and (d), were compared with the corre- 
ling original sample and with each other. The numbers in this table represent the results of these 
miparisons, and are not comparable to preceding or succeeding observations. 


2. THE EFFECT OF LIGHT UPON RESINATES OF SODIUM, 
ALUMINUM, AND IRON 


The compounds of rosin which are usually met with in paper- 
making processes are the ‘‘resinates”’ of sodium, aluminum, ferrous 
and ferric iron. In order to protect the compounds used in these 
‘xperiments from the action of light, the whole process of preparing 
these resinates was conducted in a photographic dark room. 
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These resinates of aluminum and of ferrous and ferric iron wer 
made by transposing sodium resinate* with the sulphates of the 
metals just mentioned. Aqueous solutions, using triple-distille 
water, were used in these reactions. The method of operation was as 
follows: The sodium resinate and one of the three solutions of metallic 
sulphates were each placed in separate 250 ml Erlenmeyer flasks, 
previously swept free of air by a stream of purified nitrogen. The 
two flasks were connected with each other in the same manner as , 
train of gas wash bottles, so that the liquid contents of each could 
be agitated by bubbling, transferred by blowing the contents of one 
flask into the other, and mixed by bubbling. The gas used for these 
purposes was the purified nitrogen just described. After blowing 
this stream of purified nitrogen through the sy stem for five minutes, 
the metallic sulphate solution was blown over into the sodium resinate 
solution, and mixed with it by bubbling. The precipitated resinates 
were washed in this manner with distilled water, filtered by suction, 
in an atmosphere of purified nitrogen gas, and dried in a desiccator 
filled with purified nitrogen. All these operations were conducted in 
a photographic dark room. Samples of these resinates were sealed 
in pyrex glass tubes, which contained an atmosphere of purified nitro- 
gen, and stored in the dark for use. Two tubes of each resinate wer 
prepared. One set of these tubes was placed in the light of the carbon 
arc lamp, while the other set was preserved in total darkness as 4 
reference sample. The colors of these compounds before exposure 
were: 

Aluminum resinate: White to pale yellow. 

Sodium resinate: White to pale yellow. 

Ferrous resinate: Light tan. 

Ferric resinate: Yellow-brown to red brown. 
Those resinates which were exposed to the adopted source of light for 
experimental purposes, were exposed for a total period of 25 ‘hours. 
The results obtained are shown in Table 2 


TABLE 2 — of light on resinates of sodium, aluminum, and iron 


Time of exposure (in hours) 





Substance exposed 


| 





Sodium resinate from un- | (¢).__-- 
bleached rosin. 


bleached rosin. 
Aluminum resinate from | (¢)-_. ae en i i 
unbleached rosin. | | 
Aluminum resinate from | (®)...---.---|} ee | aah sen Coen adel 
bleached rosin. | 


| 
| 
| 
Sodium resinate from | (*) ‘ : _) See | 
| 
| 


x (*) 
| 


Ferrous resinate from un- | No change_-| Slight darkening-_| Increased dark- Distinctly | (>) 
bleached rosin. | ening. darker. | 
Ferrous resinate from Noticeable dark- | Very dark color. a --| (*) 
bleached rosin. i | 
Ferric resinate from un- ’ aaa er eeey (4) 2 | ee 
bleached rosin. 
Ferric resinate mem 2 0)... cca a a Se eee. | © ee ig SMO Be | 
bleached rosin. 

| 





* No apparent color changes + No further change. 
” Prepared as per section 1 (d) 
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» With the exception of the rosins, the chemicals used—sodium 
» hydroxide, sodium carbonate, and the sulphates of aluminum, ferric 
' and ferrous iron—were of the ‘‘C. P.” variety. Since iron salts were 


to be added to various compounds which contained no iron during 


: the course of the experimental work, the first three of these chemicals, 


as well as the rosins and solvent acids, were tested for and found 


- free from iron.” 


The water used in making the solutions used in this series of experi- 


| ments was triple distilled water, which had been boiled for six hours, 
during which time a stream of purified nitrogen was bubbled through 


it, for the purpose of sweeping out dissolved oxygen, since ferrous 
iron was to be put into solution. This water was kept in stoppered 
flasks until ready for use. 


| 3, THE EFFECT OF LIGHT UPON ROSIN SIZE CONTAMINATED WITH 


FERROUS AND FERRIC IRON” 


Since Schoeller *® and others*! have claimed that the rosin size 
used in making papers is responsible for the discoloration of 


| papers, and even their subsequent deterioration, and since it has 


been claimed * that small quantities of iron hastened the discolora- 


' tion and deterioration of both size and paper, it was thought advisable 


to study the behavior of rosin sizes contaminated with known small 
quantities of ferrous and ferric iron. The ‘‘size milk” used was 
obtained from the Bureau of Standards paper mill, and contained 
30 per cent of “‘free’” (colloidal) rosin. Samples of this size milk 
were contaminated with solutions of ferrous and ferric iron in con- 
centrations ranging from 1 part of iron in 100,000 to a concentration 
of 1 part of iron in 1,000,000,000. These contaminated samples of 
“size milk” were prepared in duplicate, one set being sealed in 
quartz tubes in air for exposure to the action of light, and the other 
set sealed in hard glass tubes in air and preserved in total darkness, 
as reference samples. The results of a run of 192 light hours are 
shown in Table 3. 





* For Fe++: Potassium ferricyanide, in presence of HCl. \q., ron: y 7 Che 3 
For Fet++: Potassium ferrocyanide, in presence of HCI. }(See Koningh, Z. Angew. Chem., 1898) p. 463. 
, bg Fet++; Potassium thiocyanate. (See Wagner, Zeit. Anal. Chem., 20, p. 350; 1881; also Natanson, Ann 

30, p. 246; 1864. 

By these means a detection of 1 part of iron per 500,000 to 1,600,000 parts of water is claimed by the above 
authors. The chemicals were tested in 10g samples, as above, with negative results; the rosins, after 
decomposition of 40g samples by hot concentrated sulphuric acid, also yielded negative results. 

* In order to avoid ambiguity of meaning and awkwardness of expression, “‘ ferrous and ferric iron” as 
used in this discussion means ‘‘taken singly and together,’’ and “together” means ‘‘equimolar concen- 
trations of each.” 

® See footnote 7, p. 820. 

‘| See footnote 3, p. 819, and 13, p. 820. 

* See footnote 3, p. 819, and 5 and 13, p. 820. 
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TaBLe 3.—Effect of light on rosin size contaminated with ferrous and ferric iron: 
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1! At the end of each 24-hour period each set of tubes was compared with the corresponding original sample 
and with each other. The numbers in this table represent the results of these comparisons, and are not 
comparable to preceding or succeeding observations. 

2 No apparent color change. 

3 In equimolar concentrations. 


Note.—Color is expressed by a comparative numerical scale: 1 (darkest) to 5 (lightest). Before ex- 
posure to light, the original sample was lighter in color than ‘‘5.’’ 


III. DISCUSSION OF EXPERIMENTAL RESULTS 
1. EFFECT OF LIGHT UPON ROSIN 


All of the samples of rosin which were exposed to the action of the 
adopted source of light darkened in color to a considerable degree, 
regardless of the artificial gas atmosphere used. However, the 
bleached and unbleached specimens differed as to comparative 
behavior with respect to light in every artificial atmosphere used, 
with the exception of air and, in a single instance, carbon dioxide. 
In the case of the unbleached samples, maximum darkening occurred 
in those tubes which contained oxygen or air. Some color changes 
were observed in the absence of oxygen, indicating that the unbleached 
rosin can be affected by light alone. The increased effect in the 
presence of oxygen would indicate a photochemical oxidation. On 
the other hand, the bleached rosins darkened most in those tubes 
which contained nitrogen or air, which would tend to show that 
oxygen is not essential to a considerable change in bleached rosins 
upon exposure to light. Abietic acid, found in ordinary rosins, con- 
tains one double carbon bond, and, therefore, may be expected to 
add two atoms of chlorine when subjected to bleaching as described. 
Such chlorinated compounds are, in general, unstable to light.” 
These observations are at variance with the recommendations of 
Schoeller ** and Haas as to the use of ‘‘bleached” rosin for the 
manufacture of paper size. The fact that rosin, bleached as recom- 
mended, was observed to darken in color under the influence of light 
eliminates it from consideration as a substitute for unbleached rosin 
on the basis of comparative light sensitivity. 





*3 Plotnikow, Lehrbuch der Photochemie, p. 547; 1919. 
4 See footnotes 7 and 8, p. 820. 
35 See footnote 12, p. 820, 
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9. EFFECT OF LIGHT UPON RESINATES OF SODIUM, ALUMINUM, 
is AND IRON 


Ferrous resinate was the only one of these compounds which yielded 
a visible color change upon an exposure of 25 hours. Regardless of 
whether bleached or unbleached rosin was used, a distinct change in 
color was observed after an exposure of two hours. This effect is 
probably due to a photochemical change in valence of ferrous to 
ferric iron. 


3. EFFECT OF LIGHT UPON ROSIN SIZE, CONTAMINATED WITH 
FERROUS AND FERRIC IRON, SINGLY AND TOGETHER * 


The samples of ‘‘size milk” containing ferrous iron alone all dark- 
ened in color upon exposure to light, exhibiting a maximum darkening 
in the samples containing 1 part of ferrous iron per 1,000,000 parts of 
“size milk.”? Apparently this is a critical concentration. None of 
the samples to which ferric iron alone had been added exhibited any 
change in color. The evidence indicates that the ferrous iron cata- 
lyzes the action of light upon the rosin size examined. Catalysis of 
photochemical reactions by traces of ferrous iron is a well-known 
phenomenon.*” The color of this intimate mixture of ferrous resinate 
and rosin after exposure to light is much darker than is the color of 
rosin similarly exposed. 


IV. SUMMARY 


As a result of this study of the effect of light upon rosin and rosin 
size, the following conclusions may be drawn: 

1. Rosin is light sensitive, and darkens in color even in an atmos- 
phere of an inert gas. 

2. Both bleached and unbleached rosin show approximately the 
same degree of change in color. 

3. Ferrous resinate is light sensitive, and the color change it under- 
goes upon exposure to light is probably due to photochemical change 
in valence of ferrous iron to ferric iron. Ferric resinate, aluminum 
resinate, and sodium resinate are not light sensitive. 

4. An intimate mixture of ferrous resinate and rosin, as found in 
rosin size contaminated with quantities of ferrous iron smaller than 
those usually found in industrial waters, is light sensitive to a greater 
degree than rosin alone. The evidence indicates that the ferrous iron 
catalyzes the action of light upon the rosin size examined. 

5. Since rosin-sized paper contains rosin, ferrous resinate, and rosin 
size contaminated with ferrous iron, and since these substances have 
been shown to be sensitive to light, some relationship must exist 
between these facts and the discoloration (‘‘yellowing’’) of papers 
upon exposure to light. 

6. This article reports the first of a series of investigations as to the 
deteriorative agencies and their effects upon the components of paper, 


WasHINGTON, February 10, 1931. 





% See footnote 22, p. 821. 
* Plotnikow, Lehrbach der Photochemie, p. 320; 1919. 
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DIMENSIONAL CHANGES IN THE MANUFACTURE OF 
ELECTROTYPES 


By N. Bekkedahl ! and W. Blum ? 


ABSTRACT 


This investigation was conducted by cooperation between the International 
Association of Electrotypers and the National Bureau of Standards. Electro- 
types made under typical but carefully controlled conditions were measured after 
each step in the process. The results showed that each operation is likely to 
produce some change in dimensions. The largest and most variable changes are 
produced in curving. By controlling all conditions it is possible to produce plates 
that will register with each other within about 0.05 per cent, though they will 
usually not have the same dimensions as the originals, 
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I. INTRODUCTION 


Electrotyping is extensively used for the reproduction of printing 
plates. In this process the details of the original, whether type or 
half tone, are faithfully reproduced. However, it has long been 
known that, unless special precautions are observed, the electrotypes 
are not exactly the same size as the originals. This difference, which 
in curved plates may be over 1 per cent, is usually not significant if 
the printing is done in one color, as then only one plate is involved.® 
In multicolor printing, however, two or more plates are employed; 
that is, one for each color used. For high-grade printing, these plates 
must have the same corresponding dimensions in order to register 


' Research associate of the International Association of Electrotypers. 
* Chemist, National Bureau of Standards. 
In printing certain charts, scale diagrams, or cross-section paper, accurate reproduction may be required 
even for 1-color plates. 
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properly. Even in this case it is not essential that the different plates 
should be exactly the same size as the originals, but they must, bp 
alike in size. The difficulties sometimes encountered by commercia| 
electrotypers in meeting this latter requirement led the Internationa 
Association of Electrotypers to select this subject for study by thei 
research associate at the National Bureau of Standards. 

This paper is a report of an investigation of this subject, in which the 
dimensional changes produced in each step of the electrotyping process 
were measured. In addition, the effects of changes in electrotyping 
conditions upon the changes in dimension were studied in order to 
determine how closely such conditions must be controlled to insur 
any given degree of uniformity in the size of the finished electrotypes, 


II. OUTLINE OF THE ELECTROTYPING PROCESS 


Electrotypes are made by depositing metal either upon wax molds 
or lead eal, of which the former are more extensively used. In wax 
molding, an impression of the original form is made by pressing the 
latter into a layer of wax that is attached to a flat metal plate. Graph- 
ite is then applied to render the wax surface electrically conducting, 
and metal is deposited on the graphite, forming a thin ‘shell.’ 
Sometimes the graphite surface is treated with a copper sulphate 
solution and iron filings, which produce a thin film of copper on the 
graphite. This process, known as “oxidizing,” causes the metal to 


subsequently deposit almost at once over the entire surface, instead o! 
““srowing”’ over from the point where metallic contact is made with 
the graphite. In copper electrotypes, the shell consists entirely of 
copper, while in nickel electrotypes the first layer is nickel and the 


rest of the shell is copper. The nickel is used especially for large 
editions and for high-grade printing. At present chromium is often 
applied to finished nickel or copper electrotypes to increase their weat 
resistance. 

A lead mold is made by pressing the original form into a thin sheet 
of lead. The lead mold is cleaned and treated with a chromate soli: 
tion to permit separation of the subsequently deposited shell, which 
is generally of nickel and copper. Lead molds are used especially 
for half tones. Type-metal forms are not usually reproduced by lead 
molding, as the high pressure used is likely to deform the type. |! 
is a general practice, however, to make lead molds from specially pre- 
pared electrotype patterns of a type subject. 

The back of the shell is treated with a solution of zine chloride which 
serves as a flux; after which “tin foil,” consisting of about 35 per cent 
tin and 65 per cent lead, is laid on the surface. The shell is then 
placed face down on a flat cast-iron pan which is heated, usually by 
floating it upon molten lead or electrotype metal. When the tm 
foil has melted, molten electrotype metal is poured over the shell to 
the desired thickness. After the resultant plate is cool it is cleaned, 
trimmed, shaved on the back to approximately the desired thickness, 
“finished”? to produce a plane printing surface, shaved to the exac! 
thickness, and for flat presses is beveled or mounted on wood. _ If the 
plate is to be used on a rotary press, it is curved to fit the cylinder ©: 
that press. Any large nonprinting areas on the plate are usually 
“routed out,’ in order to prevent them from taking up ink and 
yielding it to the paper. 
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III. OUTLINE OF THIS INVESTIGATION 


| General observations indicated that each of the major steps in the 
above process might affect the dimensions, namely, the molding, 
' depositing, casting, finishing, and curving. In a preliminary study, 
' plates were made under typical commercial conditions, and the dimen- 
- sions were measured before and after each of the above operations. 
’ These tests indicated the magnitudes of the changes produced in 
each step. The process was then carried out under more carefully 
controlled conditions, each of which was varied over definite limits, 
to determine the corresponding variations in the dimensional changes. 

This investigation was facilitated by the active cooperation of the 
Potomac Electrotype Co., of Washington, D. C., and the Govern- 
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FicgurE 1.—Design of “‘plain”’ pattern plate, showing lines used as gage marks 
The plate was 11 by 14 inches (28 by 26 cm). 


ment Printing Office. In general, the molds were made in one of 
these plants, the shells were deposited at the National Bureau of 
Standards (in order to measure and control the current used on each 
plate), and the casting, finishing, and curving were done in the plants. 
lhe data, therefore, represent closely the results of actual electro- 
typing conditions. 

In this study two pattern plates were used. ‘The first, designated 


ée 


as the “plain” form, consisted of a polished copper plate about 11 
by 14 inches (28 by 36 cm). Upon one side were ruled ‘ fine lines, 
about 0.001 inch (or 0.025 mm) in width which formed a large rec- 
tangle divided into eight small ones, as shown in Figure 1. These 
lines served as gage marks, between which the distances were meas- 
ured after each step in the process. 





‘ Through the courtesy of the United States Coast and Geodetic Survey. 





Bureau of Standards Journal of Research (Vol 


The other pattern, designated as the “‘type”’ form, was an electro. 
type plate, the surface of which consisted mostly of type, with some 
“rules.”” This plate was about 9 by 12 inches (23 by 30cm). Meas. 
urements were made between certain well-defined points on this plate 
and its reproductions, 

These two patterns were selected to represent the extremes likely 
to be encountered in printing plates. Halftone plates are in very 
low relief, and hence approach the plain form. Commercial plates 
that include both type and halftone may reasonably be expected to 
yield results that are intermediate between those obtained with the 
two patterns used in this investigation. 

No difficulty was experienced in measuring the distances between 
given points on flat plates by means of a graduated steel scale. It is 
however more difficult to measure accurately the distance between 
two points on a cylinder, especially if between these points a part of 
the cylindrical surface has been removed by “‘routing.” If a flexible 
steel tape were used, it would tend to span the depression as a chord, 
instead of as an are of the specified radius, and would hence not 
yield the true distance on the cylindrical surface, which latter is the 
actual printing surface. This difficulty was overcome by using the 
special gage shown in Figure 2, which was designed by one of th 
authors (W.B.). ‘This gage was made from a polished flat steel plate 
about 15 by 8 by *i5 inches (38 by 20 by 0.5cem). The arc had a radius 
of 6.250 + 0.002 inch (158.8+0.05 mm). This are and one straight 
edge were graduated in linear units, the smallest division being 0.01 
inch (equal to 0.25 mm). With an optician’s lens, readings could 
readily be made to +0.002 inch (0.05 mm). The gage was used for 
measuring both the flat and curved plates. 

This particular gage was designed for the curved plates used in 
this investigation, which were made to fit the cylinder of a 3-color 
Sterling press. In order to make similar measurements on plates 
curved for different size cylinders, a separate gage would be required 
for each radius of curvature of the printing surface, which is one-hali 
the so-called ‘‘bearer diameter.” In a large electrotyping plant, this 
might necessitate the use of a large number of such gages if the di- 
mensions of the curved plates were to be measured. 

The dimensional changes are expressed in terms of (a) changes 
produced by a given step in the process, and (b) the net change 
from the original, produced by all operations performed up to and 
including that step. The latter values, or more strictly the final 
dimensions, are significant in the use of the plates; while the former 
show how variations in the latter may arise. 

The register of two or more printing plates depends upon the 
actual differences in distance between corresponding points on the 
plates. Experience indicates that for fine printing; for example, of 
color process half tones, the extreme difference of corresponding 
points in the final impressions should be not more than about 0.003 
inch (0.08 mm). To secure this accuracy of register, it is necessary 
that the dimensions of the plates should correspond only within 
about 0.005 inch (0.13 mm); as the effort is always made to have the 
plates register exactly at their centers, which causes half of the error to 
appear on each end. For practical purposes we may then consider 
that the object of this investigation was to define conditions that would 
insure an equality of length within + 0.005 inch (0.13 mm). 
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Such a variation in length would, however, represent entirely 
different percentage discrepancies, according to the dimensions of 
the particular plate; for example, 0.1 per cent on a 5-inch (12.7 cm) 
plate, or 0.02 per cent on a 25-inch (63.5 cm) plate. Evidently, 
therefore, a measurement of the percentage variation is not in itself 
a criterion of the accuracy of register. As printing plates, including 
those made in this study, vary in dimensions, the only practicable 
way to express general conclusions is in terms of percentage changes. 
For any given size of plate the corresponding actual differences 
may then be calculated. Thus for a 10-inch (25 em) plate, which 
is probably about the average length used commercially, each one- 
hundredth of 1 per cent is equal to 0.001 inch (or 0.025 mm). The 
general conclusions of this paper are based on a plate of about that 
length. 

All the measurements reported in this paper were the averages of 
at least two plates, and in some cases of 10 or more. Except in a 
few cases that will be mentioned, the individual results were con- 
cordant to within about +0.03 per cent. The changes were gener- 
ally about equal in both directions, except in curving, when the 
significant changes occur only in the direction of curvature. In 
lead molding, changes in the two directions may be slightly different 
if the lead overlaps the pattern in one direction and not in the other. 
In the measurements here recorded, the lead projected beyond the 
pattern on all four sides, which is the common commercial practice. 
A positive sign attached to a change, for example, + 0.03 per cent, 
indicates an increase in length (‘‘stretch’’) and a negative sign a 
decrease (‘“‘shrinkage’’). All measurements were made at room 
temperature, which was approximately 75° F. (24° C.). 


IV. RESULTS OF EXPERIMENTS 


The conditions used in each operation and the dimensional changes 
thereby produced, may conveniently be considered under each of the 
principal steps involved. 

1. MOLDING 


(a) WAX MOLDING 


The cases used for wax molding were prepared by pouring 
molten wax consisting largely of ozokerite, to a depth of about 9.15 
inch (4 mm) on a flat plate made of electrotype metal (containing 
about 93 per cent lead, 4 per cent antimony, and 3 per cent tin). 
The surface of the wax was shaved to make it plane, and brushed 
lightly with molding graphite to prevent the form from sticking. 
Molding was done in a hydraulic press at the slightly elevated tem- 
perature that is commonly used, and also at room temperature. The 
results in Table 1 show that in neither case was there any measurable 
difference in the dimensions of the original form and the wax mold. 


TaBLE 1.—Dimensional changes in molding 
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This result is somewhat surprising in view of the fact that the 
molding wax has a much higher coefficient of expansion than that of 
copper. If the wax were free to expand and contract laterally the 
wax mold made at a high temperature would become smaller thar 
the copper pattern when both were again at room temperature. 
The expansion coefficients of some of the materials used in electro. 
typing are shown in Table 2. 


TABLE 2.—A pproximate linear coefficients of thermal expansion 


( Parts per 
Material million per 
degree C, 





Aluminum. Bd 
PT oe 


(0 to 50° C.) if 
— ate 300° C.) 18 





Zinc 
Molding wax (ozokerite base) - - 








The values in this table for the metals are derived from those values 
in the literature that represent approximately the temperature ranges 
involved in this investigation. Aluminum 1s included because it or 
its alloys are sometimes used in the case on which the wax is poured. 
Zinc is included because zine plates often serve as patterns for electro- 
typing. No values were found for electrotype metal. The value 
found ° for a type metal containing 85 per cent lead, 12 per cent 
antimony, and 3 per cent tin is 25 parts per million per degree centi- 
grade. The available data indicate that the expansion ‘coefficient 
of ordinary electrotype metal, which contains 93 per cent of lead, is 
intermediate between this value and that of pure lead, or about 
28 parts per million per degree centigrade. The expansion coefli- 
cient of electrotypers’ wax varies with its composition. The value 
given in the table is the approximate average of a few samples on 
which purely preliminary measurements were made in connection 
with this investigation. 

If then the wax expanded freely when heated from about 75° to 
120° F. (24° to 49° C.), its length would increase by about 25x 
0.000250 = 0.00625, or 0.63 per cent. The copper form would expand 
in this same temperature range, by about 25 X 0.000016 = 0.000400, 
or 0.04 per cent. Hence, we would expect that on cooling, the 
wax mold would contract more than the copper and at 75° F. (24° C 
would be smaller than the original copper by the difference in dheee 
two expansions (or contractions); that is, by 0.63 minus 0.04 per 
cent, or practically 0.6 per cent. As no such difference was observed, 
it is evident that the lay er of wax that is attached to the metal base, 
in this case type metal, is not free to contract or expand laterally, 
and that it retains very closely the size of the underlying plate. 
If this is true, the only 5 Monee that would be measured after molding 
would be that due to the difference in expansion of the copper pattern 
and the type metal case. This would be equal to about 25 (0.000028 - 





5 Geo. F. Taylor, Phys. Review, 26, p. 841; 1925. 
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9.000016) =25 (0.000012) =0.000300, or 0.03 per cent. This dif- 
ference, which would result in —0.03 per cent change, is just about 
at the limit of observation, especially on wax molds. If, however, 
the copper pattern is expanded slightly by pressure, even at the 
moderate pressure used in wax molding, this will tend to compensate 
for the tendency of the mold to contract. At least the results show 
clearly that expansion of the molding wax does not lead to an appre- 
ciable change in lateral dimensions, and that differences in the 
composition or expansivity of the wax are not likely to affect the 
dimensions of the mold. 
(b) LEAD MOLDING 


As lead molding is done at ordinary temperature, there is no oppor- 
tunity for dimensional changes to be caused by differences in thermal 
expansivity. Nevertheless it was found, as shown in Table 1, that 
the lead molds are slightly larger than the originals. This difference 
is undoubtedly due to the fact that at the higher pressures used in 
lead molding, the pattern expands slightly. If this expansion is with- 
in the elastic limit of the pattern metal, in this case copper, the pat- 
tern returns to its original length when the pressure is removed, while 
the lead, which has a very low elastic limit, retains the slightly larger 
dimension. 

That the pattern plates do expand slightly in lead molding, even in 
extreme cases to the extent of a permanent change, is indicated by 
the fact that the plain pattern was 0.02 per cent larger in both length 
and width after about 100 impressions had been made from it. The 
type pattern was about 0.07 per cent longer and 0.02 per cent wider 
after about 50 impressions had been made. The greater expansion in 


length of the type plate was probably due to the fact that the type 
lines ran across the plate. Hence, the lateral pressure perpendicular 
to these lines was greater than that parallel with them. These very 
small changes in the size of pattern plates as a result of repeated 
molding have been reported by electrotypers. 


2. DEPOSITION 


It is well known that all electrodeposited metals are in a state of 
strain and tend to contract during their deposition. This contraction 
is most marked with the harder metals, and with those, such as nickel, 
that absorb much hydrogen. For a given metal the tendency;to con- 
tract is greatest at a low temperature and a high current density, and 
increases with the thickness of deposit. 

These well-known principles are illustrated by the data shown in 
Table 3. In these experiments the copper was deposited from a solu- 
tion containing 250 g/l (33 oz./gal.) of copper sulphate (CuSO,.5H,0), 
75 g/l (10 avoir. oz./gal.) of sulphuric acid (H,SO,) and 1.0 g/l (0.13 
oz./gal.) of phenol, added as phenolsulphonic acid. Preliminary 
observations showed that the use of phenol as an addition agent pro- 
duced no measurable change in the dimensions of the shells. (As 
will be noted later, it may affect the changes that occur in curving.) 
lhe nickel was deposited from a solution containing 70 g/I (9 oz./gal.) 
of nickel sulphate (NiSO,.7H,O) and 5 g/l (0.7 0z./gal.) of ammonium 
chloride (NH,Cl) at a pH of about 6.4 (with brom cresol purple, 
uncorrected for salt error).® 


* W. Blum and N. Bekkedahl, Trans. Am. Electrochem. Soc., 56, p. 291; 1929. 
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TaBLE 3.—Dimensional changes in deposition 
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In most experiments the shells were smaller than the molds, and 
the contraction of the nickel-copper shells was always greater than of 
the copper alone. With the latter, the contraction was significant 
only with shells deposited at a relatively low temperature (experi- 
ment 5). The absence of a measurable change when the copper was 
deposited at a high temperature, even at a high-current density, was 
propably due to the fact that, as previously explained, the wax mold 
on a type-metal base was slightly larger (about 0.03 per cent) at the 
elevated temperature. Hence, the small contraction produced by 
deposition may have just compensated for this tendency to expand. 
Strictly speaking, the contraction caused by the deposition itself on 
either the wax or lead molds at the higher temperature was, there- 
fore, about 0.03 per cent greater than that shown in Table 3.  Prac- 
tically, however, the actual changes in dimensions of the shells are 
significant. 

The shrinkage produced in depositing nickel-copper shells was 
relatively large (from 0.12 to 0.15 per cent) but fairly uniform, regard- 
less of whether a plain or type form was used, or whether the shells 
varied in thickness by as much as 100 per cent. In the latter case, 
the contraction was slightly greater, but not sufficiently different to 
affect the register. 

The last column in Table 3 shows that the changes in molding and 
in depositing tend to compensate each other. Hence, the dimen- 
sions of the shells are not muc +h different from those of the patterns. 
As will be seen, the changes in casting and curving are relatively 
large. If these ‘operations could be eliminated, there would be little 
difficulty in obtaining plates of uniform dimensions. ‘This elimina- 
tion is not practice able in regular electrotyping as now practiced, but 
it has been accomplished by making the entire printing plate of 
deposited metal; for example, of iron, at the Bureau of Engraving 
and Printing,’ or of copper as is being done for special purposes in & 
few elec ‘trotyping plants. 








. Thomas and W. Blum, " Trans. Am. Electrochem. Soc., 57, p. 59; 1930. 
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TaBLE 4.—Dimensional changes in casting 





| 


8 Change— 
composition | | 


Electrotype metal Shell thickness 
A | 
Experi- ——---— — ee a Casting Rate of |—— 
‘ Pattern | | rature| cooli | 
\temperature) cooling | 1, From 
Sb | Nickel | Copper cast- | orig- 
| | | ing j|inal! 








P. My IP. ct. .|P.a.| Inch | mm | Inch LO | ; | Piet. 1 P.c. 
4 3 |....-.|------0.006 | 0.15} 635} 335 | Fast__..|—0.18 |—0.15 

93 4 006 | .15 335 | 335 : — 15 |— .2: 

93 4 3 | .001 | .025 | .006| .15 oH 20 |...do. - 114] 

93 | 4 5 | .006 -15| 635 335 | Slow — .16 |—- 


93 | 4 3} .002 | .05 012 

88 | ¢ Be | 025 | .006 

93 d 9 Fes . 006 

93 ‘ 3 | .001 | .025 | .006 

oa. - 93 3 | .00 . 02! . 006 
_do. aha 93 | 4 Si 025 | . 006 


93 012 


ae : é 
_do. a 8S { Si« 1 | .025 | .006 


! These values were obtained by assuming lead ial see copper deposition at 105° F. (40° C.), and nickel 
deposition at 85 F. (30° C.). 


3. CASTING 


The data in Table 4 show clearly that the casting process always 
causes shrinkage of the plates. ‘The most reasonable explanation is 
that the solid electrotype metal has a higher expansivity th an the 
copper shell. Hence, when the cast plate cools, the type metal tends 
to contract the copper. This effect can not start until the electro- 
type mets al has solidified, as above the melting point of the type metal 

e latter and the copper are free to expand and contract independ- 


coil, As the electrotype metal freezes at about 527° F. (275° C.) 
(see fig. 3) we may assume that the di ifferential contraction of the 18 
type _ al and copper is significant from about 275° C. to 25° C. 

that is, for about 250 centigrade degrees or 450 Fahrenheit degrees. 


The independent contractions in this range would be about as follows: 


Copper = 250 X 0.000018 = 0.0045 =0.45 per cent 
Type metal = 250 x 0.000028 =0.007 =0.7 per cent 


From these figures we might expect that the cast plate would be 
about 0.25 per cent shorter than the shell, instead of only about 0.15 
per cent as shown in Table 4. The difference is no doubt due to the 
fact that the copper (or copper-nickel) shell has too high a compres- 
sive strength to be fully contracted by the electrotype metal. Hence, 
the latter must either slip or be stretched beyond its elastic limit 
during the cooling. The net result is a shortening of the plate by an 
amount that is less than the difference between the two coefficients. 

If the above explanation is correct, the casting temperature should 
have no effect on the shrinkage, as the latter can not start until the 
metal solidifies. This prediction is borne out by the data in Table 4. 
(Compare experiment 18 with 20 and 25 with 26.) 

Differences in the rate of cooling produce no marked differences in 
the dimensions. In the “slow” process, the cast plates were allowed 
to stand on the pans without any artificial cooling. In the “fast” 
cooling, cold air was blown against the bottom of the pans. The 
different rates of cooling are shown by the ‘cooling curves” in 
Figure 3. These were made by imbedding ‘‘chromel-alumel’’ 
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thermocouples in the type metal as the latter was poured on the pans, 
and reading the temperatures at regular time intervals by means of q 
millivoltmeter. In all cases it was observed that the “break” jp 
the curve; that is, the freezing point, occurred at about 527° Ff, 
(275° C.), even with type metals of different compositions. It is 
recognized that an accurate measurement of the melting or freezing 
point was not obtained under these conditions, as the metal was not 
stirred while cooling. The results represent, however, the approxi- 
mate cooling rates of the type metal under the conditions used. 
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Fic. 3.—Cooling curves for electrotype metal in the casting pans 


The ‘‘break”’ in each curve represents the freezing or solidification of the metal. 

A, 88 per cent lead, 9 per cent tin, and 3 per cent antimony. Cooled rapidly. 

B, 93 per cent lead, 4 per cent tin, and 3 per cent antimony. Cooled rapidly. 

C, 93 per cent lead, 4 per cent tin, and 3 per cent antimony. Cooled slowly. 

(For convenience the initial readings are spaced at 5-minute intervals. To compare B or C with A, 
5 minutes should be deducted from the values on B, and 10 minutes from those on C.) 


An increase in the tin content of the type metal produced a some- 
what smaller contraction of the plain plates (compare experiment 
21 with 23), but no similar difference with the type pattern (com- 
pare experiment 27 with 29). The higher tin content presumably 
reduces slightly the expansivity of the electrotype metal, which 
would result in a smaller contraction of the cast plate. The larger 
contraction with the type pattern may be due to the fact that the 
latter presents a rough surface and, therefore, retards any slippage, 
and also that the corrugated surface of the shell is more easily con- 
tracted than is a plain surface. Any warping of the plates is likely 
to — their measured dimensions less than those when the plate 
is plane. 
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The other factor that caused a variation in the contraction was the 
thickness of the shell. The smaller contraction of a thick shell (com- 
pare experiment 21 with 22) is to be expected since the thicker copper 
and nickel would offer more resistance to shrinkage. Here again the 
effect is negligible with the type pattern, probably for the reasons 
stated above. 

In general, it may be concluded that the casting process produces 
a contraction of about 0.15 per cent, which value is not appreciably 
affected by differences in the temperature of casting or rate of cooling; 
but it may be changed slightly, especially on plain plates, by a change 
in the composition of the type metal or in the shell thickness. As the 
deposited shells have nearly the same dimensions as the originals, the 
contraction in casting always results in a cast plate that is smaller 
than the originals. 


4. FINISHING 


The cast plates are not likely to have a plane printing surface, as 
both the deposition and the casting tend to produce distortions. 
The purpose of the finishing process is to straighten the plate so that 
the entire printing surface is in one plane. This operation generally 
tends to cause an increase in the dimensions of the plate. The extent 
of this stretch depends upon the amount of finishing required and the 
method used to eliminate defects that may have been present in the 
original patterns or produced in the electrotyping process. It is there- 
fore difficult to obtain data that represent ‘‘average”’ conditions. 
The results in Table 5 are simply typical values that show the probable 
magnitudes of the changes. 


TABLE 5.—Dimensional changes in finishing 





Change— 
Experi- ™ 


mant No Pattern | Shell metal Method of finishing 


|In finish-| From 
| | ing original ! 





| | Per cent | Per cent 

| Copper.........-------| Hand hammering.-.......| -+0.14 —0. 

Pe an eee Nickel-op0pet..........--|-.+«-00..-~.-- : +.09 —.% 
do ol Se is” | | +. 03 ~~. 1 
es ted aaa Nickel-copper-.--.-.-.-...|..--.do.---.. +. 03 —.17 
-----| Copper-.---- Hot pre -.11 


—. 08 


Nickel-copper-_-----. tts dvaw'si , 
—.09 


| Copper..--- : ei... 


+.10 
-+. 08 


satu! AZODMOE hacnd-.c.~0-.-.) Ole DEOBa. 
.--| Nickel-copper-------- eee és 








; : Eee 
See OS eee ee: | — 06 


1 These values were obtained by assuming lead molding, copper deposition at 105° F. (40° C.), and 
nickel deposition at 85° F. (30° C.). 


In the common practice of hammering the back of the plates at the 
points that are too low on the printing surface, the amount of stretch 
is variable. In some cases at least the expansion thereby produced 
is just sufficient to bring the dimensions back to the original. At best, 
however, hand finishing is troublesome and expensive, and in recent 
rs many efforts have been made to perform this process mechani- 
cally, 
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Straightening in a press at ordinary temperature also causes expan- 
sion, the extent of which depends largely on the pressure applied, 
This expansion by pressure is similar to but greater than that pro- 
duced by molding, in which lower pressures are generally used. The 
results obtained indicate that for similar plates and designs and at 
a controlled pressure, the expansion is fairly uniform. Type plates 
expand more readily than plain plates. In all cases the net result is 
to make the plate dimensions approach the originals. 

When the plates are straightened in a press at a temperature of 
approximately 500° F. (260° C.); that is, just below the melting point 
of the type metal, there is a tendency to contract. Apparently in 
this second cooling of the electrotype metal, its high coefficient of 
expansion (or contraction) causes it to produce additional shrinkage 
of the copper. Plates straightened by this method are likely to be 
appreciably shorter than the originals. 

Tests not recorded in the tables showed that increasing the thick- 
ness of the shells had no appreciable effect upon the changes produced 
by finishing with any of the above processes. 





5. CURVING 





If a flat homogeneous plate is bent around a cylindrical mandrel, 
the outer or convex surface tends to expand from tension, and the 
inner or concave surface tends to contract from compression. If we 
assume in the ideal case that the effect of the tension is exactly equal 
to that of compression, the proportional stretch of the outer surface 
is equal to the thickness of the plate divided by the outside diameter 
of curvature. This assumption implies that there is in the middle of 
the plate a ‘“‘neutral plane”’ which neither expands nor contracts. 

The plates used in this investigation, had a thickness of 0.183 inch 
(4.65 mm), and were curved to yield an outside diameter of 12.500 
inches (317.5mm).° According to the above formula, the percentage 
stretch would be 






Thickness X 100 0.183 x 100 
Per cent stretch =—— we? fa. = 1.46 per cent 
Diameter 12.5 








Actually the average stretch for all the plates included in Table 6 is 
1.42 per cent. This close average agreement is, however, largely a 
coincidence, as under different conditions the stretch varied from 1.02 
to 1.65 per cent, for reasons that will be discussed. Even under 
specified conditions the results were not as reproducible as those 
obtained in other operations. Small differences in the tabulated 
results are, therefore, not significant. 
















* As these plates were made for use on a cylinder with a diameter of 12.125 inches ‘308 mm), the desired 
outside diameter would correspond to a plate thickness of 0.1875 inch (4.76 mm). A smaller plate thickness, 
in this case 0.183 inch (4.65 mm), is generally used in order to allow for an ‘‘under lay” of paper in the 


“‘make-ready”’; that is, the process of bringing all parts of the printing surface into the position needed to 
give it the desired “printing value,” 
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TABLE 6.—Dimensional changes in curving 








| Electrotype metal | Shell thickness | Change— 
— tr f tae Hours| Method 
Experl- } | | Rate of | till of | 
ment | Pattern | | cooling | cur- | curv- In |From 
No. Pb Ss Vicke Copper ved ing curv- | orig- 
} ing | inal! 
40 Sheena lta 4 
| 
P.ct.| P.ct.| Inch | mm | Inch | mm | P. et. | Pas 
PES ERS a ne ----|+1.41 |41. 29 
0. 002 |0.05 | .012 ; ‘ . 0: 85 
.O01 | .025 | .006 ou a: a 1 i+1. 30 |+1.13 
001 | .025 | . 006 | J . 33 | .16 
001 | .025 | . 006 ga lowe ened ‘ . 29 +1. 12 


| tH 
93 





.001 | .025 | .006} .1f -do....| 24- |+1. 2% . 05 
| .001 | .025| .006] .1é -do....|  24- +1, 58 |+1. 43 
| 002} .05 | .012] .: ...-| 24+ |...do. +1.41 |+1. 26 
3} .002| .05 | .012] .30|..-do._..| 24+] C --- +1. 35 |+1. 20 
.001 | .025 | .006| .15| Fast.-.| 24+ |-..do...../+1. 5 . 38 





.001 | .025 | .006] .15|_-.do....| 244 wonul ke . 50 
.001 | .025 | .006| .15! Slow ..| 2 Yold ....'!+1. 51 . 36 
.001 | .025 | .006| .15|---do....| 244 _..--| +1. 60 |+1. 45 
: 006} .18:|...do... : old ...-|+1. 49 . 34 
93 ; | .006} .15 }-..do....| 24+ |--- +1. 56 51 

















1 These values were obtained from assuming lead molding, copper deposition at 105° F. (40° C.), nickel 
deposition at 85° F. (30° C.), and cold-press method of finishing. 


An electrotype plate is not homogeneous, and the outer shell is 
much stronger than the electrotype metal that constitutes over 90 per 
cent of the plate. Hence, the neutral plane as above defined is nearer 
to the printing surface, and on curving, the plate stretches like a 
thinner plate; that is, the expansion is less than that calculated as 
above. This is illustrated by the data in Table 6 for plain plates, 
for which the average stretch is only 1.26 per cent. 

On the other hand, the type pattern represents a surface corrugated 
in the direction perpendicular to that of curvature. When such 
plates are bent, the average increase in length is 1.52 per cent; that is, 
slightly greater than that calculated. In bending a type plate the 
increase in length of the printing surface is caused not only by stretch- 
ing the shell, but also by bending the latter where the vertical portions 
meet the depressed horizontal portions. This is illustrated in Figure 
4. The depressions in the face have the same effect as if the diameter 
of the plate were decreased, while the thickness (as defined by the 
printing surface) remains constant. 

This same effect is observed to an even greater degree if certain 
areas of the plate are routed out before the plate is curved. Experi- 
ments not recorded in the table showed that the stretch of a routed 
plate is much greater than of an unrouted plate. This difference is 
due partly to a decrease in the effective diameter as above noted, and 
also to the reduction in the strength of the shell by removing portions 
of it. It is therefore preferable to rout the plates after curving. It 
was found in these tests that this subsequent routing does not sated 
the dimensions of the curved plate. When it is practicable to reduce 
or eliminate the routing by ‘“‘building up” the wax mold, the shell is 
continuous, and the stretch is more nearly uniform than that of a 
plate routed before curvature. 

_ With either kind of plate, the stretch was decreased by any increase 
in strength of the shell such as may be produced by using both nickel 
and copper, by increasing the thickness, or by increasing the strength 
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of the copper; for example, by the use of addition agents. <A few 
tests not recorded in the table showed that plates made from shells 
deposited in a copper solution that contained no phenolsulphonic acid 
stretched slightly more than those here recorded. The smallest 
stretch (1.02 per cent) was obtained with a nickel-copper shell, plain 
form and of double thickness. 

A slight decrease in the stretch of plain plates was caused by ip. 
creasing the tin content of the electrotype metal (compare experiment 
43 with 45), but any such effect with the type pattern was very small 
(compare experiment 46 with 52). As the addition of tin hardens the 
type metal, a decreased stretch would be expected. 

No appreciable difference was observed in the stretch whether the 
cast plates were cooled rapidly or slowly (compare experiments 5! 

















Ficure 4.—Diagram showing types of curved plates 


A. A homogeneous plate. When curved the convex surface stretches approximately in accordance with 
the equation 


Per cent stretch= eee RS. 
er Diameter of curvature 


The neutral plane is in the center. 

B. A plain electrotype plate, with a shell of copper or nickel, backed-up with electrotype metal. The 
stretch is less than that calculated in A. The neutral plane is closer to the outer surface. 

C. An electrotype ofa type pattern. This stretches slightly more than that calculatedin A. The neutral 
plane is closer to the inner surface. 


with 49, and 52 with 50). Similarly it made no difference whether 
the plates were curved soon after casting (within about two hours), 
or after 24 hours or more. 

The cold curving was done at room temperature with the ordinary 
3-roll bender. In this method the plates are rolled to fit the mandrel, 
a hemicylinder which is known in electrotyping as the ‘‘turtle.” 
In hot curving the plate is heated to about 480° F. (249° C.) and is 
then bent directly to fit the appropriate cylinder. In general, the 
hot curving produced a slightly greater stretch than cold curving. 
The difference was not significant with the plain plates (compare 
experiment 42 with 43), but was appreciable with the type plates 
(compare experiments 51 with 52, and 49 with 50). The greater 
stretch in hot curving is probably caused by the fact that at the higher 
temperature the electrotype metal softens relatively more than does 
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the shell. The results with hot curving were somewhat more uniform 
than with cold, and the hot-curved plates retained their curvature 
better than did the cold-curved plates. 

At best it is difficult to produce curved plates of equal dimensions. 
When a given method of curving is used, the results are fairly repro- 
ducible if the shells are of the same thickness, strength, and type of 
pattern and if the type metal is of uniform composition. This con- 
clusion agrees with the common experience that it is much easier to 
make several plates of a series with accurate register if they are pro- 
duced simultaneously than at different times. For example, when a 
single plate is required to replace a spoiled one of a set, it is often 
difficult to make 1t register, because then the conditions may not be 
the same as when the first plates were made. 

As previously noted, it is occasionally desirable to produce elec- 
trotype plates that have exactly the same dimensions as the originals. 
The results described above show that it may be possible, though 
difficult, to produce fiat plates that meet this requirement. Inspec- 
tion of the data in Table 5 shows that the stretch produced by cold 
finishing tends to compensate for the contraction in casting. In 
special cases the finishing operation may be deliberately extended in 
order to restore the original dimensions. 

In the ordinary methods of curving, however, the stretch is usually 
1 per cent or more. If uniform, this stretch is not objectionable for 
regular printing. For “dimensional” printing, however, such plates 
are useless. One way to overcome this difficulty is to make the 
pattern plates shorter in the direction in which they are to be curved. 
This involves a preliminary measurement of the stretch. In another 
ingeneous method, commonly known as the ‘‘nonstretch’’ process, 
electrotype metal is poured upon the face of a finished flat electrotype 
which has been treated with chalk or graphite to prevent the metal 
from sticking permanently. As a first approximation it might be 
assumed that when such a plate is cast so as to have the printing face 
practically in the middle, and is then curved by the usual methods, 
the printing surface would’ become the “neutral plane’’; that is, it 
would neither stretch nor contract during the curving. 

In order to determine the effectiveness of this procedure a few 
experiments were made with the same patterns used above. The 
results showed clearly that if the printing surface is in the center of 
the plate there is still an appreciable stretch. This can be reduced 
by making the added type metal layer thicker than the layer on the 
under side of the shell. For example it was found that if the plate 
(including the shell) had a thickness of 0.094 inch (2.35 mm) and 
type metal was cast on the face to a thickness of 0.166 inch (4.2 mm), 
thus making the total thickness 0.260 inch (6.6 mm) and the plate 
was subsequently curved on the 3-roll bender, the stretch of either a 
plain or type plate was only about 0.15 per cent. When, however, the 
thickness of. the added type metal was increased to make the total 
thickness 0.354 inch (9.0 mm) and the plate was then curved, there 
Was an actual contraction in length of the printing surface. 

_ To eliminate the stretch the shell must evidently be closer to the 
inside than to the outside surface. It is necessary to determine by 
experunent what thickness of type metal must be applied to the 
surface to reduce the stretch to approximately zero. When this is 
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once determined for a plate and shell of specified thickness and 
curvature, fairly uniform results can be obtained. It is readily 
possible thereby to reduce the stretch to about 0.2 per cent, which 
tends to counteract the shrinkage caused in casting. (This compen. 
sation occurs only in the direction of curving, and not in the other 
direction.) At best it is very difficult to produce finished curved 
plates that reproduce the originals within about 0.2 per cent in 
both directions. 
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V. SUMMARY 
The different operations in the electrotyping process produce the 
following changes in dimension and must be controlled accordingly 
to produce plates with uniform dimensions. re 
A 


1. MOLDING 






Wax molds are almost exactly the same size as the originals. 
Lead molds are slightly larger (about 0.05 per cent) than the originals. 
No special precautions are required to produce uniform plates by 
either method of molding, but plates made by the two methods are 
not likely to match each other exactly. 


2. DEPOSITING 





When copper shells are deposited from cold solutions there is a 
slight contraction (about 0.05 per cent). The use of warm solutions 
at a high-current density produces the desired hardness with no 
appreciable change in dimensions. Nickel-copper shells always 
contract by about 0.15 percent. This value is not affected by increas- 
ing the thickness of copper and nickel. Such a change in thickness 
does, however, mateelalie affect the stretch in curving. It is there- 
fore important to keep the conditions for nickel and copper deposition 
constant, and to produce shells of uniform thickness and strength. 


3. CASTING 






Casting always causes a shrinkage of about 0.15 per cent, which 
is not much affected by the temperature of casting or the rate of 
cooling. Any large change in composition of the electrotype metal 
will affect the shrinkage. Shrinkage is decreased by increasing the 
thickness of the shells. 

4. FINISHING 






Hand finishing always stretches the plates, by amounts that vary 
with the extent of the finishing operation. Straightening under 
pressure at ordinary temperature also stretches the plate, but more 
uniformly than does hand finishing. Straightening in a hot press 
causes a contraction of about 0.10 per cent. 


5. CURVING 






The stretch of plain plates in curving is usually less than that 
calculated from the formula for a homogeneous plate; and a type 
form stretches more than calculated. Hot curving may produce a 
slightly greater stretch than cold curving. With a given method, 
most uniform results are obtained when the shells are of the same 
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© thickness and strength and have the same type of pattern, and when 
> the electrotype metal is of the same composition. The temperature 
F of casting and rate of cooling of the electrotype metal have little, if 
| any, effect on the stretch in curving. 

In curving by the ‘“‘nonstretch”’ method, it is possible by adjusting 
the thickness of metal cast upon the face of the plate to reduce 
the stretch to approximately zero. If too much metal is applied, 
shrinkage may be produced in curving. 
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THE HEAT OF IONIZATION OF WATER 
By Frederick D. Rossini 


ABSTRACT 


The heat of ionization of water has been calculated by combining the data on 
heats of neutralization given by Richards and Rowe and the data on heats of 
dilution recently compiled by the present author. The value so obtained is cor- 
roborated by using the data on heats of neutralization reported by Gillespie, 
Lambert, and Gibson. For the reaction, HOO=H++OH-, at infinite dilution 
in water and for the temperature range 10 to 35° C., the heat absorbed is given 
by the equation: AH=13,721—57.9 (¢—18)+0.15(¢—18)? g-calis per mole. 
The uncertainty is estimated to be +16 calories per mole. Using the factor 4.181 
for coverting g-cal,g to absolute joules gives AH= 57,370 absolute joules per mole 
at 18° C. 
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I. INTRODUCTION 


The ionization of water may be expressed by the equation 
H,O=H*t+OH- (1) 


Although the heat absorbed in this reaction is not susceptible of direct 
measurement, its value can be determined in several ways. 

Data on the change with temperature of the equilibrium constant 
for reaction (1) permits calculation of the change in heat content by 
means of the thermodynamic formula:' 


dinK AH 
~=—- 
dr) &# 
This method was employed by Noyes, Kato and Sosman,? who 
measured at various temperatures the hydrolysis of ammonium 
acetate and calculated the ionization constant for water. 
A summation of the heat changes associated with appropriate 


other reactions may lead to the desired heat value. The following 
reactions compose such a series: 


HCl-100H,O + NaOH-100H,0 = NaCl-201H,O (3) 
HCl: © H,O = HCI1-100H,0 + ( —100)H,O (4) 





The nomenclature of Lewis and Randall (“Thermodynamics,’’ McGraw-Hill Co., New York, 1923) 
will be used throughout this discussion. 
? Noyes, Kato, and Sosman, J. Am. Chem. Soc., 32, p. 159; 1910, 
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NaOH: ~ H,O= NaQH-100H,0 + (o — 100)H,O 
NaCl-201H,O + (o —201)H,O0 = NaCl- » H,O (6 
HCl. © H,0 + NaOH: ~ H,0+ © H,0 = NaCl-oH,0+ (2 +1)H,0 


(H+ + Cl- + Nat +OH-)-3 © H,O= (Nat+Cl-+H,0):30H,0 (g 
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or 








Ht +OH-=H,0O (9) 


The reaction expressed by equation (7) is obtained by summing 
equations (3), (4), (5), and (6). Neglecting the somewhat. philo- 
sophical question as to the heat of mixing NaCl-H,O with 2 @H,0, 
etc., which is negligibly small, equation (8) is equivalent to (7) and 
leads directly to the desired reaction (9). In other words, the hea‘ 
of neutralization at infinite dilution may be taken as equivalent to the 
heat evolved when H* combines with OH~ to form H,O, which reaction 
is the reverse of the ionization of water. 

Utilizing the then best available data for the heat changes associated 
with reactions (3), (4), (5), and (6), Richards and Hall * calculated 
a value for the heat of ionization of water in the above manner, 
Richards and Rowe‘ obtained values for the heat of ionization of water 
by plotting heats of neutralization against the concentration and 


pa A 


extrapolating the curves to infinite dilution. Richards and Mair’ 
employed a like procedure. 

In their calculation of the heat of ionization of water, Richards 
and Hall were limited to the use of one series of reactions because of 
the lack of suitable data on heats of dilution for electrolytes other 
than HCl, NaOH, and NaCl. Since that time, however, addi- 
tional data on the heats of dilution of electrolytes in aqueous solution 
have been published, which permits the use of six different series of 
reactions. 

It is the purpose of this paper to combine the data on heats of 
neutralization obtained by Richards and Rowe,’ and incidentally the 
data of Gillespie, Lambert, and Gibson,’ with the values for heats 
of dilution recently calculated by the present author,’ and obtain 
a value for the heat of ionization of water by utilizing six different 
series of reactions similar to the one employed by Richards and Hall. 


II. THE AVAILABLE DATA 
1. HEATS OF NEUTRALIZATION 







































































































A review of the data on the heats of neutralization of uni-univalent 
strong electrolytes indicates that the only reliable measurements- 
especially with regard to purity of materials, proper statement of 
concentrations, and recognition of the proper temperature of the 
reaction—are those of Richards and Rowe,’ and some few by 
Gillespie, Lambert, and Gibson.’ 

















’ Richards and Hall, J. Am. Chem. Soc., 51, p. 736; 1929. 

4 Richards and Rowe, J. Am. Chem. Soc., 44, p. 703; 1922. 

6 Richards and Mair, J. Am. Chem. Soc., 51, p. 739; 1929. 

® Richards and Rowe, J. Am. Chem. Soc., 44, p. 684; 1922. 

7 Gillespie, Lambert, and Gibson, J. Am. Chem. Soc., 52, p. 3806; 1930. 

§ Rossini, B. 8. Jour. Kesearch, 6, p. 791; 1931. 

§ See footnote 6. 
10 See footnote 7, 
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' Richards and Rowe measured the heats of neutralization of HCl, 
' HBr, HI, and HNO; with LiOH, NaOH, and KOH, respectively. 
Lack of suitable data on heats of dilution precludes the use of their 
‘data on HBr and HI in the present calculations. For each salt 
formed, Richards and Rowe performed 6 to 11 experiments, some 
/ having an initial temperature of about 15.5° C. and the remainder 
» 16.5° C. The temperature rise in each experiment was about 4° C. 
| These authors oy the heat capacities of the factors in calculating 
‘the heat capacity of their calorimeter system, and, consequently, 
' obtained the heat of the reaction at the final temperature. The 
heats of neutralization given by Richards and Rowe were reduced 
F to 20° C. from the average values for the experiments having final 
' temperatures of about 19.5° and 20.5° C., respectively. The con- 
' centration of each acid and base was 1 mole in 100 moles of H,0O, 

and the salt solution formed was MX-201H,0. 
' The heats of neutralization calculated by Richards and Rowe 
depend directly upon the accuracy of the values chosen by them for 
' the heat capacities of aqueous solutions of HCl, HNO;, LiQH, NaOH, 
and KOH. Because of the fact that in most cases the heat ca- 
| pacities of aqueous solutions of acids and bases are not known as 

accurately as are those of the corresponding salt solutions, the present 
» author has recalculated the heats of neutralization from the exper- 
imental data of Richards and Rowe, using the heat capacities of the 
formed salt solutions and thus obtaining the heat of neutralization 
at the initial temperature. Inasmuch as the amount of excess 
alkali solutions used by Richards and Rowe was only about 1/1000 
of the resulting salt solution this latter can be assumed to be the 
pure salt solution in computing the heat capacity, without an error 
ereater than 1 or 2 in 100,000. An example of the method of cal- 
culation employed for the data of Richards and Rowe is given in 
Table 1. 

TABLE 1.—Sample calculation 

Data given by Richards and Rowe for experiment No. 13: 

610.22 g HCl-100.06H,O and 615.42 g NaOH-100.19H,0. 

Heat capacity of calorimeter=18.98 g-calis per °C. 

Initial temperature = 15.66° C. 

Final temperature = 19/50" C. 

Average temperature =17.58° C. 

Corrected temperature rise = 3.841° C. 

Calculations by the present author: 

Heat capacity of formed solution (1,225.64 g NaCl-201H;,O) at the 
— temperature, 17.58° C.= (1,225.64) (0.97908)=1, 200.00 g-calis 
ver °C, 

Total heat capacity = 1,218.98 g-calis per °C. 

Total heat evolved = (1,218.98) (3.841) =4,682.1 g-caljs. 

Number of moles of HC]=0.33180. 

Heat evolved per mole of HCl=14,111 g-caljs. 

This gives the heat of reaction at the initial temperature, 15.66° C. 

In the case of NaCl, Richards and Rowe performed five experi- 
ments at an average initial temperature of 16.79° C. and four at 
an average initial temperature of 15.66° C. A value for 16.00° C. 
was interpolated from the two average values, and the small correc- 
tion for differences in concentration was made. In like manner, the 
data for the other five salts were calculated. A summary of the data 
of Richards and Rowe for 52 experiments is given in Table 2. 

he values used by the present author for the specific heats of the 
various salt solutions at 18° C. were computed from values for the 
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apparent molal heat capacities of these solutes in water at 18° () 
These values, given in Table 3, are slightly different from those of 


Richards and Rowe.” 


TaBLE 2.—Data on heats of neutralization 
[From the data of Richards and Rowe] 














Number 

of experi- Pe ,| AH, 

Reaction ments | Pe of re heat 
per- a | absorbed 

formed | Teaction | ss 

| 9-Cals 

°C. | per mol 
HCl- 100 H20+LiOH-100 HzO =LiCl- 201 H20_..---.--.-___-- PERE 6 16.00 | —14,173 
HCl-100 H20+Na0H- 100 H20=NaCl- 201 H20----- aoe s) 16.00 | ~—14 0 
HCl- 100 H20+KOH- 100 H20=KCl- 201 H2:0-_-_--______-- ‘ saad 7 16. 00 | —14 189 
HNOs- 100 H2O+LiOH- 100 HxO=LiN0Os3-201 H:0-_-- a ‘ 10 16.00 | —14, 0% 
HNOs3-100 H20-+NaOH- 100 H2O=NaN0Os3- 201 H20- 4 11 16.00 | —14 03 
HNOs3-100 H20+KOH-100 H2O=KNO3-201 H20__-_-__- ¥ 9 16. 00 —14, 27) 








TasLe 3.—Heat capacities of salt solutions at 18° C. 

















: Heat capacity 

Solution of solution 
g-calis 

g—° C.—! 
Fos tS eres oe» 0. 9842 
NaCl- 201 H20__.---- Aint F | . 9792 
cs { )* =e . 9730 
LiNO- 201 H20.-.--.---- . 9807 
NaNOs3- 201 H20___. | . 9760 
KNOs3- 201 H320.__---- . 9702 





Gillespie, Lambert, and Gibson reported the results of eight experi- 
ments on the heat of neutralization of NaOH and HCl, and one for 
KOH and HCl. Their data, excluding two experiments at 50° C., are 
givenin Table 4. The next to the last column gives the heat capaci- 
ties of the salt solutions as computed from the values for the ap- 
parent molal heat capacities of these aqueous salt solutions." These 
can be compared with the values actually used by Gillespie, Lambert, 
and Gibson. The last column gives the heat of neutralization for 
the given temperature and concentration as corrected by these nev 
values for heat capacity. 


Tas LE 4.—Data on heats of neutralization 
[From the data of Gillespie, Lambert, and Gibson] 
































Calculations by 
| present author 
Tempera- Heat ca- : 
Experi- ture of | 4 verag ply Heat ca- | 
os . eee ok reaction | /+Verage | Sarl Solu-| 4 77 heat | pacity of | 
ae | Solution formed bo tet tempera-| tionat |°y- oe gf | 
No = initial ture erage | 2bSorbed | salt solu-| 4 77 jo 
i tempera- ° Pre tion at | apsorbed 
ture — ‘afl average | SU50rS 
} - tempera- | 
| ture 
| g-caloo | g-calo 
1 A PR, °C, = \g-1° C1) per mole 
eR es NaCl- 1080.7 H30...----.-_-_-- 20. 00 20. 35 0. 9961 
, &e --| NaCl- 566.5 H20____- | 20. 00 20. 68 . 9924 
3... .| NaCl- 1101.4 H2O_.______-- 24. 99 25. 33 - 9957 
4 ...--| NaCl-568.8 Hs0_____-- : 32. 27 32. 92 . 9922 
au. ..-| NaCl- 566.0 H:O_-..___- wu 32. 27 2. 95 - 9922 
= | NaC)- 520.5 H30...._....._-- 32. 30 33. O1 . 9916 
y KC]l- 847.8 Hs30-......-......- 20. 00 20. 71 | . 9896 





1! Rossini, B. S. Jour. Research; 1931. 
8 Richards and Rowe, J. Am. Chem, Soc., 41, p. 770; 1921. 
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2. HEATS OF DILUTION 


The changes in heat content associated with reactions (4), (5), 


» and (6) can be obtained from the values of &,—4,°, the relative ap- 
© parent molal heat content of the solute, given for aqueous solutions 
© of HCl, HNO;, LIOH, NaOH, KOH, LiCl, NaCl, KCI, LiNO;, NaNO;, 


» and KNOs, at 18° C., by the present author in a previous publication.” 


For example, in reaction (4) 
HCl- oH,O=HCI- 100 H,O+ (o — 100) H,O (4) 


AH = 4, — $;° (10) 


» where &, is the apparent molal heat content of HCl in the solution 
' HCl-100 H,O, and ®,° is the apparent molal heat content of HCl 


at infinite dilution in water. 
III. THE CALCULATED RESULTS 


The values cited above for heats of dilution are given for 18° C.; 
the data of Richards and Rowe on heats of neutralization, shown in 
Table 2, are for 16° C.; and those of Gillespie, Lambert, and Gibson 
on heats of neutralization are given at temperatures of 20°, 25°, and 
32.3° C. (Table 4.) 

Because the changes in heat content for these reactions have large 
temperature coefficients, it is necessary that these latter be known 
accurately for the proper conversion of the various data to 18° C. 
Once the values for the heats of neutralization are calculated to 18° 
C., it becomes a relatively simple matter to sum the changes in heat 
content for reactions (3), (4), (5), and (6), and obtain thereby values 
for the change in heat content for reaction (9) 


H+ +OH-=H,0 (9) 


For reaction (3) 
the temperature coefficient of the change in heat content is given by 


AC, = 201 Om + Penacn — (100C%, + Beazer + 100C%, +®enaon)) (11) 


AC, - Cn + Bewacn — Pecacrn — Pecwaom (12) 


where Penacns P71) and ®.naon) are the apparent molal heat 
capacities of NaCl, HCl, and NaOH in their respective solutions and 
(x is the molal heat capacity of pure H,O. In correcting a heat of 





8 See footnote 8, p. 848. 
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neutralization from t° C. to 18° C., values of ®, for the average te. 


on be 
perature, °C., are needed. For converting the data of Richar 


and Rowe “i 16° to 18° C., values of &, for 17° C. were used; an) 
for converting the data of Gillespie, Lambert, and Gibson from 323° 
25°, and 20° to 18° C., values of #, for 25°, 21.5°, and 19° C., respec. 
tively, were used.” 

A comparison of the temperature coefficients given by Richards 
and Rowe” for the six different heats of neutralization “with thos 
calculated by the present author is given in Table 5. 

The data of Richards and Rowe on heats of neutralization at 16°C. 
as given in Table 3, can now be converted to 18° C., and combined 
with the appropriate heats of dilution, to give values for the heat o 
ionization of water. ‘The results of these calculations are summarized 
in Table 6. 

The calculations from the data of Gillespie, Lambert, and Gibson 
are summarized in Table 7. 


TaBLeE 5.—Temperature coefficient of heat of neutralization, at 18° C. 


[For the reaction: HX-100 H:,0+MOH- 100 H2.O=MX- 201 H20] 





| AC;, in cal, mole~!°C,-1 





| 
| Richards Present 
and Rowe author 


Solution formed 





Fle bP (ere 55. 42. ¢ 
NaCl-201 H20..-.....-. §2. 46. ( 
KCl. 201 | STi 8. } 46. 

LiNOs;- 53. 6 40. 8 
NaNOs-2 2 50. ! 45.7 
KNOs3- 201 H30 .| 5. 47. § 














TABLE 6.—The heat of neutralization at infinite dilution 


[Combining data of Richards and Rowe with data on heats of dilution.] 





HX-100 H,.0+MOH- 100 Heats of dilution at | Ht+ 
H20=M X-201 H20 18° C., g-cal.is per mole | OH-= 
H20 





AH at AC, at | AH at 


Solution formed | 
16° C, 17°C is° C . 

Acid Base 

g-cal.is g-cal.1g g-c alts 18 |e. 6,%p,—a,° 

pet mole~!° per villas igdilidals 
mole C-1 mole | 





LiCl.- 201 H2O.-.-..--- ...| —14, 173 42.6 | —14, 088 305 24: —174 | —13, 
NaCl- 201 H20_..------ | —14, 099 3. 3 | —14, 006 305 36 —48 | —13,7 
KCi-201 HsO.......... —14, 189 b —14, 096 305 37 —46 | —13,7 
LiN O3- 2 | —14, 024 .1 | —13, 942 +86 24; —132 | —13, 745 
NaNO3-201 H20O —14, 031 3.0 | —13, 939 1 3s +83 | =13, 731 
KNO3-201 H2O f — i4, 271 47.6 | —14, 176 37 +229 | —13, 724 | 





Mean 








14 See footnote 11, p. 850. 
6 Richards and Rowe, J. Am, Chem. Soc., 44, p. 699; 








Heat of Ionization of Water 


; Rossini] 
TaBLE 7.—The heat of neutralization at infinite dilution 


{Combining data of Gillespie, Lambert, and Gibson with data on heats of dilution] 





HX-nH20+MOH-nH20O=MX- Heats of dilution at 18° | H++OH- 

(2n+1) H2O C., g-cal.1s per mole =H20 
| Devia- 

| tion 

| AH at from 

Base | Salt 18° C, mean 
>n—Paol|PpO—P,|  g-cal.ss 
| | per mole 


| 


|— : 


luti AH at | AC» at 
Solution formed there Kick ll tls), 
temper- temper- 2 id 
lature. | ture. | C.g-cal.is er seated 
= g-calis | mole! | mole | 
per mole oC =-1 





vaCl- 1080. 7H2O . - . 20. 3, 77 52.6 — 13, 879 +140 
1C1- 566. 520 -_.-- 20. 3,7 50.5 |—13, 895 | +190 
NaC}-1101. 4H20---...| 24. 9 3, 56 .8 |—13, 931 | +139 
YaCl- 568. $H20...-- at eas 3, 192 8.8 }—13,890 | +190 
1—13,951 | +190 
3} 1-13, 958 | +197 
—13,959 | +193 


+84 | 
Ler | 
+-83 
+85 
+85 
+84 | 
+123 —6 —13, 703 | +15 





| 
| 
1 | 
\d 
| 
| 
$$$ —— -+|— ——|- —-|—— 
| 
| 
| 
| 
| 
ol 





Mean neccchalon - — 2 | wall ~13, 718 | +27 





THE TEMPERATURE COEFFICIENT 
The temperature coefficient of the heat of reaction (9) 
H*++OH-=H,0O (9) 


can be calculated from the values of the heat capacities of the sub- 
stances involved. That is 


AO, = ©, %. ee ee (13) 


(HO) 


The heat capacity of H,O is known, and the sum of the partial molal 
heat capacities '® for H* and OH™ at infinite dilution can be deter- 
mined as follows: 

, the partial molal heat capacity of the solute at infinite dilu- 
tion, is known for NaCl, HCl, and NaOH. Now 


ce } 4 YO | CO ° nl re, ° 
P (NaOB) ie DP (Nat) ' »P (OH-) 


(14) 


(HC) 


TF dan hig iin tte bit’ (15) 


P (NaC)) D (Nat) 


Subtracting equation (15) from equation (14) gives 


ry -_( ° C fe) Cc ° 
’ P (NaOH) P (NaC}) p (Ht) 


417 0 N 
(CHC) P C, (OH-) (16) 


From the values for ®°, (which is to C,.°) compiled by the 


present author,” one finds, for 18° ¢ 


For an extended discussion of apparent and partial molal heat capacities in aqueous solution see (a) 
ndall and Rossini, J. Am. Chem. Soc., 51, p. 323; 1929; and (6) Rossini, B. 8. Jour. Research, 4, p. 313; 





’ See footnote 11, p. 850. 
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Ge... +0,° —39.9 g-cal. mole! ° C.- (17) 


Oon- 








And for reaction (1) 
H,O = Ht +OH- 


AC, 





= —39.9—18.0= — 57.9 g-cal. mole ° C.“! 


The temperature coefficient!” of C,° a1) + Cp°oq-) Can be taken as 0). 












g-cal. mole! °C.-*. For any temperature ¢, near 18° C., the ten. 
perature coefficient of the heat change associated with reaction (1) is 


AC, = —57.9+0.3(¢—18) g-cal. mole °C.~! (19 





V. DISCUSSION 


The value for the change in heat content for reaction (9), as obtained 
from the data of Richards and Rowe on heats of neutralization, is 











AH = — 13,721 g-cal.,, per mole at 18°C, (20) 


or for the ionization reaction 





AH = 13,721 g-cal.jg per mole at 18°C. (21 





















The average deviation of the six individual results from this value is 
+13 calories per mole. 

Using the data of Gillespie, Lambert, and Gibson on heats of nev- 
tralization, one obtains for the ionization of water 


AH = 13,718 g-cal.is; per mole at 18°C, (22 








Here the average deviation is +27 calories per mole. Inasmuch 
as the experiments of Gillespie, Lambert, and Gibson consist of 3 
at 20° C.,1 at 25° C., and 3 at 32.3° C., as compared with the 52 
experiments of Richards and Rowe at 16°C., the value given by equa- 
tion (21) will be selected as the “‘best’’ value for the change in heat 
content which accompanies the ionization of water at infinite dilution. 
The value obtained from the Gillespie, Lambert, and Gibson data 
serves as an excellent verification of the selected value. The agree- 
ment also serves to establish the reliability of the heat content values " 
as well as the apparent molal heat capacity values’® which were used to 
convert the data of Gillespie, Lambert, and Gibson to infinite dilu- 
tion and to 18° C., respectively. 

The error in the value calculated for the heat of ionization of water 
by the present method depends upon the accuracy of the determins- 
tion of the various heats of neutralization and of the heats of dilution. 
The absolute accuracy of the final values for heats of neutralization 
calculated from the data of Richards and Rowe, using the new heat 





\§ See footnote 8, p. 848. 
1% See footnote 11, p. $50. 
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FigurE 1.—Plot of the data on the 
heat of ionization of water 


The ordinate scale, in g-cal.1s per mole, gives 
the heat absorbed in the reaction, Hx0=H+t-+- 
OH-, at 18° C, and infinite dilution. The 
circle shows the estimated error in the selected 
“best” value, 13,721. The values obtained 
from the data of Richards and Rowe on heats 
ofneutralization are indicated by @; and those 
from the data of Gillespie, Lambert, and Gib- 
son by ©. The latter values are from single 
experiments on the heat of neutralization, 
while the values of Richards and Rowe for 
heats of neutralization are each the average 
of 6 to 11 experimental values. 
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capacity values, can be placed at +8 
calories per mole. If the uncer- 
tainty in the heats of dilution for 
reactions (4), (5), and (6) is estimated 
to bein each case + 8 calories, the error 
in the final selected value for AH may 
betakenas + -/8?+ 8?+8?+8?= +16 
calories per mole. 

A graphical representation of the 
results obtained for the heat of 
ionization of water from the data 
of Richards and Rowe, and of Gilles- 
pie, Lambert, and Gibson, on heats 
of neutralization is shown in Fig- 
ure 1, 

It is interesting to note the values 
obtained for the heat of ionization 
of water by the investigators cited 
at the beginning of this paper. 

For the ionization reaction 

H,0 = H*+OH- (1) 
the values, converted by means of 
the proper temperature coefficient 
to 18° C., from whatever tempera- 
ture they were given are: 


AH, at 18° C., in g-cal.is per mole 


Noyes, Kato,and Sosman_ 14,070. 
Between 13,730 
and 13,790. 
13,760. 
13,760. 


Richards and Hall 
Richards and Mair-_-_-_-__- 
This paper 


VI. CONCLUSION 


By combining the data on heats 
of neutralization obtained by Rich- 
ards and Rowe with appropriate 
heats of dilution, there has been 
found for the reaction 


H,O= Ht +OH- 


AH = 13,721 +16 
g-cal.,; per mole at 18°C, 3) 
and 
ACp= — 57.9+0.3(¢—18) 
g-cal.,, per mole per degree (24) 
For the temperature range of about 
10° to 35° C., the heat absorbed 
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in reaction (1), at infinite dilution in water, can be expressed by the 
equation 


AH = 13,721 —57.9 (€—18)+0.15 (¢—18)? g-cal.ig per mole (25 
Here ¢ is in °C. and the unit calorie is the actual heat capacity of | 
g of water at 18° C. Using the factor 4.181 for converting g-cal,, 
to absolute joules, gives, per mole at 18° C., 


AH = 57,370 absolute joules (26) 


WASHINGTON, February 14, 1931. 
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FACTORS AFFECTING THE STRENGTH OF MASONRY 
OF HOLLOW UNITS 


By Douglas E. Parsons 


ABSTRACT 


A study was made of the data reported in University of Illinois Bulletin No. 
97 and Bureau of Standards Technologic Papers Nos. 238 and 311 in order to 
determine the importance of some of the factors affecting the strength of masonry 
f hollow units. For end-construction walls a simple empirical equation was 
found to express fairly well the relation between the strength of the walls and cer- 
tain easily determined properties. For both end and side construction, the 
geometrical properties of the units appeared to have an important effect, and 
the strength of the units was not a close measure of the effect of their charac- 
teristics on the strength of the walls. 


CONTENTS 


I. Introduction 
II. The test data 
III. Studies of the data_- 
1. General _ _ _- 
End-construction masonry - 
3 Side-construction masonry ‘ 
IV. Concluding remarks..c 2... gos5- eee 5 --u cx ee 


I. INTRODUCTION 


Allowable working stresses for walls of hollow units for buildings 
usually are the same for units having equal compressive strengths. 
These stresses are different with different mortars, but appear to be 
based upon the assumption that the compressive strength is the only 
property of the unit affecting the strength of walls. The strength 
values obtained in tests of large masonry specimens indicate strongly 
that there are other properties of the units which have an important 
effect on the strength of masonry, and it seemed worth while to study 
further the test data and to find, if possible, a better measure of 
masonry strength. 


II. THE TEST DATA 


The reports of Talbot and Abrams', Whittemore and Hathcock? 
and Stang, Parsons, and Foster’ give the results of tests on large 
masonry specimens ‘of hollow clay units in which the arrangement 
(bonding) of the units was similar to that used in ordinary building 
construction. Moreover, these reports contain rather complete data 
on physical tests of the materials and descriptions of the method of 
construction and the workmanship. A summary of the data per- 
taining to compressive tests under concentric loading of end-con- 
struction masonry (cells vertical) with broken joints is given in Table 
1. A similar summary of the data pertaining to side-construction 
walls (cells horizontal) is given in Table 2. 





' Tests of Brick Columns and Terra Cotta Block Columns, Univ. of Illinois Bul. No. 27; 1908. 
‘Some Compressive Tests of Hollow-Tile Walls, B. S. Tech. Paper No. 238, 1923. 
‘Vompressive and Transverse Strength of Hollow Tile Walls, B. 8. Tech. Paper No. 311; 1926. 
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No other data‘ on tests of large masonry specimens of either holloy 


clay or hollow concrete units with as complete information on mate. 
rials and workmanship have come to my attention. Indeed, there 
appears to be no published information of this character on the 
strength of masonry of hollow concrete units. 

Values of the strength of the mortar specimens are lacking in q 
few cases as indicated in Table 1. And the strength (2,150 Ibs. /in.’) 
given for the 1C:1/4L:3S mortar at the age for testing the masonry 
was estimated from the values of strengths at the ages of 7 days and 
6 months. 

The ratios 6 in Table 1 of bearing area to gross area of masonry 
are estimates based upon the dimensions of the units and descriptions 
(and observations in the case of the}tests reported in Bureau of Stand- 
ards Technologic Paper No. 311) of the workmanship and method of 
laying the units. 


III. STUDIES OF THE DATA 
1. GENERAL 


It is expected that the compressive strength M of a masonry speci- 
men would depend chiefly on the compressive strength of the units 

, the compressive strength of the mortar m, other strength proper- 
ties and the elastic properties of the mater rials of the units and of th 
mortar and on the various dimensions describing the shapes and sizes 
of the masonry and its component parts. If these quantities are the 
only ones having a significant influence, the relations between them 
may be expressed according to the principle of dimensional homo- 
geneity® by the following equation: 


oi eel Oe 
M=KuF (=. eo tn) ( 


K=an unknown numerical constant. 

F=an unknown function of the dimensionless ratios within 
the brackets. 

R,, R..-R»=ratios of other strength properties and of the modul 
of elasticity of the materials to u. 

T1, T2--’, =ratios of other linear dimensions of the masonry to 
a chosen linear dimension L, or ratios of areas to a 
chosen area. In any case, geometrical ratios describing 
the shape and size of the masonry and its component 
parts. 

It is, of course, optional which of the strength properties or elastic 
moduli of the materials is chosen in the place of u, provided the /’s 
contain all of the quantities of the same dimensions that have an 
influence on M. The form of the function F and the value of the 
constant AK can be determined only from test data. 

The obvious direct way to determine the form of the function F 
by a set of new experiments would be to vary the dimensionless ratios 





‘ Data on the strength of walls 3 feet long and 3 feet high ¢ are rey orted by Ingberg and Foster in the paper 
“Fire Resistance of Hollow Load-Bearing Wall Tile,’ B. 8S. Jour. Research, 2 (RP37); 1928. 

5 Also known as the principle of similitude or when a; wniied to mechanical problems as the law or principle 
of dynamic similarity. Those not already familiar with the subject will enjoy reading Dr. Edgar Bucking- 
ham’'s excellent presentation in the paper ‘“‘ Model Experiments and the Forms of Empirical Equations, 
Trans, Am, Soc, Mech, Engrs,, 37, p, 263; 1915, 
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of F, one at a time, holding all the others constant. Then writing 
the equation (1) in the form 


—KF ( i? aa ra) (1a) 


M : ; 
the relation between - and each of the ratios could be determined. 


When using the common materials for hollow walls, however, it is 
sometimes impracticable to vary one ratio over a wide range without 
affec ting the values of some of the others. For example, mortars 
having a wide range in compressive strengths also will vary in other 
strength properties and in elastic properties. The present purpose, 
however, is to find what useful information may be derived from 
existing data. 

[t is apparent from the outset that, with so few data and so many 
variable factors, the form of the presumably complete equation (1) 
or (la) can not be determined at present. In order to obtain useful 
results it becomes necessary to simplify the expression by omitting 
some of the ratios. In doing this it is necessary to omit none that 
has a large influence if the resulting equation is to be a fair approxi- 
mation, and, if the results are to be useful, to include only those 
terms based on values that may be readily determined in practice. 


2. END-CONSTRUCTION MASONRY 


The test results do not provide sufficient data of the right sort to 
determine the manner in which the ratios (rn, 12--Tn) “describing 
the sizes and shapes of the masonry, masonry units, and mortar 
joints enter into the equation. If as a first approximation one ratio 
were to be selected, which by itself would serve to indicate roughly 
the effects of all of the geometrical ratios, it seems likely that the 
one best suited would be the ratio b of the bearing area to the gross 
area of the masonry. Numerically the value of b is equal to the 
ratio of the total area, at the plane of a joint, of the material of the 
units given a bearing by the mortar, to the gross area of the masonry 
at the same plane. If as a further approximation all the R’s are 
omitted, equation (1) may be written 


M=kuf(, b) 
aaW(2, b) 


where f is some unknown function of the ratios within the brackets 
and £ is an unknown numerical constant. 

By comparing the distributions of the values in Figures 1 and 2,° 
t appears that the value of 6 has an important effect on _ the strength 


(2) 





e proportions for and the properties of the 1C:3S and the 1C:1/4L: 3S - mortars were so near alike that 
symbol was used in the figures to represent values relating to both. 


Pat 
i 





Bureau of Standards Journal of Research 





Ce a | 


Her Tro SYMBOLS 


/4L:38 Morrar 

1C: 58 ” 

10:38 ano IC:4L:33S Morrars 
1C:lEL: 45 Mortar 

IC: |EL:68 ; 

10:25 





| 





8 
8 
8 





X@peebpo 


| 





: 








S 
S 








S 
S 
S 





® ® 




















Compressive strength of rnasorry, M, lb. per sq. itt 























oO 





oO 1000 2000 3000 4000 5000 6000 
Compressive strength of units, (gross area) u, 1a per sq. in 


Fiagure 1.—Relation between compressive strength (gross area) of end-con- 
struction masonry and compressive strength of units 


The proportions for the mortars are by volumes, with the abbreviations C for cement, L for lime, 
and 8 for sand. 
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Figure 2.—Relation between compressive strength (gross area) of end-con- 
struction masonry and the ratio of net bearing area in masonry to gross area 


The proportions for the mortars are by volumes, with the abbreviations C for cement, L for lime, 
and § for sand, 
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fPorsons] 
Bot the masonry M. Since the data are inadequate for determining 
ihe relation between af and b while . remains constant, all of the 
values are plotted in Figure 3. By considering, one at a time, the 


' ie 0 , 
Fyalues for which z varied only over a small range there seemed 


F M . 
‘be a rough relation between ir and 6. Such meager data do not 


Fiustify any refinements in the selection of a cumbersome form for 
representing the relationship. <A linear relation probably is close 
enough for the purpose and on that basis the approximation becomes 


by ki (5 (3) 


. . m™m ye ° 
‘where f is an unknown function of = and A; is an undetermined 


' numerical constant. 
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5 ; M 
FiaurE 3.—Relation between wn and b, where 


M=compressive strength of masonry (gross area) lbs./in.? 
u=compressive strength of units (gross area) lbs./in. 
b=ratio of bearing area in masonry to gross area 


The proportions for the mortars are by volume, with the abbreviations C for cement, L for lime, 
and 8 for sand, 


Without burdening the reader further with the details of developing 
the approximate equation, it was found that the following relation 
fits the data for end construction masonry as well as any other 
equally simple form investigated : 7 

M_y,(m 
or =K 4 
bu Nau, (4) 


M= K,b vn imu 


’ Othe r > forms of equations were investigs ited, t but the n main point to be emphasized is the necessity for 
sidering other factors than wu and m, and this is done equally well by the more simple form (4). 
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The value of K, is approximately unity. A comparison betwee 
values of masonry strength calculated by means of equation (4) with 
test results is shown in ‘Fi igure 4. Different empirical equations {o; 
expressing the relations between the quantities of equation (2) may 
be found readily, but with so many factors (quantities) left out of 
consideration it is useless to begin splitting hairs here. 

The scattering of the plotted ‘points of Figure 4 gives an indic: ation 
of the incompleteness of equation (4). <A part of the scattering ma 
be attributed rather definitely to known causes. For example, the two 
values for which the ratio of test to estimated strength are lowes 
pertain to specimens 12 feet high and only 6 inches thick. It seems 
likely that the large ratio of heizht to thickness of the walls accounts 
in part for the unusually low strengths of these specimens. Of the 
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Ficure 4.—Relation between compressive strength (gross area) of end-con- 
struction masonry and by mu, where 
b=ratio of bearing area in masonry to gross area 
m=compressive strength of mortar specimens, lbs./in.? 
u=compressive strength (gross area) of units, lbs./in.? 


three values from the data of Talbot and Abrams, one pertains to 
masonry specimens built with units specially selected because of 
their uniformity of size and shape, another to units having warped 
bearing surfaces and the third to a ‘poorly laid”? masonry specimen. 
Obviously, such factors as these are not taken into account by 
equation (4). 

Among the several other possible causes of the scattering, one 
difficult to evaluate is that the properties of the mortar in the joints 


probably are different than those of the mortar specimens. Water 


absorbed from the mortar by the units and differences in curing con- 
ditions would be likely to affect the relations between the properties 
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‘of the joints and specimen mortars. Another possible cause for devia- 
Ftions is the fact that the size and shape of the mortar specimens were 
not the same in the three investigations. Such variables also suggest 
Freasons Why it would be quite useless to attempt a great refinement in 
‘selecting a form of equation to represent the available data. 


3. SIDE-CONSTRUCTION MASONRY 


| The relation between the comprehensive strengths (gross area) of 
‘ <ide-construction masonry and of the units is shown in Figure 5. 
‘Here the relation seems to be somewhat more regular than that for 


‘end construction masonry illustrated in Figure 1. 


A part of the lack of regularity in the distribution of the plotted 


points in Figure 5 is assignable to known causes. The masonry 
' specimens designated Nos. 33 and 35 in Table 2 were built with units 





, s+, . Fs 


— Hey To SYMBOLS 
ce) IAL: 3S Mortar 
@ /C: 35 ano 10:41:38 Monrars 
4 /C:/4L:48 Morrar 
® /C:/ZL:6S5 ” 








8 





§ 











streng’?) of rrasolury, 11, 


1b, per sq. 177, 








: 
: 
: 
8 



































0 

QO 200 400 600 800 1000 1200 1400 1600 1800 
Compressive strength of urits (gross area) u, 
1b. per SQ. 177. 


Figure 5.—Relation between comprehensive strength (gross area) of 
side-construction masonry and compressive strength of units 


The proportions for the mortars are by volume with the abbreviations C for ce- 
ment, L for lime, and S for sand. 


having irregular bearing surfaces. A portion of the bearing shells 
near the vertical mid plane of the walls was concave toward the 
mortar joint. It is understood that these tiles were formed in this 
manner purposely in order to lessen the likelihood of there being a 
continuous mortar joint through walls. The recesses in the bearing 
shells resulted in there being a continuous bearing only for the outer 
shells; for ordinary tile of this design about 70 per cent of the net 
sectional area would be given a direct bearing. 

Differences in the ratios between the thickness of the bearing 
(horizontal) shells and the horizontal distances between the supports 
provided by the vertical shells and webs probably is another major 
cause for the lack of a closer relationship between the variables shown 
in Figure 5. Comparing No. 36 with No. 37 in Table 2, it is seen that 
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the masonry with 3-cell tiles was stronger than that with the 2-cj 
units. In fact, all of those having relatively closely spaced vertics| 
shells and webs (Nos. 30, 37, 38, 39, and 40) lie above an average lin. 
through the points of Figure 5. It ‘has been observed ® that in many 

cases the cracking of the horizontal shell in contact with the mortar 
was the first sign of the failure of side-construction walls in c ompressiye 
tests. It seems apparent, therefore, that the bending strength of th; 
shells is important with some unit designs. 

However, the available data relating to side construction masonry 
are inadequate for a determination of the form of equation (1). I 
fact the two important factors mentioned in the preceding paragraphs 
can not be well evaluated at present. 

For the end-construction walls account was taken of the workman. 
ship by the ratio b, this ratio being a function of the design of the units 
and the proportion of the shells and webs given mortar bearing by the 
mason. No such simple way of evaluating the effect of workmanship 
in side-construction walls is apparent. It is probable, however, tha: 
there were marked differences in the workmanship of the different 
masons. Walls Nos. 32 and 42 of Table 2 were of identical construc. 
tion except for the mortar and workmanship. ‘The strength (M = 4); 
lbs./in.?) of No. 32 with 1C:1—-1/4L:4S mortar was apprec iably greater 
than the strength (IZ =360 Ibs./in.2 ) of No. 42 with the much stronger 
1C:3S mortar. The tests of Ingberg and Foster® gave, on the aver- 
age, higher strengths for walls with cement mortar than for those wit! 
1C:1L:45 mortar. Hence, it is logical to assume that with the sam 
workmanship No. 42 would have been stronger instead of weaker tha 
No. 32. 

Without having some way for evaluating the effects of workman- 
ship, direct comparisons to determine the relations between wall 
strength and properties of the mortars are not easily made. About 
the only general statement warranted is that the use of stronger 
mortars usually results in stronger masonry. 

In so far as the units are concerned, the data indicate that the three 
following te cpimaan have an important effect on the strength of 
masonry: (1) Compressive strength, (2) ratio of thickness of bearin 
shells to maximum span between vertical supports (shells and we bs), 
and (3) features of design affecting the proportion of area given : 
bearing at bed joints. 
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IV. CONCLUDING REMARKS 





Equation (4) must be considered as only a rough approximation. 
It is obvious that no such simple expression could be correct since it 
does not take account of many factors, such as the elastic properties 
of the materials and the size and shape of the masonry and its parts 
In addition to the cause mentioned previously, the methods of testing 
masonry, the units and the mortar would affect the values de “igi 
from the tests and, hence, the agreement between any predicted and 
determined values. 

In addition to the general size and shape of the masonry and thi 
units, two other geometrical properties, namely, the thickness of the 


‘B. S. Tech. Paper No. 311, p. 339. 








*See Table 2, B. S. Research Paper No. 37. 
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Fortar joints and the regularity of the bearing surfaces of the units 
sare of importance. A discussion of the effect of the latter is given by 
iTalbot and Abrams, but the effect of neither of them has been con- 


sidered here. 


Nevertheless the test data indicate that properties of the units 


Fother than the compressive strength must be considered for even 
Syough estimates of the compressive ‘strength of masonry. The values 


shown in Figures 2, 3, and 4 indicate that the ratio b of the net bearing 


‘area to the gross see tional area of the masonry is an important prop- 
| mp with end-construction masonry. For side-construction masonry 


t appears that the two factors (1) the ratio of thickness of bes ing 


i shells to maximum span between supports and (2) the proportion of 
Farea given a bearing, should be considered. 


WasHincTon, February 24, 1931. 








RP 311 


DETERMINATION OF THE BENZENE AND THE NORMAL 
HEXANE CONTENT IN A MID-CONTINENT PETRO- 
LEUM * 


By Johannes H. Bruun? and Mildred M. Hicks-Bruun 


ABSTRACT 


The benzene content of fractions of a mid-continent petroleum has been deter- 
mined by converting the benzene to m-dinitrobenzene. All of the benzene was 
found as @ constant boiling mixture with the hexanes in the fractions boiling 
below the boiling point of pure benzene. The large fraction boiling between 69° 
and 70° C. contained the greatest amount of benzene. The greatest concentra- 
tion (6.5 per cent of benzene) was found in the 68° to 69° fraction. Based upon 
the crude petroleum the benzene content was found to be 0.08 per cent. A 
sample of pure benzene was isolated from its constant boiling mixture by extract- 
ing with aniline and successive equilibrium melting. 

The n-hexane content of the crude petroleum has been estimated to be about 
ne-third per cent. This is more than the combined amounts of the other isomers 
of hexane. A sample of pure n-hexane has been isolated by treatment with 
chlorosulphonie acid. 

Time-temperature cooling curves, as well as the following physical constants, 
have been determined for the isolated hydrocarbons: Freezing point, boiling 
point, density, and refractive index. All of these constants agree well with those 
of the synthetic hydrocarbons. 

The benzene-hexane ratio of a petroleum oil is the determining factor with 
respect to the temperature at which the benzene distills. 

The boiling ranges have been accurately determined by means of a continuous 
boiling-point apparatus which has been placed in the distillation system between 
the top of the column and the condenser. The advantages of this method are 
discussed. 
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I. INTRODUCTION 




























The presence of benzene in American petroleum was established 
as early as 1863 by Schorlemmer’*, who found that some nitrobenzey 
was formed, when a fraction from Pennsylvania petroleum a 
boiled with nitric acid. Mabery,* in his pioneer work on Amer. 
can petroleum, found that, in the distillation of an Ohio crude oj! 
benzene concentrated in the fractions boiling between 75° and 86° ( 
The greatest concentration of benzene (15 per cent) was found jy 
the fraction boiling between 79° and 81° C., or, in other words, around 
the boiling point of pure benzene (80.2°C.). Based upon the crud 
oil, Mabery estimated that the Ohio petroleum contained 0.017 p: 
cent of benzene. Mabery’s method consisted of nitrating the ber. 
zene to mononitrobenzene and weighing it as such. Mabery and 
Hudson® found that the bulk of the benzene from a California pe. 
troleum also collected in the fractions distilling between 79° and 81°( 

On the other hand, Young,® in his work on Pennsylvania petroleum 
found that all of the benzene distilled over far below its norm: 
boiling point and that it concentrated largely in the fractions boiling 
at about 65° to 66° C. The reason for this is that constant-boiling 
mixtures (with minimum boiling point) are formed between benzen: 
and the hexanes. Young determined the benzene by nitrating ii 
to its dinitro derivative. 

The isolation of the isomers of hexane has been discussed in 
previous paper.’ In the present investigation some additional in- 
formation has been obtained concerning the approximate amounts 
in which these isomers of hexane are present in the petroleum. 


II. EXPERIMENTAL PROCEDURE 
1. DISTILLATION * 


The material for this investigation was obtained from an Okl- 
homa oil field. 600 gallons of crude oil from a single well ® was sub- 
jected to a preliminary distillation in an oil refinery. The fractions 
thus obtained were first redistilled twice in a rectifying still! (20- 
liter capacity) with 20 bubbling-cap plates of steel." The lower boi 
ing fractions from these distillations were then subjected to a fina 
distillation through a 30-plate rectifying column" of Pyrex labors 
tory glass. This distillation was carried out at a rate of about? 
ml per minute with a reflux ratio of about 10:1. The final dis 
tribution of the fractions with respect to their boiling range is shown 
by Graph I, in Figure 2. 

Instead of the usual method of reading the condensation tempers- 
ture at the top of the column the following method has been de- 
veloped, in accordance with a suggestion by E. W. Washburn. The 





4 


C. Schorlemmer, Ann., 127, p. 311; 1863. 
C. F. Mabery, Proc. Am. Ac., 31, p. 34; 1894. 
5 C, F. Mabery and E. J. Hudson, Proc. Am. Ac., 36,, p. 259; 1901. 
6S. Young, J. Chem. Soc., 73, p. 914; 1898. 
7 J. H. Bruun and M. M. Hicks-Bruun, B. 8. Jour. Research, 5, pp. 933-42; 1930. 
* The final distillation fractionation was carried out by 8. T. Schicktanz and a staff of still operators 
* Well No. 6 of the South Ponca Field, Kay County, Okla. For properties, see E. W. Washburn, J. # 
Bruun, and M. M. Hicks, B. S. Jour. Research, 2, p. 469, Table 1; 1929. 
10 E, W. Washburn, J. H. Bruun, M. M. Hicks, B. 8. Jour. Research, 2, pp. 470-3; 1929. 
uJ, H. Bruun, Ind, Eng. Chem., Anal, ed., 1, p, 212; 1929, 
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distillate, after leaving the top of the column, is allowed to pass 
down into a Cottrell” “boiler as shown in Fi igure 1. The hot liquid 
lished from the top of the retifying column enters the boiler at A, and by 
nzey [ie means of an electric heating coil, H, in the bell, B, is kept boiling 


wai continuously, A mixture of vapor and boiling liquid s ascends from the 


mer. 


| 

le oll - | | + 
a | 

> ( Ae a 

nd ir = 

‘ound 


crude | 
ha he igs FROM DISTILLING 
ie | COLUMN 


and 
a] pe- 
iat 
eum 
rm 
iling 
iling 


Zen: 
























































condehce| 


Figure 1.—A continuous boiling point apparatus 
The lower part of the apparatus in Figure 1 is insulated by a vacuum jacket, as shown. 








bell, B, through the narrow tubes, N, and is discharged against the 
thermometer, T. The liquid over flows continuously ‘through O into 
the condenser and receiver unit of the still. The main adv antages 
of this method of reading the distillation temperatures are: (1) 
The method gives the true boiling point rather than the ‘condensa- 





2 F, G, Cottrell, J. Am. Chem. Soc., 41, p. 721; 1919 
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tion” point, which is the usual one obtained at the top of a Tectifying 
column; and (2) greater precision as well as convenience in reading 
the thermometer, because the Cottrell boiler may be placed at any 
convenient level in the laboratory. 
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Figure 2.— Mass data plotted against boiling range of fraction 
2. NITRATION OF BENZENE 


As a result of the distillation, 37 one-degree cuts were obtained 
between 44° and 81° C., as shown in Graph I, Figure 2. The benzene 
in each of these cuts was then converted into m-dinitrobenzene by 
a nitrating method similar to that previously used * for the detern min- 


8 J, H, Bruun, R. T, Satie @. 7, Sdulektons, B.S Seen esteech, pp. 363-367; 1931, 
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ation of toluene. However, as it is well known that the reactivity 
of benzene toward nitric acid is less than that of its homologues, it 
was necessary to use a large excess of nitrating acid and very vigor- 
ous agitation for a long time in order to obtain quantitative conversion. 
During the nitration the temperature was kept below 10° C. by means 
of a continuous stream of tap water through the cooling bath. 

Before nitration, the absence of olefine hy drocarbons was proved 
by testing with bromine. A check run on an unnitrated benzene- 
free fraction indicated that the nitrating treatment did not remove 
or attack the aliphatic hydrocarbons. A nitration was also carried 
- on a synthetic mixture of olefin and aromatic-free petroleum 

ther to which known amounts of benzene were added. A quan- 
titative yield of dinitrobenzene was obtained. 

Because small amounts of dinitrobenzene remain in solution in 

he petroleum fraction after nitration, a corresponding correction 
was made in calculating the benzene content. The correction was 
found by determining the solubility of m-dinitrobenzene in the frac- 
tion. ‘This was in every case less than 0.2 g dinitrobenzene per 100 
g of hydrocarbons. 

After each nitration the melting point of the dried nitro product 
was determined. All of the melting points were close to that of 
m-dinitrobenzene and by recrystallization of the material they could 
easily be raised to the accepted melting point of this compound 
(90.7" (3. 

From the weight of the nitro products the amount of benzene in 
each fraction was calculated. The results of the nitration are shown 
in Figure 2, Graphs IJ and IIT. 


3. EXTRACTION OF BENZENE 


From some of the fractions benzene was removed by extraction 
with aniline at about 0° C. This was done as illustrated in Figure 3.'* 
The apparatus consists of a narrow brass box,A, (80 em long 8 em high 
and 4 em wide) divided into 20 sections by means of 19 brass plates, 
B. A shaft, C, with 17 propellers, D, passes through the entire 
length of the apparatus. The shaft is rotated by means of a motor, 
M. The plates are provided with narrow vertical slits (not shown in 
the figure) to allow liquid to pass from one section into another. 

The apparatus is operated on the counter-current principle, aniline 
entering at F and leaving at G. The hydrocarbon to be extracted 
enters at H and overflows through J. The two liquid layers are thor- 
oughly mixed by the propellers (about 300 r. p.m.) so that equilibrium 
is established betw een the phases. The apparatus is covered with a 
glass plate, J, and is placed in a cooling bath as shown. — A sight 
glass in each end of the apparatus (not shown in the figure) indicates 
the height of the interface between the two layers. The right-hand 
side of the apparatus may be lowered or raised by means of the level- 
ing screw, K. This permits regulation of the overflow ratio of the two 
liq uids. As indicated by the broken line, the hydrocarbon phase 
may be recirculated for futher extraction. This can be done either 
continuously by means of a circulation pump or intermittently by 
pouring. 





1 AC nieaiiallaien nt is me " to F. W. Rose, jr., for the drawing of Fig. 3 
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» After extraction the hydrocarbon phase is washed with hydro- 
I hloric acid and with water until it is free from aniline. The aniline 
‘phase is distilled and the first distillates, containing the extracted 
‘hydrocarbons, are freed from remaining aniline in the same way. 
; The extraction method was applied to the constant boiling mixture 
of benzene and n-hexane. In this experiment 1,600 ml of the 68° to 
69° C. cut (n?2=1.395) was extracted once as described. As a 
Fresult 2 fractions were obtained, 1 of 200 ml with n?=1.419 and 1 
Fabout 1,300 ml of material with a refractive index nf =1.391. 

By subjecting the 200 ml fraction to equilibrium melting © a small 
purified sample of benzene (n= 1.501) was obtained (See p. 877.) 

Removal of the benzene by extraction with liquid SO. was also 
iried. The apparatus and technic developed by R.'T. Leslie, of this 
‘bureau were used. A 500 ml portion of the 68° to 69° (nf =1.395) 
cut was split up as follows: 7 ml with nj =1.492, 20 ml with np = 
1.427, and 425 ml with n? = 1.392. 


4. ISOLATION OF N-HEXANE 


In order to isolate the n-hexane, the fraction boiling at about 69° 
(. was treated with chlorosulphonic acid.’ This acid reacts com- 
paratively slowly with the normal hydrocarbons, while the isomers 
and cyclic compounds are attacked vigorously and destroyed. A 
1,000 ml sample of the 69° to 70° C. cut was agitated with an equal 
volume of chlorosulphonic acid (technical) at 40° to 60° C. for about 
12 hours. The acid sludge was then replaced with fresh acid and 
the agitation continued for another 12 hours at about 60° C. After 
separation from the heavier acid sludge, the hydrocarbon layer was 
washed in a separatory funnel with 5 per cent NaOH solution and 
with water finally dried over CaCl,. After this treatment the material 
' consisted of almost pure hexane which was distilled over sodium in a 
\)-plate column.'* The result of this distillation is shown in Table 1. 


TABLE 1.—Dhistillation of n-hexane 


| 
Refractive 
index nj) 


Boiling 
point 
corrected 


Volume of 


Tracti N 
Fraction No. fraction 


ml | °C. 
3 .6 





1, 3750 





. 3751 
59 





"Residue. 

Fractions Nos. 4 and 5 were mixed, after which the freezing be- 
havior and the physical constants of the sample were determined. 
As shown later, these constants agree well with those of the best 
samples of synthetic n-hexane prepared by other investigators. 





Ml. M. Hicks, B. 8. Jour. Research 2, p. 484; 1929. 
YU, Aschan, Ber., 31, p. 1801; 1898. 
'A. F. Shepherd and E. L. Henne, Ind. Eng. Chem., 22, p. 356; 1930. 
See footnote 11, p. 870. 
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III. DISCUSSION OF THE RESULTS 


1. BENZENE 


Graph II in Figure 2 represents the concentration of benzene (i 
per cent) in the different cuts. The fraction boiling between 6¥ 
and 69° C. has the highest concentration of benzene, 6.5 per cent, 

Graph III in Figure 2, representing the distribution of the benzene. 
in the different fractions, shows that the greatest amount of benzene 
(433 g) is present in the large fraction boiling between 69° and 70°C. 
Practically all of the be nzene of this petroleum is found in the frac. 
tions boiling between 65° and 75° C. 
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Fiaure 4.—Freezing curve for benzene 


It is interesting to note that while the peak for Graph ITT (mass o! 
benzene in fraction) is exhibited by the fraction boiling between 6% 
and 70° C., Graph IT (cone entration of benzene in the fraction 
reaches its maximum for the fraction distilling between 68° and 69° C. 
This is to be expected in view of the fact that the fraction 69° to 70° 
consists mainly of n-hexane and the fact that the constant-boiling 
mixture of n-hexane and benzene has a minimum boiling point. Ha 
this been a maximum boiling point, the peak of Graph II would have 
fallen to the right of the peak for Graph III. 

Another interesting feature of Graph II is the fact that while the 
right branch rapidly falls toward the abscissa, the left branch falls 
off at first rapidly and then exhibits a rather irregular behavior. For 


Bor 
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Foncentrations of benzene below 1 per cent the slope of this left side is 
nonuniform and only slowly approaches the abscissa. The ex- 
; planation of this may be found in the fact that the isomers of 
‘n-hexane’ also form constant-boiling mixtures with benzene. 

' A total yield of 1,530 g of benzene was obtained. Based upon the 
B original material, 600 gallons of crude oil (specific gravity 0.9), it can 
She concluded that the crude oil contains not less than 0.08 per cent of 
'benzene. 
; Tasie 2.—Physical constants of benzene 





Normal | Freezing —" 


Benzene boiling | point 20 


Density 
20 
point | (in air) | ay | a4 





°c. | °C, 
Salata Tih SRNR ols ctiamnnt denn amavcinsan sya am 80. 1 | 5.3 1.501 | 0. 878 
Pile DERUNID << ccasaatnepneds ce éndscdnous snes ausatdwesnueees 1 80. 2 | 15.5 Miosianes 2. 878 











J.T 23, p. 751; 1926. 
1Int. Crit. uD abies, 1, p. 198. 


The physical constants of the purified sample of this benzene are 
compared with those of pure benzene in Table 2. The time-tem- 
'perature freezing curve of the sample of benzene isolated i is shown in 
‘Figure 4. As indicated by an arrow, the benzene froze solid after 
eight minutes, at which time the stirrer had to be stopped. The 
pbenzene near the thermometer bulb then cooled more slowly than that 
near the walls of the tube, and this caused the curve to fall gradually 
at first, then to drop suddenly. The dotted curve indicates the ap- 
‘proximate slope, which would have been expected, if the stirrer could 
have been kept in operation. The initial freezing point of 5 
corresponds to an estimated purity of 99.8 mole per cent. 


2. n-HEXANE 


The physical constants of the n-hexane isolated from petroleum are 
shown below. 
TABLE 3.—Physical constants of n-hexane 





| Refractive | Density | 


n-Hexane ony ; | Aniline 
760 mm i n20 | d; point 





——————| 





n-Hexane isolated from petroleum 
Synthetic 


1, 3750 | 0. 6599 
. | 
B. J. Mair ! (99.94 mole per cent) } 


7 
pe ¢! 
Timmermans ? 80 














'B. J. Mair, unpublished work. 
‘J. Timmermans, J..chim. Phys., 24, p. 411; 1928. 
'G, Chavanne, Bull. soc. chim. Belg., 31, p. 331; 1922. 


A sensitive criterion for purity of a hydrocarbon is its behavior 
during freezing.” A comparison was made by B. J. Mair * between 
the er behavior of a sample of very pure synthetic n-hexane 


nee 


‘J. H. Bruun and M. M. Hicks-Bruun, B. 8. Jour. Research, 5, pp. 933-42; 1930. 
ne W. Washburn, Ind. fan Chem., 22, p. 985; 1930. 
1 See footnote 1, Table 3 
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prepared by him, and that of the n-hexane obtained by us. 1, 
c omparison was made i in the same apparatus, on the same day, and by 
using exactly the same technic in both cases. The time- -temperatuy 
freezing curve of the n-hexane from petroleum is shown in Figure ; 
The initial freezing point of the sample (— 95.67° C.) is 0.33° C. lowe 
than that of the synthetic hexane (—95.34° C.). If the value 34, 
cal./g.” is used as the heat of fusion of hexane, a depression of ¢\ 
freezing point of 0.33° C. corresponds to a purity of 98.3 mole per cent 
As it is not within the scope of this project to prepare hydrocarbons! 
the very highest purity, further treatment with chlorosulphonic a¢ii 
of the hexane sample was omitted. However, the n-heptane fry. 
tion * from petroleum, which has been treated for a longer time wit) 
chlorosulphonic acid, has freezing behavior (F. P. —90.68° C)}: 
which indicates a high degree of purity. 

Graph I in Figure 2 shows that normal hexane is present to a mue 
greater extent than any of its isomers. All of the fractions boilix: 
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Figure 5.—Freezing curve for n-hexane 


between 67° and 72° consist mainly of n-hexane. The minor consti: 
uents of these rections are probably 3- methylpentane,” benzene, 
and methyleyclopentane.” From a conservative estimate of t! 
fractions distilling between 67° and 72° C., we may conclude that: 
least 10 kg of normal hexane, or one-third of 1 per cent, is present! 
the crude petroleum. 

Graph I in Figure 2 gives us the following information as tot th 
quantities of material distilling at the respective boiling points of 
isomers of n-hexane: 














Mass of fraction (g) | Boiling range (°C.) 
i | 
DE tbh dnd huduicinimonenacbol 44 to 55 (boiling point of 2.2-dimethylbutane is about 50° C.). 
800....-.--------------------------| 55 to 58 (boiling point of 2.3-dimethylbutane is about 58° C.). 
0) WEE Ee CCR Eee ee ett. | 58 to 61 (boiling point of 2-methylpentane is about 60° C.). 

Seah eS | 61 to 64 (boiling point of 3-methylpentane is about 63° C.). 

” Parks and Todd, Ind. Eng. Chem., 21, p. 1236; 1929, 25 See footnote 7, p. 870. 

* To be discussed in a later paper. 26 See p. 873. 

% See footnote 20, p, 877. 27 To be discussed in a later paper. 
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Without further investigation these data do not justify any state- 

| ment concerning therelative amountsin which the isomersof hexane are 
present. In a previous investigation on the same petroleum by the 
authors, however, it was found that all the isomers of n-hexane 
(except the 2.2-dimethylbutane which was not found) were found 
in the fractions boiling between 54° and 75° C. Based upon this 
fact as well as the Graph I in Figure 2, we might draw the conclusion 
that the isomers of n-hexane are present in the following approximate 
quantities: 

3-methylpentane 0.18 (or less) per cent of the crude oil. 

2-methylpentane 0.09 (or less) per cent of the crude oil. 

2.3-dimethylbutane 0.04 (or less) per cent of the crude oil. 

The relation between the content of hexanes and the percentage of 
benzene in a given petroleum is important in that it determines the 
temperature at which the bulk of the benzene will distil. In the ab- 
sence of the hexanes, the benzene will normally be present in the 
fractions which boil around 80° C. When relatively small amounts 
of the hexanes are present some of the benzene will distil (as a constant 
boiling mixture) with the hexanes below 70° C. after which the bulk 
of the benzene will come over at about 80° C. If the hexane-ben- 
zene ratio is large most of the benzene will be found (as an azeotropic 
mixture) in the distillates boiling at about 70° C. 
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CALORIMETRIC DETERMINATIONS OF THERMAL PROP- 
ERTIES OF METHYL ALCOHOL, ETHYL ALCOHOL, 
AND BENZENE 


By E. F. Fiock, D. C. Ginnings, and W. B. Holton 







ABSTRACT 






The method and calorimetric equipment previously used in the determination 
of thermal properties of saturated water and steam were applied, without essential 
modification, in an investigation of corresponding properties of methyl] alcohol, 
ethyl alcohol, and benzene. Precautions were taken to obtain these liquids in 
a very pure state. Determinations were made in the range from 40° to 110° C., 
yielding values of latent heats of vaporization and changes in heat content. 
Empirical equations fitted to these properties aid in the interpolation, extropola- 
tion, and application of the data. An estimate is made of the precision of the 
final results. A skeleton table of the more commonly used thermal properties 
of the saturated liquids and vapors has been prepared. Comparisons with the 
values of previous investigators are included. 
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I. INTRODUCTION 


The calorimetric equipment which was designed and used at the 
National Bureau of Standards ! for the measurement of thermal prop- 
erties of saturated water and steam became available for a limited 
survey of similar properties of other liquids upon completion of the 
initial steam program. Reliable thermal data on methyl alcohol, 
ethyl alcohol, and benzene were needed to supplement calorimetric 
measurements on heats of combustion which are now in progress at 
this bureau. On account of the physical and chemical properties of 
these liquids, the apparatus could be used without essential modifica- 
tion. They can be obtained in comparatively pure form, and their 
purity can be estimated by rather simple physical measurements. 


1B. 8. Jour. Research, 5, p. 411; 1930. 
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II. METHOD, APPARATUS, AND PROCEDURE 


Each of the three organic liquids already mentioned was subjected 
to two distinct calorimetric processes. The first of these consisted in 
heating a known mass of the fluid from one saturation temperature to 
another, observing the total energy added and the temperature rise. 
Such measurements yield essentially the change in heat content or 
enthalpy of the calorimeter and its contents over the observed tem- 
perature interval. The second process consisted in removing satu- 
rated vapor from the calorimeter at practically constant temperature. 
These measurements yield essentially the value of the latent heat of 
vaporization at the temperature of the experiments. 

To obtain values of heat content and of latent heat from the meas- 
urements as they were performed, a correction is necessary to take into 
account the change in the ratio of the masses of liquid and vapor in 
the calorimeter from the beginning to the end of an experiment. This 
correction has been shown in the published theory of the method ? to 
be the same for both types of experiment, and, where desirable, can 
be determined by a third type of experiment in the same apparatus. 
For temperatures near the normal boiling points this correction is 
small compared to the quantities with which it must be combined, 
and its value can be calculated from specific volume and vapor pres- 
sure data more accurately than it can be measured calorimetrically. 

The present series of experiments was undertaken primarily for 
obtaining values of latent heats, but, since it is not possible to start 
and finish each latent-heat experiment at identically the same tem- 
perature, the heat-capacity experiments were necessary for the evalua- 
tion of end corrections. The values of heat content obtained from 
the present work are probably as precise as any now available, but 
are not as well established as would have been desirable had more 
time been available for this work and had there been a specific need 
for more precise values. 

The theory of the method, the apparatus, and the procedure were 
essentially the same as those developed in the work on steam. The 
two reports to which reference has already been made are given in 
considerable detail, and it is deemed unnecessary to include a lengthy 
repetition here. 

The apparatus is built to provide for accurately controlling, observ- 
ing, and accounting for the amount, change in state, and change in 
energy of a sample of fluid. A single calorimetric equipment yields 
results which establish values of latent heat and changes in heat 
content of the saturated fluids in the range covered by the experiments. 
A schematic representation of the calorimeter is given in Figure 1. 

A quantity of fluid, part liquid and part vapor, is contained in a 
metal calorimeter shell (C). The liquid is circulated rapidly over an 
electric heater (#7) and about the interior of the shell by a small cen- 
trifugal pump (P), to distribute heat and promote close approximation 
to thermal equilibrium. The sample at some chosen saturation state 
is heated to some other chosen saturation state, or else is withdrawn 
as saturated vapor. Energy is added and measured electrically. __ 

Outlets with valves provide for the introduction and withdrawal o! 
both liquid and vapor. The line at the top (L’) is for removing vapor, 
and is provided with a throttle valve (7) for controlling the rate o! 


2 B. 8. Jour. Research, 4, p. 609; 1930. 
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vapor flow. Following the throttle valve in this line are three other 
valves (S) for diverting the flow into any one of three receivers. 
These detachable receivers, suitable for weighing, are connected to 
the outlets to hold the samples of fluid transferred. 
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Figure 1.—Schematic representation of the calorimeter 























The calorimeter is surrounded by an envelope (£) and the inter- 
vening space is filled with nitrogen. The envelope temperature is 
controlled and regulated by an oil bath which surrounds the entire 
apparatus. 
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Absolute temperatures are measured by means of a platinum resist- 
ance thermometer and small temperature differences by means of 
thermoelements. 

In the present series of experiments the operators had the advantage 
of previous experience with the calorimeter. This proved useful not 
only in the actual manipulation, but also served as a guide in the 
simplification of the procedure. 

During the steam research it was demonstrated that the pump 
energy was independent of the height of the liquid in the calorimeter, 
within the limit of error involved in its experimental determination, 
For this reason all measurements of pump energy with organic liquids 
were made with low fillings of the calorimeter. 

The experimental values of the latent heat of water were shown to 
be independent of the rate of withdrawal of vapor. Hence the present 
series of vaporization experiments was performed at a single con- 
venient rate. The work on steam, it is believed, demonstrates that 
only saturated vapor is withdrawn under normal operating conditions. 

The quantity of energy passing by conduction, convection, and 
radiation between the calorimeter and its envelope is independent of 
the fluid in the calorimeter. The values of this thermal-leakage 
coefficient were determined during the work on steam, and were not 
remeasured for the present series. Because of the favorable range of 
temperature covered in these experiments and the greater skill which 
had been acquired by the operators, the actual values of the thermal 
leakage were so small in all cases that they could have been entirely 
neglected without appreciable effect on the final results. 

The rate at which mechanical energy was supplied by the circulat- 
ing pump was measured for each of the three compounds at several 
different temperatures. Curves through the experimental points 
made it possible to interpolate at the places where there were no 
determinations. 

The resistance thermometer was recalibrated late in the steam pro- 
gram and showed no significant change from the previous calibration. 

The potentiometers and electrical standards have been calibrated at 
frequent intervals at this bureau, and the thermal quantities are 
reported in international electrical units referred to the standards as 
maintained here. 

III. NOTATION 


The notation employed in this report is the same as that used by 
Osborne * in the published theory of the method. For convenience 
his notation will be repeated in part. 

The quantity obtained directly from a heat-capacity experiment 1s 
called a and from a vaporization experiment y. The correction which 
takes into account the change in the proportion of liquid and vapor 
masses during an experiment is called 8. Each of these quantities, 
when referred to in this paper, will be used specifically; that is, refer- 
ring always to 1 g of the fluid. 

Other quantities may be defined as follows: 
H=heat content or enthalpy per unit mass of saturated liquid 

(e+7u). 
H’=heat content or enthalpy per unit mass of saturated vapor 
(e’ +7u’), 


8 See footnote 2, p. 882, 
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e=internal energy per unit mass of saturated liquid. 
.«’=internal energy per unit mass of saturated vapor. 
L=latent heat of vaporization per unit mass. 
g=temperature on the international scale. 
@=temperature on the Kelvin scale. 
7+ =saturation vapor pressure. 
u=specific volume of saturated liquid. 
uw’ =specific volume of saturated vapor. 
Q=total energy added to the calorimeter in any experiment. 
Z=a constant of the calorimeter depending on its temperature and 
the liquid it contains. 
M=total mass of fluid in the calorimeter. 
As shown in the published theory, the following thermodynamic 
relations may be used to reduce the quantities actually observed to 
those more commonly employed. 


H=a+ 86 (1) 
Bree 
; u 
u’—U 


/, 


y=L 


U 
=L——— =ou= 
B=L u’—U “710 


A’=H+L 
Q=Z+Ma (6) 
ae w|e-lZ + Mra} + f2(H— H,)dM (7) 


The values of H and L were obtained from the experimental values 
of a and y and the values of 8 calculated from other data with the aid 
of the second form of equation (4). 

The absolute value‘of H can not be determined, but any arbitrary 
datum may be chosen within the temperature range covered by the 
experiments, and values of H at other temperatures in the range can 
be established relative to the datum. Adoption of a datum for H 
simultaneously fixes a reference value for a. 


IV. PREPARATION AND PHYSICAL PROPERTIES OF 
COMPOUNDS 


1. METHYL ALCOHOL 


Mallinckrodt’s methyl alcohol, specially ordered as chloroform free, 
was distilled by Dr. J. H. Bruun, of the chemistry division of this 
bureau, in a nonsiphoning, 30-plate, bubbling-cap still * to remove 
water and other impurities which can be separated by fractional 
distillation. The middle fraction of 1,000 cm* was used from a charge 
of 2,500 em’ which was introduced into the still pot. The distillation 
was performed in an atmosphere of dry carbon dioxide at approxi- 
mately atmospheric pressure. The distillate was collected in a bottle 
filled initially with dry air and was stored in a dessicator. Provision 
was made for transferring the liquid without further contact with the 





‘Ind. Eng. Chem., Anal. ed., 1, p. 212; October, 1929. 
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air to a second still, which could be used both for refluxing and dis. 
tiling. Here the sample was freed from dissolved gases by repeated 
refluxing, cooling, and evacuation. After thus pumping off the dis. 
solved gas, the sample was distilled under its own vapor pressure at a 
temperature of about 50° C., and collected in evacuated silver 
containers. From these it was transferred to the calorimeter. 

The specific gravity of the gas-free product, before its use in the 
calorimeter, was measured by the volumetric section of this bureau, 
To avoid contact with the atmosphere the alcohol was introduced 
into an evacuated picnometer through an evacuated line. The 
volume of the picnometer was about 87 cm®, and the samples used 
therein were always discarded. The observed specific gravity at 
20° C. compared to water at 4° C. was D?=0.79133. According to 
the International Critical Tables, vol. 3, p. 27, Klason and Norlin give 
as their value for absolute methyl alcohol D? =0.79134. 

An accurate combustion analysis of the original material was made 
by Dr. F. D. Rossini, of the chemistry division of this bureau. This 
analysis showed that the carbon-hydrogen ratio was that required by 
the formula CH;,OH, within the limit of error of the analysis, which 
was + 0.03 per cent. 

After the completion of a vaporization experiment the sample was 
in four containers, and had to be reassembled before being weighed 
into the calorimeter again. This was accomplished by transfer 
under its own vapor pressure through an evacuated line, without any 
contact with the atmosphere. Such a procedure made possible the 
performance of the whole series of calorimetric experiments on an 
original sample of about 350 g. 

At one stage in the experiments with methyl alcohol a leak devel- 
oped in the lower tube leading to the calorimeter. This permitted 
contact of a portion of the sample with the atmosphere. This por- 
tion was again freed from dissolved gas, but contained a small amount 
of water. For the remainder of the experiments the entire sample of 
alcohol was reunited, since the amount of water which could have been 
taken from the air was known to be small. 

After the completion of the calorimetric experiments the specific 
gravity of the assembled sample was again measured and the value 
obtained was D*?=0.79142. 

The increase in specific gravity was thus 0.00009, corresponding, 
if the impurity was all water, to about 0.026 per cent of water by 
weight. This fraction of water would produce no measurable effect 
on the values of a and would correspond to an increase in y of only 
about 3 parts in 10,000. Since this is well within the experimental 
uncertainty involved, it was deemed unnecessary to repeat the experi- 
ments with another sample of alcohol. 


2. ETHYL ALCOHOL 


Through the courtesy of Dr. F. M. Hildebrandt, the United States 
Industrial Alcohol Co. furnished an especially selected sample of 
ethyl alcohol (192 volume proof) having a very low aldehyde content. 

Most of the water present was removed by refluxing and distilling 
from an excess of freshly ignited lime at atmospheric pressure. 
About 2,500 em* of alcohol, whose specific gravity indicated a water 
content of 0.15 per cent by weight, were prepared by this method. 
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This 2,500 cm* sample, together with 40 cm*® of Timmermans’ 


}enzene, which will be described later, was introduced into a still 
.jmilar to that used for fractionating the methyl] alcohol. 


The remainder of the water was removed from this mixture by 
Yractional distillation. The first fraction to distil consists of a con- 


' stant boiling mixture of benzene, alcohol, and water. When all the 
| water has been removed, a constant boiling binary mixture of benzene 
| and alcohol continues to distil as long as any benzene remains. Since 
"the indices of refraction of benezne and alcohol are so widely different, 
| this property was used as a criterion for the complete removal of the 


‘benzene. About 250 cm® were rejected after the index of the distil- 


F late became constant. The next 900 cm* which distilled were used 


as the best product of the fractionation. 

The method of collecting the sample and of removing dissolved 
vases from it was the same as that described for methyl alcohol. 
The specific gravity of the gas-free product was found to be D?= 
0.78946. According to the International Critical Tables, vol. 3, p. 
27, Osborne, McKelvy, and Bearce give as their best value for 
absolute ethyl alcohol D #? =0.78934. If water was the only impurity 
the initial sample contained 0.038 per cent by weight. The specific 
gravity of the aleohol which had been used in the entire series of 
observations was D2? =0.78949, indicating 0.047 per cent of water by 
weight. The effect of this quantity of water on the measured values 
of latent heat could not have been greater than 0.6 joule per gram, 
or 1 part in 1,400 of the total. 


3. BENZENE 


About 1,000 em* of benzene prepared by Timmermans and Martin® 
at the International Bureau of Physico-Chemical Standards at 
Brussels were refluxed for three hours with about 15 g of a 50-50 mol 
per cent sodium-potassium alloy. This alloy was furnished by Dr. 
0. C. Bridgeman, of this bureau, who had prepared it in very pure 
form for use in determining quantitatively the water content of 
gasolines. 

The benzene was introduced into an apparatus which could be used 
both for refluxing and distilling. After it had been frozen and the 
permanent gas present removed by pumping, the alloy was intro- 
duced. After refluxing for a time the benzene was again frozen, and 
the gas evolved on boiling was pumped off. This process was re- 
peated four times, the amount of noncondensible gas coming off the 
last time being practically nil. 

_ After three hours of refluxing, during which numerous droplets of 
liquid alloy were in vigorous circulation within the liquid benzene, 
the benzene was distilled under its own vapor pressure at a tempera- 
ture of about 50° C., and collected in evacuated silver containers. 

_ The specifie gravity of the air-free product was not measured before 
the calorimetric experiments because a sufficient quantity could not 
be spared. After the measurements were completed, specific gravi- 
lies were measured on two samples which had been distilled from the 





‘J. chim. phys., 23, p. 18; 1926, 
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calorimeter, and on the residue left in it after a series of Vaporizatio, 
experiments. ‘The values of D ?? obtained were 0.87906 and 0.879(7 
for the distillate and 0.87906 for the residue. 

The recorded determinations of the specific gravity of benzene 4 
20° C. vary from 0.8787 to 0.8795, with those which seem mop 
reliable averaging about 0.8790. From specific-gravity measur. 
ments alone it is difficult to estimate the purity of benzene for wan; 
of a well-established reference value. 

The product used in these experiments was prepared from the source 
which, according to our advices, was the best available. ‘The water 
content was probably reduced to well below 0.01 per cent. Th 
presence of homologs of benezene in minute quantities would have 
produced no measurable effects on the thermal properties observed. 


V. RESULTS OF MEASUREMENTS 


In the following record of these measurements the results have 
been assembled in the form of tables. An explanation of the mech. 
nism of the reduction of the quantities observed to the properties 
desired is given in the report on the work on steam.® 

Three complete sets of a experiments were performed with a larg 
mass and three with a small mass of fluid in the calorimeter. Th 
experiments did not start and stop at exactly the even 10° temper- 
tures so the measured values of Q were corrected to the even intervals 

; ‘ AQ 
for convenience. From curves of aN plotted against temperature 
the small corrections to the even 10° points were determined, and 
these corrections, applied to the measured values of Q, yielded the 
values of Q]} corrected to the even temperature intervals. 

Choosing for example the interval 40° to 50°, the value of the 
change in @ for the liquid was computed for this interval as follows. 
For each observed value of Q]%§ there is an observed mass M of fluid 
in the calorimeter system. To determine the values of the two 
constants Z and a@ in the equation. 


Qi = Z}i8 + Mal’? 


the method of least squares was employed. As an indication of the 
consistency of the results, the value of Z]3§ so obtained was used for 
the computation of the value of @ indicated by each individual 
experiment. The tabulated values of a, are the values obtained 
by the method of least squares, and the other values of a were 
computed for the separate determinations as just described. The 
results of the a experiments are shown in Table 1. 


6 See footnote 1, p. 881. 












40-5 






BFiock, Ginnings | Thermal Properties of Organic Liquids 889 


{olton 
TaBLE 1.—Results of a experiments 


2 uU 2 
aft-[a-rti); 


ee a 








| 
Methy! alcohol | Ethyl alcohol | Benzene 
| 
Temperature interval (°C.) “7 ahi 
i : | . 2 | ge ae 
al} a|i—am | al; a|{—am al; | a|{—am 





9 
. 025 



















































































Int. J./9 Int. J./g | Int.. Int. J./9 Int. J/g | Int. J/9 
26. 075 —0. 005 26 —0. 012 |-- : ; 
26. 075 —. 005 26. 030 | —. 007 ee 
: 26. 089 . 009 | 26. 057 |  " Rees Tssibiseds 
40-50..-----eeneonnnnnnrrnnrnr= 26. 087 . 007 26. 023 | —.014 bs 
26. 068 —. 012 | 26. 061 | 024 
26. 084 . 004 26. 028 | —. 009 |. 
pits cos oaeeeodeaanee ae eee 26. 037 
26. 799 —019| 26.90|  —.016| 18.340 0. 037 
26. 847 029 26. 918 —. 004 | 18. 290 —.013 
bs 26. 808 —.010 | 26. 944 . 022 | 18. 279 —. 024 
50-60..----n-neneennonennconce= 26. 807 —o11 | 26. $08 | —.014 | 18. 280 | —.023 
26. 816 | —. 002 26. 920 | —. 002 18. 319 016 
O11 26. 940 | 018 18. 310 | 007 
at 26. 922 Ce 4 ee 
021 | 27. 986 004 18. 546 | 010 
—.034| 27. 967 —. 015 18. 554 000 
0-90 O11 | 27. 992 “010 18. 546 —. 008 
a —. 026 | 27. 990 . 008 18. 582 . 028 
007 | 27. 960 —_ 022 18. 513 —. 041 
018 | —-27. 995 | “013 18. 568 | 014 
| eal 27. 982 18. 554 |_- 2 
—.004 | 29. 045 006 18. 782 014 
“O11 | 29. 033 —. 006 18. 764 —. 004 
“n-R0) —. 008 | 29. 026 —, 013 18. 757 ~. O11 
. 050 | 29. 017 | —. 022 18. 726 —. 042 
—. 036 | 20. 073 034 18. 792 024 
—.015 29. 041 | 002 18. 785 ‘017 
ek | 29. 039 | Pt 16,/768' 1..---=< 
—003 | 30.026 | 005 19. 012 . 009 
| 005 | 29. 996 —. 025 18. 998 —. 005 
“ie |) 28.868 —. 001 | 30. 042 021 18. 999 —. 004 
a aa cia aa } 28. 892 | . 023 | 29. 999 —. 022 18. 968 —. 035 
28. 854 | —015 | 30.058 ‘037 19. 022 019 
98. 862 | —007| 30. 607 | —.014 19. 019 016 
ipa eeecunawasd= 98.869 |...---------| 30. 021 |----------- cic: 4 ee 
39.704;  -«. 034| +~—«31. 316 021 19. 284 024 
29.651 | —. 019 | 31. 309 “014 19. 262 . 002 
29. 656 | —.014 31. 259 —. 036 19. 235 —. 025 
O- “ 
0-100..--nconeneon-nonennn-e= 29. 693 | ‘023 | 31. 327 "032 19. 295 035 
29. 667 | —. 003 | 31. 288 —. 007 19. 239 —. 021 
29. 650 } —. 020 | 31. 269 —, 026 19. 247 —.013 
Ga™ ..ce-ceceeeeeaee-ee- 29. 670 |. ..| 31. 295 19. 260 
~~ 30. 338 | -020| 32. 208 —. 021 19. 501 004 
30. 313 | —. 005 | 32. 220 —. 009 19. 501 004 
100-110 30. 302 | —.016 | 32.258 ‘029 19. 488 —. 009 
piacere nea sn eames 30. 340 | 022} 32. 209 —. 020 19. 504 007 
30. 302 5 002 19. 490 ni OUT 
30. 311 ‘017 19. 496 —. 001 
et REPEC Ri Pe eee 30. 318 7 19. 497 





Before continuing further with the reduction of the data from the 
heat-capacity determinations, the results of the experiments will be 
tabulated. The values of y were calculated from equation (7), in 
which 1; is the value of at the initial temperature and AM represents 
the mass of vapor removed. The quantity (Z-+ Mz] i is that energy 
quantity involved in a change in calorimeter temperature from the 
initial to a different final value. The values of Z and a are available 
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from the results of the @ experiments, and that of M,, the mags ». 
maining in the calorimeter at the completion of a y experiment, \ 
known from the weighings. The quantity /,°(H7— H1)dM represen 
the energy involved in a change in the temperature of the vapy 
leaving the calorimeter during an experiment. It is negligibly sma 
under thesteady temperature conditions which prevailed during the pra 
ent work. It is mentioned only for consistency with the exact theory 

The results of the y experiments are assembled in Table 2. Tip 
corrections of the observed values of y; to the desired temperature 
dy 
do 
small because the values of 0, were so near to the temperatures 
recorded. 


involved graphical determinations of The corrections are very 


TaBLE 2.—Results of y experiments 


Iu’ 
, 


7 


u’—Uu 





ee —_—$—$————__. 


Methyl alcohol Ethyl aclohol Benzene 
Temperature 4 he Temperature, : We Temperature, 
°C.) = (°C.) ron (° 0.) 
| 


J/g | Int. J./9 Int. J./g Int. J./9\| 
2 0.1 886.5! —0.4 | 
| 


) 
885.6 | —l.% 


886, | —, 4 50. 





119 


1129. 3 
1101. 6 
1102. 2 
1102. § 


1101. 5 





1102. $ 
1102, § 
1102. 7 
1102. 9 | 


1102. 
1102. § 
1103. 
1103. 








1102. 5 


1049. ¢ 
1048. 8 | 
1049. 4 | 
1050. 
1049. 
1048. 6 
1050. 7 


1049. 5 


1001. 7 
1000. 

1001. § 
1001. : 


1002. 
1000. 
1002. 
1002. ¢ 





1001. 
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‘3 To reduce the values of a and y to the desired quantities H and 
“ls ay it is necessary to evaluate the correction term 8 throughout the 
"por i ange of the experiments. Since 8 is small in this range, it is believed 
Mal Bhat its values can be calculated more accurately from v apor pressure 
_ a nd specific volume data than they could be measured calorimetri- 
| ally. 





Ures y All terms in the right-hand member of the equation B=Ou 


very were evaluated for each of the substances concerned from the results 


uns (of measurements made by Young.” Va 


_ by differentiation of his vapor pressure equations. ‘Table 3 shows the 
quantities taken from the work of Young and the values of 6 com- 
puted therefrom at those temperatures w here they are needed in the 


| further reduction of the data. 


TABLE 3.—Calculation of values of B 


{1 liter- ieenaintainll 101.29 international vieaiaanl 


dr 
S=ON ss 


Ethyl] alcohol Benzene 

| uX10 | 8 | uxie | a | 2 
| ie do | 
} 

























dr 
de 


lues of = were calculated 























_— | Atm./de 


1 | 0.00111} 0.04] (1.1110) 

1. 2534 . 00200 | .07 1. 1243 . 00309 -10 
im. 

1. 27 . 00443 





Temperature Methy] alcohol 
| dr 
) S) ux10 
49 B 

°C, °K. | Liters/ /g Atm./deg.| Int. J./g 
0 | 273.1 | 1.2345] 0.00245 0.08 | 
10 | 283.1] 1.2489] .00416 15 
20 | 208.1} 1.2637] .00677 25 
25 | 2081] 1.2713} .00852 . 33 
30 | 303.1] 1.2790] .01061 .41 
4) | 313.1] 1.2948] .01606 65 
5 | 323.1] 1.3111] .02353 1. 00 
60 | 333.1) 1.3240] . 03355 1. 48 
64.7, 337.8] 1.3362] .03917 1.7 
70 | 343.1] 1.3405] .04651 2.14 
80 | 353.1] 1.3600]  . 06302 3.03 
%) 368. 1 | 1.3790 | . 08359 4.19 
100 | 373.1] 1.401 . 1088 5. 69 
110 | 3831] 1.424 . 1391 7. 59 
120 | 393.1] 1.449 .1751 2 7 | 
130 | 408.1] 1.477 . 2174 11 | 





1.2808 | . 00570 -22| 1.1520] 00688 |. 24 
1.2953 | . 00907 | .37| 1.1664] .00978 | . 36 
1.3104 | . 01394 | .59| 1.1813} .01350! .52 
1. 3261 - 02073 | 92 | 1.1967 | .01814 | .72 
| | 
1.3430 | .02989| 1.38| 1.2126] .02383| .99 
1.3609 .04192/ 201! 1.2278] .03064| 1.33 
1. 3791 05729} 2.87) 1.2436} 03865) 1.75 
} 
1. 397 .07648 | 3.99 1, 262 04794 | 2.26 
1. 417 . 0999 | 5. 42 1. 281 05857 | 2.87 
1. 444 -1280 | 7.27 1. 300 07059 | 3.61 
1. 473 -1612 | 9.69 1.321 7| 4.53 


2669 8A :13| 1.1380} 00469] . 16 




























g.| Int. J.Jg| Liters/g | Atm./deg.| Int. J./g 
(0. 0020) | (0. 08) 


17} 1.1449] .00568| 20 
| 














Table 4 gives the values of L computed from the mean experi- 
mental values, ym, by the relation L=y—8. 


TABLE 4.—Computation of latent heats 




































L=y—6 
) 
Temper- | Methy] alcohol Ethyl alcohol Benzene | 
: ; | 
meek EO TS L y | B L | vy | 6 | ZL | 
°C. Int. J/g | Int. J./g \Int. J.{g| Int. J.J | Int. J.Jg Int. J.Jg\ Int. J.Jq | 
50 1,123.3) 886.9 0.6} 886.3] 415.9 05! 4154 


wee eww ee we eww eeeoee 


Wee eee ene eee eee ne meee eww eee Cee eee ee Serene 











Proc. Roy. Dublin Soc., 12, p. 374; 1910, 
49527°—31——-9 
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To facilitate the interpolation and extrapolation of the values of 
latent heat, an empirical equation expressing J as a function of § 
was fitted to the observed values for each of the substances. The 
type chosen for these equations was the same as that which proved 
adequate to represent the experimental values for water. Each equ. 
tion requires that the latent heat approach zero at the critical ten. 
perature 6, at a negatively infinite rate, and that no real values exist 
above this temperature. Table 5 gives the values of L calculated 
from the equations which appear below the columns of figures. 


TABLE 5.—Formulation of latent heat values 
a a et cial | aed 
Methy] alcohol | 


Ethy] alcohol Benzene 


Temper- | 
ature 
6 


Latent heat of vaporization 


| i i . . } . ' 
|} Observed | Equation | Observed | Equation | Observed | Equation 


| 


Int. J./9 
1 (1, 200. 1) (447. 2) 
: | 7. ee | (936. 2) sa | (441. 2) 
, 175. 4) | . .-| (925. 4) (435. 0) 
Jeu io... (919. 6) | (431.9) 


| | 


Int. J./g Int. J./g Int. J./9 


(946. 0) 


Int. J./g Int. J./g 


, 161. 0) (913. 5) (428. 7) 


64. 7 
70 
80 
90 


100 
110 
120 
130 





, 145. 


, 128. 3 
, 109. 


. 3) 


“ 
é 


, 100. 4 


, O89. § 


855.1 


, 068. 3 


, 045. 2 


818.4 


, 020. ! 


994. 


(965 
(935 


. 9) 

) 

3 oS 
| 


777.8 | 


(900. 5) 
886. 5 
871. ¢ 


854.9 
837. ; 
818. 5 


798. 4 
777.0 
(754. 0) 
(729. 6) 





(422. 1) 
415.4 
408. 4 


401.3 
393. 9 
386. 3 


378. 
370. 
(361. § 
(353. 2) 

; 


1 Values in parentheses are outside the range of the experiments. 

For methy] alcohol 1 
L=—0.005 (240—6)?+42.60875 (240—0)+219 (240-0) 
For ethyl] alcohol \y 
L=—0.004067 (240—6)?+2.198 (240—0@)+165.83 (240-0) 4 
For benzene \Y 

L=—0.00056185 (290—6)?+-0.65028 (290—6)-+74.11 (290-0) 4 

The values of the changes in heat content as shown in Table 6 


follow directly from the data of Tables 1 and 3, through the equation 


H 2 Zz. 2 
"Siahar | Radke: Fa 
TABLE 6.—Computation of changes in heat content 
2 2 
H]j=«]} +8]; 


Ethy] alcohol 





| Methy] alcohol Benzene 
{Tem per-| 


ature | | 





9 -4+-aee | 

a |i | Bio | 40 5 Ph | H] 5 

Int. J./g |Int. J./g\ Int. J./9 
0.00 | 0.00 | 0. 


26.04 | 
52.96 | 


~ 
oe 





| 


Int. J/g) Int. J/g 


SERS 

= 

go gr 8 | 

BRSS Lasss 


80. 94 


109. 98 
140. 00 
171, 30 
203. 52 


| 
| 


od 
SEE 





# Eng po 
RLEB 

BREN Bho 
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Holion 


Empirical equations in the form of power series have been fitted 


to these values of H]j. These equations are shown in Table 7, 
together with a comparison of the experimental values with those 


calculated from the equations. 


TaBLE 7.—Formulation of neat content values 


Methy] alcohol Ethyl alcohol Benzene 





| 

| 

Tem- | 6 C7] 

perature H49 1) 5 
0 








| Observed | Equation |Observed | Equation | Observed | Equation | 


ae he ee | 


°C. | Int. J/g | Int. J./g | Int. J/g | Int. J/g | Int. J/g | Int. J/g | 
40 0. 00 0. 00 0. 00 Le ae 

50 26. 43 26. 47 26, 26 26. 26 0. 00 
60 53. 73 53. 75 53. 51 53. 57 18. 50 
70 81. 93 81. 91 81, 95 81. 99 37. 33 


| 
| 
80 | 111.04 111. 01 111. 62 111. 62 56. 43 
90 | 141. 06 141. 11 142. 50 142. 55 5. 86 
| 
| 








100 172.24 | 172.26 | 17492 | 174.86 5. 63 
110 204.55 | 204.54 | 208.57 | 208.63 73 























For methy]) alcohol 


62. ron 8 X J ’ o id 02 
Ho, [ 23.520( <7, ) +0.2583( 75) +0.010036( +7 6 
For ethy! alcohol 


9 ae | ’ o 9 ) “ 4 ) 3162 
Hp. 22.648( “) +0.3063( 10 +0.014233 10 6: 
For benzene 


uifi=[17.184( 7) +0.1008( ;7,)’+-0.002462( ;7,) "|e 


From the equations for H]? there were obtained values of a 


quantity 0,=5h, defined as the rate of change of the heat content 


with temperature under saturation conditions. Since the values of 
Hj? as recorded in Table 6 were measured under such conditions, 
equations for C, were obtained by differentiation of the equations 
already given for H]?. The derived equations and the values of 
C, calculated at several temperatures are given in Table 8. 


TABLE 8.—Computation of specific heat values 


dH 
Cr= Gg = Cy (approximately) 





| Methyl | Ethyl 
— | alcohol alcohol 
ire 


Cx Cr 








° 
Q 


Int. J/g. °C. | Int. J./9. °C.) Int. J./9.°C. 
. 608 2. 577 


to 


S8se2 xgss 
wep pyp 
errr res 
S388 822 




















For methy] alcohol . a 
On==2.35204+-0.5168 (4) +0.008012 ( tp ) 
For ethyl alcohol 
Or=2.2018+-0.06100( + ) +0.004270( +7)” 
For benze 
: Crm ri+-0.0006( #,) +a.or8n( &)? 
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From specific volume and thermal expansion data it is possible to 
transform the values of changes in heat content under saturation 
conditions to values of the changes under a constant pressure of | 
atmosphere. Such transformations have been made and equations 
for H}} at p=1 atmosphere have been fitted. 
oH 
z - 08 p=latm 
defined as the rate of change of heat content with temperature under 
a constant pressure of 1 atmosphere. Values of C, were calculated 
for each of the three liquids at temperatures below the normal boiling 
points. A comparison of these vation of C, with the recorded values 
of C, showed that the difference between the two was less than 1 part 
in 500 at every temperature. Since this difference is believed to be 
within the error involved in deriving values of specific heat from 
measured values of a, the recorded values of C, are not sensibly 
different from the values of C, as indicated by these experiments. _ 

The specific volumes of the saturated vapors can be calculated 
from the relation y = Ou’ = to the degree of precision of our knowledge 


dx 
of do 


From these it is possible to get derived equations for C, -( 


The values of wu’ thus computed are shown in Table 9. The 


dr ' 
values of qo #te calculated from Young’s vapor pressure equations 


and the values of y result from the present investigation. 


TABLE 9.—Computation of specific volumes 


uy =—— 


dx 


de 





| Temperature Methyl alchohol Ethyl alcohol Benzene 


| 


| de dr | de 
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Table 10 has been prepared to summarize the thermal properties 0! 
these three saturated fluids. The values of 7 are from the data of 
Young. The values of uw at temperatures below the normal boiling 
points are from density determinations as recorded in the Interms- 
tional Critical Tables, volume 3, pages 27 to 30. Values of u above 
the boiling points are those determined by oo The values of v 
are computed as shown in Table 9. The results of the present 
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investigation are extrapolated down to 0°C. and up to130° C. The 
value of H has been arbitrarily chosen as 0° at 0°C. for each liquid. 
For the benzene which freezes at about 5° above this temperature 
the datum H7=0 is for the undercooled liquid. The heat of fusion 
is therefore not included in the tabulated values. 

It is hoped that this table may prove useful as a working table for 
those who are concerned with the thermal properties of any one of 
these saturated fluids. 


TaBLE 10.—Thermal properties of methyl alcohol, ethyl alcohol, and benzene at 
saturation 


METHYL ALCOHOL 


























| | Specific volume 
Temper-| Vapor pres- ——- Poy bw Latent heat — ve 
ature 6 sure 3 a Z ” 
Liquid Va , H por 
} uxX108 Pat 
°é. Atmospheres Liters/g Litersa/g Int. J./g Int. J./9 Int. J./g 
0 0. 0389 1. 234 # 0.0 1, 200. 1 1, 200. 1 
10 .0713 1, 249 9. 97 23. 80 1, 188. 5 1, 212.3 
20 - 1251 1, 264 5. 851 48.17 1,175.4 1, 223. 6 
25 - 1632 1. 271 4. 544 60. 59 1, 168.4 1, 229.0 
30 - 2109 1. 279 3. 567 73. 18 1, 161.0 1, 234.2 
40 3427 1. 295 2. 249 98. 89 1, 145.3 1, 244.2 
50 . 5388 1.311 1. 465 125. 4 1, 128. 2 1, 253. 6 
60 . 8255 1,324 - 9816 152. 6 1, 109. 7 1, 262.3 
7 1, 219 1, 340 - 6755 180. 8 1, 089. 8 1, 270.6 
80 1. 764 1.360 ~ 4755 209. 9 1, 068. 3 1, 278. 2 
90 | 2. 493 1.379 . 3413 240. 0 1, 045. 2 1, 285. 2 
100 3. 451 1, 401 . 2495 271.2 1, 020. 5 1, 291.7 
110 4. 686 1, 424 - 1856 303. 4 994. 1 1, 207.5 
120 6. 252 1, 449 - 1401 336. 9 965. 9 1, 302. 8 
130 | 8. 209 1. 477 - 1089 371.6 935. 7 1, 307.3 








ETHYL ALCOHOL 








0 | 0. 0162 1.240 | 30.81 0.0 946.0 946. 0 
10 | . 0312 1, 253 16. 32 22.97 936. 2 959. 2 
20 0579 1, 267 9. 062 46. 63 925. 4 972.0 
| 25 | 0771 1. 274 6.877 58. 75 919.6 978. 4 
30 | 1028 1. 281 5. 221 71. 08 913. 5 984. 6 
40 | . 1756 1. 295 3. 132 96. 39 900. 5 996. 9 
50 | . 2892 1.310 1. 944 122.7 886. 5 1, 009. 2 
60 | . 4608 1. 326 1, 247 150.0 871.3 1, 021.3 
70 | .7117 1.343 . 8244 178. 4 854. 9 1, 033. 3 
| 80 | 1. 068 1, 361 . 5599 208. 0 837.3 1, 045.3 
| 90 | 1. 561 1. 379 . 3899 239. 9 818. 5 1, 058. 4 
| 100 | 2.227 1.397 . 2776 271.2 798. 4 1, 069. 6 
| 10 | 3. 105 1,417 . 2019 305. 0 777.0 1, 082.0 
12 4. 241 1.444 . 1494 340.3 754. 0 1, 094.3 
130 5. 682 1. 473 . 1123 377.3 729. 6 1, 106. 9 





BENZENE 





0 0. 0349 1,111 8. 805 0.0 447.2 447.2 
10 . 0598 1. 124 4. 980 17. 29 441.2 458.5 
20 . 0982 1. 138 3. 127 34. 79 435.0 .8 
25 . 1285 1,145 2. 442 43. 59 431.9 5 
30 - 1556 1, 152 2. 031 52. 52 428.7 -2 








40 - 2383 1. 166 1. 362 70. 51 422.1 492. 6 
50 - 3539 1, 181 - 9413 88. 75 415.4 504, 2 
60 - 5113 1, 197 - 6684 107.3 408. 4 508. 6 
7 + 7203 1, 213 . 4857 126. 1 401.3 527.4 
80 - 9916 1, 228 - 3607 145. 2 393. 9 539. 1 
90 1, 337 1, 244 . 164. 6 386. 3 550. 9 
100 1. 769 1, 262 - 2102 184. 4 378. 4 562.8 
110 2. 300 1, 281 - 1642 204. 5 370.3 574.8 
120 2. 945 1, 300 . 1300 225.0 361.9 586. 9 
130 8. 717 1, 321 - 1045 245. 8 353. 2 599.0 























































896 Bureau of Standards Journal of Research (Vole 
VI. DISCUSSION OF ACCURACY 


An estimate of the accuracy of the tabulated values of H and [ 
within the range of the experiments has been made on a basis similar 
to that outlined in the report on the properties of steam. The sy ste- 
matic errors in calorimetry are believed to be essentially the same jn 
both sets of experiments. 

By a study of the factors which enter into the measurements, 4 
figure was assigned to each which represents an estimate of the out- 
standing systematic error, after all known corrections have been 
applied. In addition to these, corresponding figures were deduced 
which represent the error in both the a and the y due to accidental 
errors in the measurements. An estimate was also made of the prob- 
able error in the calculated values of 8. The density determinations 
furnish a basis for estimating the error which may be present on account 
of impurities in the liquids investigated. 

These estimates were combined to obtain a value for the uncertainty 
in the final tabulated results. On this basis it is believed that more 
precise determinations on compounds of a higher degree of purity 
would yield values of H and L differing by less than 1 ‘part per 1,000 
from those of the present investigation in the range from 40° to 110°C. 


VII. COMPARISON WITH PREVIOUS WORK 
1. HEAT CAPACITY 


The results of calorimetric determinations of heat capacity have 
usually been recorded in terms of specific heats. For comparison of 
the present results with those of other observers the calculated values 
of specific heats can be compared directly, or, where equations are 
available, the values of changes in heat content can be calculated from 
integrated specific heat equations. The latter method would have 
been followed if the termperature ranges of the various experiments 
had overlapped sufficiently. However, the earlier work is confined to 
the range well below the normal boiling points, and hence covers regions 
in general, below that of the present investigation. 

The values of specific heat indicated by the present work have been 
extrapolated down to 20° C. for purposes of comparison. The extra- 
polated values of this derived quantity involve considerable uncer- 
tainties and should not be given a great deal of weight in the selection 
of the best values below 40°C. A graphical comparison of previously 
published values of specific heats with those of the present work 1s 
given in Figure 2. This summary includes the experimental results 
of Bose,® Kelley,® Williams, and Daniels,’ Forch,'!' Mills and Mac- 
Rae,” Dejardin,” Tréhin,’* and Andrews, Lynn, and Johnston.” 





8 Z. phys. Chem., 58, p. 585; 1907. 

®J, Am. Chem. Soc., 51, p. 180; 1929; and 51, p. 779; 1929. 
10 J, Am, Chem. Soc., 46, p. 903; 1924. 

11 Annalen der Physik, 317, p. 202; 1903. 

12 J, phys. Chem., 14, p. 797; 1910. 

13 Annales de Physique, 11, p. 253; 1919. 

4 Annales de Physique, 15, p. 246; 1921. 

1% J. Am. Chem. Soc., 48, p. 1274; 1926, 
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Figure 2.-—Comparison of specific heats 
0 Bose. A Mills and MacRea. 
X Kelley. ©} Dejardin. 


@ Williams and Daniels. V Trébin. 
@ Forch. © Andrews, Lynn, and Johnston. 








898 Bureau of Standards Journal of Research 


2. LATENT HEAT 


In the following comparison of the various values of latent heats 
which appear in the literature, those data based on direct calor. 
metric measurements will be given first consideration. The value: 
which have been calculated from other types of measurements throug) 
the Clapeyron relation are not considered of comparable precision. 

An attempt has been made to reduce the previously published 
values of latent heats to present international joules for purposes oj 
comparison. Such a reduction involves considerable uncertainty, 
but this is believed to be within the limit of error involved in the meas. 
urements themselves. It is not possible to make reliable estimates of 
the purity of the liquids used by various investigators because of the 
lack of published details. 
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Figure 3.—Comparison of latent heats of benzene 
O Brown. @ Nagornow and Rotinjanz. 


V Griffiths and Marshall. @ Mathews. 
— Tyrer. © Sutcliffe, Lay and Prichard. 


The data on methyl alcohol and ethyl alcohol are too meager to 
warrant exhibition on deviation charts. Figure 3 gives a comparison 
of the experimental values of the latent heat of benzene. 

The literature records two determinations of the latent heat of 
methyl alcohol at its boiling point. The first of these is by Brown " 
and the second by Mathews.’ Brown’s value is reported in calories 
and has been multiplied by 4.199 and decreased by 1 part in 512 to 
convert to present international joules. When so transformed, his 
value for L at 64.7° C. becomes 1,101.6 international joules per gram 
Mathews’ value for LZ at the same temperature reduces to 1,100.4 
international joules per gram. The latter is identical with that 
indicated by the present experiments. For methyl alcohol the agree- 
ment among the few ‘results 2 is satisfactory. 


6 J, Chem. Soc., 83, p. 987; 1903. 17 J, Am. Chem. Soc., 48, p. 562; 1926. 
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Three values for the latent heat of ethy! alcohol at its boiling point 
are recorded. Marshall and Ramsay determined the ratio of the 
latent heat of absolute ethyl alcohol at its boiling point to that of 
benzene at the same temperature. This ratio, 2.293, together with 
the value of Griffiths and Marshall for the latter, indicates a latent 
heat for the alcohol of 902.7 international joules per gram at 80° C. 

A determination by Brown,’ when transformed as already described 
vives a value of L for ethyl ‘alcohol at 80° C. of 903.5 international 


ol joules per gram. ‘The work of Mathews” gives the value 840.3 
ly, international joules per gram for LZ at the same temperature. The 
as- density which he records for his sample indicates that it may have 


contained nearly 0.2 per cent of water, The comparative value 
indicated by the All experiments is 837.3 international joules per 
eram. The alcohol samples which gave values of the latent heat of 
over 900 joules per gram were probably not dry. The agreement with 
the value of Mathews is moderately good, the difference probably 
being largely due to the small amount of water in his sample. 

Considerable data are available on the heat of vaporization of 
benzene. These will be discussed in chronological order. 

The work of Griffiths and Marshall #4 at temperatures from 20° to 
50° C. gave values which were represented by the equation L= 
107.05 -0.158 9, in terms of a calorie defined as 4.199 joules. To 
convert to present international joules, the coefficients of the right 
hand member have been multiplied by 4.199 and decreased by 1 
part in 512, because of the changes in electrical units since the work 
was performed. The transformed equation which gives values of 
L in international joules is L = 448.6 — 0.66 9. 

The work of Tyrer” gives a value of the latent heat of benzene 
at its boiling point in terms of a calorie based on indications of a 
voltmeter and silver voltameter. From the date of this work it is 
assumed that the units were those established by international agree- 
ment in 1910. On this basis the electrical equivalent of the calorie 
used by Tyrer becomes 4.188 international joules. His value for the 
latent heat of benzene at 80° C. then becomes 394.8 international 
joules per gram. 

Nagornow and Rotinjanz * report a group of determinations which 
yielded a mean value of 393.4 international joules per gram as the 
latent heat of benzene at 80° C. 

From the work of Mathews * the value 395.1 international joules 
per gram at 80° C. is indicated. 

Sutcliffe, Lay and Prichard * determined the latent heat of benzene 
from its boiling point up to 152° C., and by combining their data 
with those of Griffiths and Marshall, fitted an equation to the experi- 
mental values in the range 20° to 150° C. Transformed to yield 
values of Z in international joules per gram, this equation becomes 


L =447.72—0.61630— 6.141 X 10746? — 1.509 X 10-6 


The agreement among the values of LZ for benzene in the entire 
range is very gratifying. The equation of Sutcliffe, Lay, and Prich- 
ard, which represents the data of Griffiths and Marshall as well as 





18 Phil. Mag., 41, p. 38; 1896. 2J,C ran Soe. | 99, p. 1633; 1911. 


1° J, Chem. Soc., 87, p. 265; 1905. 23 Z. phys. C hem. : + a 700; 1911. 
*0 See footnote 17, p. 898. 24 See footnote 17, p. 


21 Phil. Mag. 41, p. 1; 1896, 25 Proc. Roy. Soc. Pow 115, p. 88; 1927. 
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their own, yields values which are everywhere in agreement with 
those calculated from our equation, within the limit of our owy 
experimental error. 

For all three of the fluids the values of L derived from the present 
measurements differ from those caluclated from the data of Young 
with the aid of the Clapeyron relation by approximately 1 part in 
50 at the boiling points and 1 part in 100 at 110° C. 


VIII. SUMMARY 


Thermal properties of methyl alcohol, ethyl alcohol, and benzene 
were investigated in the temperature range 40° to 110° C. Deter- 
minations were made of the latent heats of vaporization and the 
changes in heat content of the saturated fluids in this range. The 
results have been assembled in the form of tables and empirical 
equations have been fitted to the experimental values to aid in the 
application of the data. A working table of the thermal properties 
of these three fluids at saturation conditions has been prepared from 
the present and previously published data. 
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ABSTRACT 





This paper gives the results of an investigation on the characteristics of air- 
plane rece iving antennas to determine whether an antenna arrangement could 
be devised which would have all the desirable electrical properties of the vertical 
pole antenna and yet be free from the mechanical difficulties encountered in the 
use of the pole antenna. The antennas studied include the inclined antenna 
with both forward and backward inclination, the horizontal dipole antenna, the 
horizontal L. antenna, the horizontal V antenna, the inclined V antenna, the sym- 
metrical transverse T antenna, and the symmetrical longitudinal T antenna. A 
theoretical treatment is given which enables the voltage induced by a radio range 
beacon transmitting station to be calculated for any receiving antenna in space. 
This theoretical analysis is used to determine the received voltage, course error, 
and localizing effect for each of the antenna types studied. An experimental 
study was also made to check the theoretical analysis. The results obtained by 
experiment check very well with the theoretical predictions for each type of 
antenna. The symmetrical transverse T antenna and the symmetrical fonal- 
tudinal T antenna, with vertical lead-in portions, are both found to fulfill the 
desired requirements. Neither of these antennas shows any course errors and 
give the same received voltage as the vertical pole antenna having much greater 
actual height, thus reducing the mechanical troubles caused by vibration and 
ice formation. An appendix gives the mathematical derivation of the equation 
used as a starting point for the theoretical analysis. 
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I. INTRODUCTION 


Since the advent of the radio range beacon, the vertical pole antenng 
has been widely used because of its nondirectional properties and 
consequent freedom from course errors in radio range-beacon recep. 
tion. The vertical pole antenna, however, is subject to considerable 
mechanical trouble because of vibration and ice formation, and {o; 
this reason, an antenna free from these mechanical difficulties and 
yet having all the desirable electrical properties of the simple po! 
antenna would be preferable. 

Accordingly, a study of various types of airplane antennas was 
begun, a theoretical analysis being checked by experimental observs- 
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FicgureE 1.—Experimental installations of antennas studied 
a, Inclined antenna; 6, horizontal dipole antenna (half of this was used for a horizontal L antenna); ¢,t 


view, V antenna; d, side view, horizontal V antenna; e, side view, inclined V antenna; /, symmetri: 
transverse T antenna; g, symmetrical longitudinal T antenna 


tions in the air and on the ground. <A number of antenna arrange- 
ments were includedjin this study. For each type of antenna, th 
tests in the air included observation of the received voltage, the !ocal- 
izing effect or variation of the received voltage in the immediate 
vicinity of the beacon tower, and the course error as observed bj 
circling the beacon. These were compared directly with the results 
obtained using the vertical pole antenna. The types studied included 
the inclined antenna, with both forward and backward inclination 
(one example of the latter being the trailing wire antenna), the hon- 
zontal dipole antenna, the horizontal L antenna, the horizontal V 
antenna, the inclined V antenna, the symmetrical transverse T at- 
tenna, and the symmetrical longitudinal T antenna. Figure 1 shows 
experimental installations of these antennas on an airplane. The 
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> results of this experimental work check the theoretical analysis with 
© cufficient accuracy for the desired purpose, although exact quantita- 
tive data are almost impossible to obtain in the air because of the 
difficulties incurred in any attempt to measure the geometrical! 
quantities involved. 


II. THEORETICAL ANALYSIS 


It is found convenient to treat the subject by the use of vector 
analysis. A somewhat similar treatment, using a trigonometric 
method, has been given by Capt. William H. Murphy, of the Signal 
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Ficgure 2.—Diagram showing quantities used in the theoretical analysis 


Corps.’ In this work, as in Captain Murphy’s, relative rather than 
absolute field intensities and received voltages are desired. Conse- 
quently, all constants are considered equal to unity and omitted from 
the vector equations. The airplane is also considered to be at unit 
distance from the transmitter. “ 

The transmitting loop antenna is taken in the 7-k plane, with its 
center at o (fig. 2),so that the 7-axis coincides with the axis of the 





''W. H. Murphy, “Space Characteristics of Antennae,” J. Frank. Inst., 201, pp. 411-429; April, 1926; 
203, pp. 289-312; February, 1927. 
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loop antenna. The airplane is located at the point Q whose position 
vector is 7, and is flying in the direction indicated by the line marked 
“Line of flight.” The line of flight is assumed to lie in a plane paralle| 
to the 2-7 plane. The vector d is a unit vector in the direction of the 
antenna, pointing away from Q. a, is the angle between the axis of 
the loop antenna (i-axis) and the projection of the line of sight (or 
the vector 7) on the om) plane. a, is the angle between the line of 
flight and the projection of the line of sight on the horizontal plane 
containing the line of flight. 6; is the angle of elevation of the air. 
plane. 8 is the angle of inclination of the antenna with respect to 
the horizontal plane. These angles are taken as positive when meas- 
ured in the directions indicated in the figure. he antenna on the 
airplane is here considered to be in the vertical plane containing the 
longitudinal axis of the fuselage. This limitation may be removed, 
however, and the antenna considered in any position, by replacing 
@ in the following equations by a:.+ a’, where a’, is the angle be. 
tween the longitudinal axis of the airplane and the horizontal pro- 
jection of the antenna. 

The effective value of the electric field at Q, at unit distance from 
the transmitter, may be written, using Gibb’s notation 


EB, =? Xi (1) 


The derivation of this equation is given in the Appendix. The effec- 
tive voltage induced in the receiving antenna by this field vector is 


E, = B,-4=7 xXi-a=[F (2 


If the vectors 7 and 4 are evaluated in terms of their scalar compo- 
nents 7), 72, 73, @), @2, and a3, the value of this expression may be con- 
puted. However, equation (2) is readily simplified still further: 


Yr, T2173 
E,=[Fia)=|1 0 O|= 
1,2 Az 


revo). 
“7 = AeP'3 — T2Ag 
203 
‘* ° . . . . 
Similarly, for a second loop antenna in the 7-F plane, with its axis 
. . . . “ . 
coinciding with the 7-axis 


B,=) XT 
and 
E,=7 X7-a= (ral 


010 
E,=\r1 T21s 
1, A2Az 


= 73 — TA, 


a — [firs 


Now the components 7, 72, 73, @1, 42, and a; of the unit vectors? 
and d may be expressed in any system of coordinates, but the most 
convenient form is that of the spherical coordinates shown in Figure 
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9 In this system, and since the vectors are unit vectors (that is, 


(r|=|a|= 1) 
1 = COS a COS B, 
T2=SIN a Cos B; (7) 
r3=s1n B, 

and 
aq, = cos (ay, + Qe) cos Bo 
d_= sin (a; + a2) COs Bo (8) 
a3;=sin By 


On substitution of these values, equations (3) and (6) become 
E, = dors — Tez = SIN (a, + ag) Sin B; COS Bo— SiN a, cos B, sin B, (9) 
BE, = yr3— 7; 43 = COS (a; + ag) sin B; COS Bo— COS a; COS B, sin B, (10) 


These two equations form the basis of the remainder of the theo- 
retical analysis. 

A course is obtained whenever F,=£,. Thus, for the general an- 
tenna shown in Figure 2, a course is obtained when 


sin (a,;+ a2) sin B; cos B2—sin a; cos 6, sin Bz (11) 
=CoOs (a,+ a2) sin B; cos B2—COs a, COs B; SiN Be 
or when 
tan 6, (cos. a@2— SIN a)— tan Bp 
tan B, (sin a.+cos a2)—tan Bp 





(12) 


tan a,= 


where a is the course angle measured from axis of loop antenna 1. 

If the currents in the two loop antennas are equal, the true courses 
bisect the angles between the transmitting loop antennas. Assuming 
equal currents, then, the course error for any type of receiving an- 
tenna may be found. If the course error is here defined by e,? where 


then 
tan a — z! 
= ou I SS awe ; 
tan e=tan (a,— 7/4)= +i a (14) 
Combining (12) and (14) 
bead SIN a tan Bi (15) 


tan B.— cos a, tan fp, 


Mathematically, it is possible to make this course error zero in 
several ways: 6,=0, a2=0, or B= 5° The last is the only practicable 


method (heretofore accomplished through the use of a vertical an- 
tenna). However, a consideration of the physical significance of 
equation (11) shows that the course error is introduced “by the hori- 
zontal component of the field vec wr for equation (11) can be true for 
a =45° (that is, e=0) only if a,=0. Since a, can be considered as a 
generalized direction angle of the Tisie sum of the horizontal projec- 
tions of the antenna, it follows that a, can be zero for all directions of 





‘7 his j is the usual way of defining course error in connection with flight on a radio range ceete. 
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flight only if the vector sum of the horizontal projections of the ap. 
tenna is zero. Obviously, the horizontal projection of the antenna i 
affected only by the horizontal component of the field vector, so thg 
it is evident that the course error is introduced by this horizont,! 
component. The horizontal component of the field vector can not le 
eliminated (except by changes in the transmitting antenna, which ap 
not convenient), so that it becomes necessary to construct the receiy. 
ing antenna so that it will not respond to this horizontal component: 
that is, as noted above, the vector sum of the horizontal projections o{ 
the receiving antenna must be zero. Since the use of flat top portions 
is thus permissible so long as they fulfill these conditions, it becomes 
possible to design antennas having smaller actual height than that of 
the usual pole antenna, and yet having the same effective height, 
Such an antenna would give the same received voltage as the vertical 
pole antenna, and would likewise be immune to course errors, 


III. APPLICATION OF THEORY TO SPECIFIC ANTENNA 
STRUCTURES 


The theory is readily applied to indicate the results obtainable from 
several types of simple antennas. In each case, it is desirable to 
determine the localizing effect (the behavior of the signal in the in- 
mediate vicinity of the beacon), the course error, and the relative 
received voltage. The accompanying graphs (see fig. 3) show the 
localizing effect for the antennas considered. 


1. VERTICAL ANTENNA 


: Tr ‘ 
y= 9 (16) 


— 


(a) LOCALIZING EFFECT 


Substituting this condition in equations (9) and (10) and adding, the 
total received voltage is 


E= E, r E, _ — sin Q@, COS B, —COS a COS By (17) 


This equation is strictly correct only when the two loop antennas ar 
energized simultaneously and in the same time phase (for example, in 
the case of the two-course visual type beacon). However, the results 
secured in the following analysis are independent of this limitation, 
although the equations used in securing these results are somewhat 
different if the loop antenna currents are displaced from each other 
in time phase, or if the loop antennas are energized separately as 0 
the case of the aural type beacon. 

The value of / in equation (17) becomes zero when 


Thus, this antenna gives a zero signal zone directly over the beacon. 
(Fig. 3(A).) 
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DISTANCE FROM BEACON, ARBITRARY UNITS 





Figure 3.—Graphs of localizing effecis for antennas studied 


A, vertical antenna; B, inclined antenna (8;= 30°); C, inclined antenna (62.=60°); D, trailing wire antenna 
(8:=—30°); E, horizontal L. antenna; F, horizontal V antenna (angle between each wire and fuselage 
pepe Banal Bae inclined V antenna, angle between each wire and fuselage center line=30°, angle 
of inclination = 45°, 


49527°—31——10 
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(b) COURSE ERROR 


Substituting B= 5 in (15) 
e=0 (19) 
(c) RECEIVED VOLTAGE 


In equation (17), a, may have any value, but is constant for any 
particular line of flight over the beacon, so that EH is proportional ty 
cos 6;. Thus, at great distances, where 6, is small, the vertical ap. 
tenna gives a relatively large received voltage. 


2. INCLINED ANTENNA 


The antenna is considered in the same vertical plane as the fuselage, 
either inclined forward (8,>0), or trailing wire (6.<0). The general 
equations given above apply to this case. 


(a) LOCALIZING EFFECT 


This antenna gives both a zero signal zone and a maximum signal 
zone, one occurring before the beacon is reached and the other after 
it has been passed. (See equations (9) and (10).) If 8, is positive, 
the zero signal zone occurs first (fig. 3(B) and (C)), while if B, is negs- 
tive, the maximum occurs first. (Fig. 3(D).) The magnitude of 
8, determines the distances of these two zones from the beacon. 
The sharpness of the zero and the magnitude of the maximum are 
also dependent on the value of 82. As 6, is decreased, the zero signal 
zone recedes from the beacon, and the maximum increases in value 
and approaches the beacon. In the limiting case (8,=0), the maxi- 
mum occurs directly over the beacon and the zero signal zone disap- 
pears (recedes to infinity). 


(b) COURSE ERROR 


This is given by equation (15). Thus, if the airplane is flying 
directly toward or away from the beacon (a,=0) 


e=0 (20) 


If the airplane is circling the beacon ( that is, a~5), 


tan B, 


é = 
tan tan Bs 


If ag=7/4 and B,=7/4 
V2 tan By 


2— V2 tan B; 





tan e= (2 


In equation (15), when tan 8;—cos a, tan 8; becomes very small, tan ¢ 


us 


becomes very large and e becomes nearly equal to 9 
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(c) RECEIVED VOLTAGE 


& If 8. is large (greater than =) the received voltage in an inclined 


(19) He antenna is very nearly the same as that in a vertical antenna of the 
Fame length. However, if B2 is small, the inclined antenna is mark- 
© diy inferior to the vertical antenna at great distances from the bea- 
© on, although there is a region, in the vicinity of the maximum signal 

ry Fone mentioned above, where the inclined antenna is superior to the 

ee ‘vertical antenna. This is a small region close to the beacon, and of 
al: Bittle value practically, except for localizing effects. 


3. HORIZONTAL DIPOLE ANTENNA 


This antenna is in the horizontal plane, and perpendicular to the 


age, J line of flight. The results apply also to an L-type antenna from the 
eral [wing tip to the fuselage. In this case, the antenna does not lie in 


the vertical plane containing the longitudinal axis of the fuselage, 
and consequently a, in the preceding equations must be replaced by 


) QT =* Also 
nal 2 
ter © — a= 0 
ve, INOW 
4 . ' Tv 
pi sin (a:+5) =COS ae 
iF mi (23) 
— vis ° 
on, cos ( a+ > )~ —sINn Q 
ire 2 
ral ' ee : ; eet 
ne By substitution of these values in equation (15), it is found that 
va tan €=cot a, (24) 
p- ; 
and 
T 
e=5~ M2 (25) 
g Thus the course error is the complement of the angle at which the 


airplane crosses the beacon course. 
_ (1) Localizing effect—With this antenna, a maximum signal zone 
) & occurs over the beacon. As mentioned before, this is the localizing 
effect obtained in the limiting case of the inclined antenna, when 
B.=0. (Fig. 3(£).) 
(2) Course error.—As noted immediately above, the course error is 
the complement of the angle at which the airplane crosses the beacon 
) course. 
(3) Received voltage—From equation (24), it is seen that the 
received voltage varies as sin 6, so that it becomes very small when 
the airplane is any appreciable distance from the beacon. 


4. HORIZONTAL V ANTENNA 


One variation of the horizontal antenna is the V antenna from wing 

tips to vertical fin. The performance of this type of antenna is 
essentially the same as that of the horizontal dipole, except that the 
received voltage is somewhat greater, since both horizontal compo- 
nents of the transmitted signal are received. The localizing effect is 
shown in Figure 3(F). 
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5. SYMMETRICAL T ANTENNA 


This can be either longitudinal or transverse. The received voltag: 
of this type of antenna is due entirely to the vertical portion, the fii; 
top merely increasing the charging current at the top of the verticg| 
portion, thus increasing the effective height. Consequently, the 
antenna will have the characteristics of its vertical or lead-in portion, 
If this is truly vertical, the arrangement will act the same as the simple 
vertical antenna, while if the lead-in portion is inclined, the char. 
acteristics will be those of the inclined antenna. The T antenn 
gives greater effective height than either type simple antenna having 
the same actual height because of the more uniform current distr. 
bution resulting from the additional capacity at the top of the 
antenna. 

6. OTHER ANTENNA SYSTEMS 


The general theory outlined in Section IT is also applicable to other 
and more complicated antennas. For example, consider the inclined 
V antenna from the wing tips to the vertical fin, with the lead-in 
running from the junction of the two wires at the vertical fin. This 
is essentially a form of T antenna, the effects of the horizontal com- 
ponents perpendicular to the axis of the fuselage cancelling. The 
antenna, therefore, performs as an inclined antenna having a vertical 
portion and also having a horizontal portion in the axis of the fuselage. 
The localizing effect for this type of antenna is shown in Figure 3((), 

Similarly, consider any antenna having a vertical lead-in and a 
iat-top arrangement symmetrically disposed about the vertical lead-in 
so that the vector sum of the horizontal projections of the individual 
elements of the flat top is equal to zero, This antenna will behave 
in every way as a vertical antenna, but giving greater received 
voltage than a vertical antenna of the same actual height. It is this 
type of antenna which offers advantages of reduced mechanical 
vibration over the conventional vertical pole antenna. When the 
flat-top portions lie along the longitudinal axis of the fuselage, for 
example, the longitudinal T antenna, the problems due to ice forma- 
tion are also reduced. 


IV. EXPERIMENTAL RESULTS 


Antennas of each of the types analyzed above were installed on 
the National Bureau of Standards’ airplane and compared (in flight 
on the visual radio range beacon) with the vertical pole antenna as a 
standard of comparison. Relative effective heights were determined 
by switching the antenna under test and the pole antenna alternately 
to the same receiving set and recording the receiving set output 
voltage. Course error effects were obtained by using two receiving 
sets, one connected to the antenna under test and the other to the 
pole antenna, and comparing the course indications (as obtained on 
two reed indicators connected in the output circuits of the two 
receiving sets), while circling the beacon at various altitudes and 
various distances from the beacon. The localizing effects were 
obtained by flying at a constant altitude on a beacon course directly 
over the beacon tower. 





Diam 
Davies 


T 
tati 
mee 
wer 








Vol ¢ 





Diamond | Receiving Antennas for Aircraft 911 


Davies 





The results secured in the air were approximate only, exact quanti- 
tative measurements being very difficult because of the difficulty in 
measuring the geometric quantities involved. However, the results 
were sufficiently accurate to corroborate the theoretical analysis. 









1. RELATIVE RECEIVED VOLTAGES 








The data obtained showed that, of the antennas studied, the 
symmetrical T antenna of either type and the wing tip to vertical 
fin inclined V antenna both gave effective heights equivalent to that 
obtained with the vertical pole antenna without requiring nearly 
the same actual height as the pole antenna. For example, a sym- 
metrical transverse T antenna having a 12-inch vertical lead-in and 
a flat top extending 15 feet on each side of the fuselage, parallel to 
and 12 inches above the wing surfaces, gave slightly better effective 
height than a 5-foot vertical pole. 








he 








2. COURSE ERROR EFFECTS 






(a) INCLINED ANTENNA (8 2>=20°, INCLINATION FORWARD) 


When circling the beacon at constant distance from the beacon 
tower, the on-course indications of the reed indicator operated from 
the receiving set connected to the inclined antenna were in advance 
of the on-course indications of the reed indicator in the output circuit 
of the receiving set connected to the pole antenna. The angle of lead 
increased with decreasing distance from the beacon and also increased 
with the altitude of the airplane for a given distance from the beacon. 
This angle of lead is equal to the course error. The results are in 












. ‘ , : Ts 4 ie 
accordance with equation (15), placing a2 = 3» since, when circling 





the beacon, the line of flight is perpendicular to the line of sight. 
Equation (15) then resolves into equation (21), namely 







where £, is the angle of elevation of the line connecting the airplane 
with the beacon (the line of sight) and £, the angle of inclination of the 
antenna (20° in this case). Some idea of the magnitude of the course 
errors involved may be obtained from one specific case. Flying at an 
elevation of 3,000 feet along a circle of three miles radius, the course 
error is 27.5°. 

When flying along a beacon course directly to or from the beacon, 
no course error was obtained. This is also in accordance with equa- 
tion (1 5), since, when a, =0, equation (15) resolves into equation (20), 
hamely 











tan e=0 






(b) INCLINED ANTENNA (£:=—20°, INCLINATION BACKWARD) 


The results obtained were exactly the same as for the antenna with 
forward inclination except that the on-course indications of the 
reed indicator operated from the receiving set connected to the in- 
clined antenna lagged the on-course indications obtained on the reed 
indicator in the output circuit of the set connected to the pole antenna. 
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(c) INCLINED v ANTENNA 


The results secured were the same as for the inclined antenna wit) 
forward inclination. 


(d) SYMMETRICAL T ANTENNA (EITHER TYPE) 


No course errors were obtained under any conditions. Both the 
transverse T and longitudinal T antennas were tested. 


3. LOCALIZING EFFECTS 


The localizing effect for each of the antennas studied was found to be 
very nearly the same as shown in the graphs of Figure 3. This was, 
perhaps, the most satisfactory test of all, checking, within the limits 
of error, the exact theoretical analysis. At 3,000 feet altitude, by 
changing the angle of inclination of the inclined antenna from 20° 
to 90°, the zero signal zone could be moved from approximately | 
mile from the beacon to directly over the beacon tower. When 
setting up the longitudinal T antenna, the zero signal effect directly 
over the beacon tower was employed as a check on the electrical sym- 
metry of the flat-top arrangement. 


V. CONCLUSIONS 


For reception of signals from radio range beacons, an antenna 
arrangement of which the symmetrical T antenna is a typical exampl 
was found to give good electrical performance and to be free from 
trouble due to vibration. In addition, the performance of this type 
of antenna under conditions favoring ice formation should prove 
superior to that of the vertical pole antenna. To reduce the effects 
of mechanical vibration, this antenna employs vertical poles consider- 
ably shorter than the conventional pole antenna (10 to 18 inches 
instead of 5 to 6 feet). Equivalent effective height is secured through 
the addition of a flat top. In order to prevent directional effects, 
which introduce course errors, the flat top is made up of two, three, 
or more elements symmetrically disposed about the vertical lead-in, 
so that the vector sum of the horizontal projections of the individual 
elements is equal to zero. When but two flat-top elements are 
employed, the antenna takes the form of either a transverse or longi- 
tudinal T. The latter arrangement is preferable from the viewpoint 
of freedom from trouble due to ice formation. Furthermore, such a 
longitudinal T antenna has an aerodynamical resistance very much 
smaller than that of the transverse T antenna, and probably as small, 
if not smaller, than that of the vertical pole antenna. 


VI. APPENDIX 


This calculation of the field due to a loop antenna is based on the 
Lorentz solution of the electromagnetic field equations for a known 
distribution of currents and charges. 

The notation used (Gibbs-‘‘ Vector analysis’’) is as follows: 

H = magnetic field intensity. 

E'= electric field intensity. 
p=density of charge. 
v=velocity of motion of charge. 
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(Voli Davies 
',,c, K = constants depending upon system of units used 
with vVxX= curl. 
vy: =divergence. 
= Ww : Paar 
B Pot w= ff fi dr where the integration is extended over all space, dr 
the represents the element of volume, and r the distance from 
dr to the point at which the potential is computed. 
],= Maximum value of current in loop antenna. 
w =angular velocity of current in loop antenna. 
we t= time 
oe V = velocity of electromagnetic waves in medium surrounding loop. 
its (u]=value of wu at time t— \ where 7 is the distance from source 
by 


of disturbance to point at which [uw] is considered. 
The loop is assumed to be square and the current to be the same 
; throughout the loop. 
In this notation, Maxwell’s equations of the electromagnetic field 


are 
4n/ - KOE 
~ € ( p+ 7 7) 















A solution of this system sometimes used (neglecting a term in the 
value of # which is immaterial in this discussion) may be written 





H- 


‘— 


VX Pot [p2] = ; Lap [pv] 


s 





(2) 





ley 


mu o — >. 
B= —4,9 Pot [orl= Pot] 34000) | 





2 





Now 





ov =1=1b=Iob sin wt 







where 6 is the unit vector in the direction of the current. 





B= —4pot| 3 | (4) 
) OF oleh cos wt (5) 





rt ae. 
E=— pul Pot (B cos wt] = — pool, b cos uf t - v) 
a - — ds 


oh 
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where h/ represents the distance from current element to P of Figuy 
4, which shows the position of the loop antenna and the distances an¢ 
vectors used in the following treatment: 
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Fiaure 4.—Diagram showing quantities used in appendiz 
f PP 


~ wail _cos w(t 7, ,) ath o- cos w( t 
ee . :; q 


_ cos w ( - -] 
+9 Sa if. ds 


But 
7 ds = 2a 
, : 
2 nwa] jcos w( t-) cos w( tf) 
Oe eee Qo nite . - ee 
el | 
| cos w(t cos w(t. v) 7 
| “ies en ee 

Now 


B=i—q=rd+ (aij tnd 
fer deg eae tape 
m=r—ak=ryit rg + (r,—a)k 





n=r+ak= rit rej + (r3+a)k 
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=(r+ (r— 





Also, by the same method 


=COS w (1 — pt 


1 
cos wo (¢- yt 


~ ~20n a (4-5) 
, 7) yp /sin 
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Thus 





+ Ley 
. YF 
1_1l_an 
¢ rf # 
S58 an 
mr P 
1_1_ ar 
mRrr 





cos(r+y)—cos(rx—y)= — 2 sin zsin y 


cos(z+y) +cos(4— y) =2 cos x cos y 


the coefficients of 7 and & in (9) may be simplified. 


cos w ( (=) _ cos a 


(Lt) -(1-38) a 








r — cos w (1 - 3) | 
15h u( (1-54 +co8'a ( -7- 3) | 
® Are =+( aoe 
—— +s" cos w (1 r) 


sin « (t~ 4) +252 cos o (t- F) 


a)? +732]'? = (r?—2r.a + a)? =r— 






are 


© This is an approximation obtained by expanding the value of p by 
> the binomial theorem and neglecting terms involving a in powers 
higher than the first. 





(13) 





By the use of equations (12) and (13) and the trigonometrical for- 


yee see ~*~ 
Ve rv, 


® A> 


cos Tv 












(11) 
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Re wt 2ra 
since V = zs ’ 










: 4 : : . war 
and @ is small in comparison with \, sin V nay 
7, Mal 


9 9 
2rare 2rar, 
be replaced by “<— and cos may be considered equal to one. 








Similarly 
n m 
COS w (1 TP. ieas w (¢ y) 
n m 
_4n ars: r 2ars r 
Vr sin o (1- — 2s cos o (t—F 
Then 


ome — ReFtin w (tf) +292 cos w( 1-7) | 
+3[ S22 sin wo (t-F — Ustcos w ( t— | 
= Hetle( , d—nk)| $5 sin o (¢— *)- cos w(t- 


raj — Tok =r Xi 


= [| 8xywa’l, . r\ 4 a x , 
 Ba[ Beets sin w (t-F)- - - os w(t-4) | 


The first term represents the radiation field and the second term 
the induction field. If the radiation term alone is considered, no 
error is introduced except one of magnitude for small values of 1. 
As interest in this treatment centers chiefly on relative rather than 


absolute field intensities, the quantity Szpoee may be considered 


Now 





constant and equal to unity. Furthermore, the airplane is consid- 
ered to be at unit distance from the transmitter, so that equation 
(17) may be written 


Ey. =rXt 


However, for the calculation of curves showing the localizing effect 
when an airplane flies at constant height over the beacon, it is neces 
sary to take the variation of r into account and the equation then used 
is 


re 
WasuINGTON, January 24, 1931. 


O 





